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We present the coefficient functions for deep-inelastic scattering (DIS) via the exchange of a scalar 
(|) directly coupling only to gluons, such as the Higgs boson in the limit of a very heavy top quark 
and n^ effectively massless light flavours, to the third order in perturbative QCD. The two-loop 
results are employed to construct the next-to-next-to-leading order physical evolution kernels for 
the system (F2,F§) of flavour- singlet structure functions. The practical relevance of these kernels 
as an alternative to MS factorization is bedevilled by artificial double logarithms at small values 
of the scaling variable x, where the large top-mass limit ceases to be appropriate. However, they 
show an only single-logarithmic enhancement at large x. Conjecturing that this feature persists to 
the next order also in the present singlet case, the three-loop coefficient functions facilitate exact 
predictions (backed up by their particular colour structure) of the double-logarithmic contributions 
to the fourth-order singlet splitting functions, i.e., of the terms (1 — x) a \n k (l — x) with k = 4, 5, 6 
and k = 3, 4, 5, respectively, for the off-diagonal and diagonal quantities to all powers a in (1— x). 



1 Introduction 



Structure functions in lepton-nucleon deep-inelastic scattering are classic observables probing 
Quantum Chromodynamics (QCD), the theory of the strong interaction, and in particular the 
structure of the nucleon. Indeed, structure function measurements provide the backbone of our 
knowledge of the quark and gluon longitudinal momentum distributions in the proton to this day, 
40 years after the pioneering measurements of DIS at SLAC [1], see also Refs. [2]. In turn these 
distributions, in short referred to as the parton distributions or parton densities of the nucleon, are 
indispensable for the description and analysis of hard scattering processes at the proton colliders 
forming the high-energy frontier of particle physics, i.e., the TEVATRON and the LHC. See, e.g., 
Ref. [3] for introductory overviews and the present status of structure functions and QCD. 

The asymptotic freedom of QCD [4] facilitates a partly perturbative description of hard hadron 
processes. At all practically relevant scales, however, the strong coupling constant a s is much 
larger than its electroweak counterparts. Hence calculations to a higher order are required in per- 
turbative QCD in order in arrive at precise predictions. Also in this respect structure functions 
form a special tool, since they can be expressed in terms of the operator-product expansion [5] 
which provides a framework that renders higher-order calculations considerably more accessible. 

In a series of previous articles, three of us have exploited this fact to derive the exact third- 
order QCD corrections to the splitting functions governing the scale dependence of the parton 
distributions [6, 7] - and thus the transfer of information from DIS at moderate scales to the high 
scales involved in Higgs, gauge boson, top quark and new physics processes at Tevatron and the 
LHC - and to the partonic cross sections (coefficient functions) for the most important structure 
functions [8-10]. As all modern higher-order calculations in perturbative QCD, those calculations 
have been performed in dimensional regularization [11], and the splitting functions and coefficient 
functions have been determined in the modified [12] minimal subtraction [13] scheme, MS. 

Consequently the higher-order parton densities are not defined in terms of any physical process, 
but by a mathematical description employed for isolating the mass singularities which lead to 
their scale dependence. The MS scheme has definite advantages, e.g., it leads to a stable (order- 
independent) functional form of the dominant diagonal (quark-quark and gluon-gluon) splitting 
functions in the limit of large momentum fractions x [14], see also Ref. [15], a feature that assists 
a stable evolution of the parton densities over the wide range of scales mentioned above. 

An alternative approach has been suggested already long ago in Ref. [16], see also Refs. [17]: 
As discussed below, the parton densities can be eliminated from the description of the structure 
functions (and other hard processes). The dependence of the observables on the physical hard 
scale (in DIS Q 2 = —q 2 , where q is the four-momentum of the exchanged electroweak gauge 
boson) is then given in terms of so-called physical evolution kernels or physical anomalous dimen- 
sions. While the direct relation between different observables via the universal parton densities is 
obscured in this approach, it can be practically useful for analyses of selected observables such as 
the determination of 0C S from the scaling violations of structure functions [18]. 
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The construction of physical evolution kernels is particularly simple (and unique) for flavour 
non-singlet, i.e., gluon insensitive observables. It turns out that these kernels show an only single- 
logarithmic enhancement at large values of the scaling variable x for a wide range of non-singlet 
quantities in DIS, semi-inclusive e + e~ annihilation (SIA) and the Drell-Yan process [19,20]. This 
general behaviour is in contrast to that of all corresponding MS coefficient functions which receive 
double-logarithmic contributions. As a result the physical kernels provide an - at present not 
yet formally proven - all-order exponentiation of the highest In ( 1 — x) terms beyond the leading 
( 1 —x) ~ contributions covered by the standard threshold resummation [2 1 ] to all powers in ( 1 —x) . 

In the present article we address the flavour singlet case which requires two observables and a 
2x2 matrix evolution kernel. Unlike Refs. [17, 18] we study the ideal - for theoretical purposes - 
gluonic complement to the most important gauge-boson exchange structure function F2, i.e., DIS 
by the exchange of a scalar § coupling directly only to gluons via an additional term G^G^ in the 
Lagrangian, where G*™ represents the gluon field strength tensor. Such a term, already suggested 
as a trick in Ref. [16], is of course not present in the fundamental Standard Model Lagrangian. 
However, it does occur effectively for the Higgs boson in the limit of a very heavy top quark [22]. 

This interaction was included, in order to directly access also the quark-gluon and gluon-gluon 
splitting functions, in the calculations of Ref. [7] and, in fact, its precursors [23,24]. Hence the one- 
and two-loop coefficients functions for the resulting structure function were calculated, but not 
published, quite a while ago. The more recent, but also unpublished three-loop results have already 
been employed to determine the third-order gluon jet function in threshold resummation [25] and 
the pole terms of the three-loop gluon form factor [26] - and from the latter the next-to-next- 
to-next-to-leading logarithmic (N 3 LL) threshold resummation for Higgs production in the heavy- 
top limit [27]. These coefficient functions also provide, via an analytic continuation to timelike 
kinematics, a strong check of the next-to-next-to-next-to-leading order (N 3 LO) computation in 
Ref. [28] of the top-induced Higgs decay into hadrons [29]. The first- and second-order coefficient 
functions for F§ have now been derived independently in a completely different manner [30]. 

The outline of this article is as follows. In Section 2 we introduce the formalism used for 
calculating the coefficient functions for F§ to order oc s 3 , briefly discuss the diagram calculation and 
its checks, and provide the transformation to the physical kernels to the fourth order. In Section 3 
we then present the exact expressions for the coefficient functions and address their size and end- 
point behaviour. The one- and two-loop results are then combined in Section 4 with the known 
second- and third-order splitting functions to derive the next-to-leading order (NLO) and next-to- 
next-to-leading order (NNLO) physical kernels for the system (F2,F^). Also these kernels show an 
only single-logarithmic large-* enhancement to all powers in (1— x), a feature that we conjecture 
to hold also at the next order. In Section 5 this conjecture and the three-loop coefficient functions 
of Section 3 are then employed to predict the coefficients of the highest three powers of In ( 1 — x) in 
all four singlet splitting functions at order (N 3 LO). We summarize our results and give a brief 
outlook in Section 6. Appendices A and B contain the partly very long expressions, respectively, 
for the new third-order coefficient functions and the elements of the NNLO physical kernel matrix. 
In Appendix C we finally provide the large-* limit of the gluon coefficient functions for F2. 
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2 Formalism and calculation 



We consider inclusive deep-inelastic lepton-nucleon scattering through the exchange of a scalar 
(Higgs) boson (|), which proceeds through the reaction 

<Ktf) +nucl(p) -> X . (2.1) 

Here X denotes all hadronic final states allowed by quantum number conservation. The boson § 
transfers a space-like momentum q (i.e., Q 2 = —q 2 > 0), while the nucleon carries momentum p. 
The scaling variable of the reaction is defined as usual as x = Q 2 / (2p ■ q) with < x < 1 . 

In complete analogy to ordinary (gauge-boson exchange) DIS, the cross section for this process 
can be parametrized in terms of a structure function Through the optical theorem the total cross 
section (hence the structure function Fa) is related to the imaginary part of the forward amplitude 
Ta for the scattering process in Eq. (12.11) of a virtual Higgs boson off the nucleon, 

fy(p,q) = if d 4 ze iqz (nuclei T(J(z)J(0)) |nucl,p) . (2.2) 

This quantity contains the time-ordered product of two scalar currents (see below) to which the 
standard operator-product expansion (OPE) can be applied in the Bjorken limit (large Q 2 for 
fixed x). The relevant steps have been already discussed in previous lower-order and fixed Mellin-/V 
third-order calculations, see Refs. [23, 24], hence we can be brief in recalling some key issues. 

The OPE decomposes the current product in Eq. (12.21) in terms of the standard set of spin- 
averaged matrix elements of the (renormalized) spin-/V twist-two irreducible flavour- singlet quark 
and gluon operators 

(nucl )jP | O^'-'^lmicl^) = p^...p^A ijriucl (N,p 2 ) (2.3) 

and the respective hard scattering parton coefficient functions C§j with i = q, g, where p stands 
for the renormalization scale. The operators O q and O g in Eq. (12.31) arise from the symmetric and 
traceless part of N covariant derivatives D^ 1 , respectively acting on the quark (\|/) and gluon fields 
(G MV ). Neglecting l/Q 2 power corrections, the OPE applied to Eq. (12.21) allows us to express the 
Mellin-./V moments of the structure function F§ in terms the matrix elements A !;nuc i and coefficient 
functions Ca * as 

= f Q dxx N - l F^x,Q 2 ), (2.4) 
= E C M (^V,^,a s ) A,, nucl (^,A/ 2 ) . (2.5) 

!=q,g V A* / 

Here the first line (12.41) fixes our conventions for Mellin moments and the second relation (12.51) 
holds (due to symmetry properties for x — > —x) for even-integer values of N. All (complex) mo- 
ments, and the complete x-dependence, are uniquely fixed by analytic continuation though. 
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The renormalization of the singlet operators O q and O g in Eq. (12.31) in terms of their bare 
counterparts proceeds as 

Oi = Z lk O h k ™, (2.6) 

where the renormalization factors Z^ are matrix- valued (i.e., summation of k is understood) and 
the dependence on N has been suppressed here (and in the following) for brevity. The well-known 
anomalous dimensions ji k arise from the renormalization factors as 

1f * = -(^ z «)( z " 1 ^' <17) 

and the scale dependence of the operators Oi is given by 

d 



Oi = -J ik O k = P lk O k . (2.8) 



Jin /j 2 

Eq. (|2.8I) also recalls the conventional relation between the anomalous dimensions and the mo- 
ments of the splitting functions Pi k {x). 

Using dimensional regularization [11] in D = 4 — 2e dimensions and the modified minimal 
subtraction scheme MS [12, 13], the renormalization factors Z ik in Eq. (12.61 ) are a given by a series 
of poles in 1/e which can be expressed in terms of the perturbative expansion coefficients yW of 
the anomalous dimensions and the coefficients |3 n governing the running coupling. In a power 
expansion in oc s the former can be written as 

oo 

y ik = L4 +1 y^, (2-9) 

1=0 

while the latter are given by 

rfiy^ " Jhv = — • • (2 - 10) 

Here Pi ... 4 denote the known four-dimensional expansion coefficients of the beta function of QCD 
[4,31-33], Po = H/3Ca — 2/3«y etc, with representing the number of active quark flavours 
and Ca = 3 (and Cp = 4/3) denoting the usual SU(n c = 3) colour factors. 

The N'LO anomalous dimensions yO are related to the 1/e single poles of Z/fc in Eq. (12.61) at 
order a s /+1 . The D-dimensional coefficient functions C^/, on the other hand, have an expansion in 
non-negative powers of e, viz 

if JO , „„(0 , ,2,(1) 



C*,i = S ig + E 4 {cf} + eo$ +e 2 ^ + . . . (2.11) 

where z = q, g as before, and we have again suppressed the dependence on N (and Q 2 //j 2 ). Recall 
that contributions of order £ k>0 in Eq. (12.1 II) enter the extraction of the anomalous dimensions 
and coefficient functions at higher orders in a s . For example, as shown in Eq. (12.191 ) below, the 
determination of the third-order coefficient functions c^] requires the two-loop e 1 quantities a^j. 
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For the perturbative determination of splitting functions and coefficient functions the nucleons 
in Eqs. (12.21) and (12.31) can be replaced by partons. Thus we need to compute all Feynman diagrams 
for the forward scattering amplitudes which contribute at leading-twist accuracy to the singlet 
structure functions i 7 ^ in partonic DIS, 

T^ k (p,q) =ifd A ze^{k,p\T{J{z)J{Q))%p), * = q, g , (2-12) 

i.e., we are considering the reactions, 

parton (p) + scalar (q) — > parton (p) + scalar (q) (2.13) 

where the scalar couples directly only to gluons via a §G MV G MV contribution to the Lagrangian. 
Such an interaction is effectively included in the Standard Model for the Higgs boson, in the limit 
of a heavy top quark and negligible Yukawa couplings to all other quark flavours [22]. Actually, 
the resulting coefficient functions for the Higgs boson differ from those of a scalar with a generic 
tyG^G^v coupling by a perturbative prefactor, known to N 3 LO [34], which is however irrelevant 
for our considerations. We will return to the limitations of the heavy-top limit in Sections 3 and 4. 

The above scalar current in Eq. (12.121) requires an additional, well-known renormalization [35], 

( G -g, v ) = Zc2 ( G ^r + - ■ %• = T^wm^) • (214) 

which is calculable in terms of the coefficients [3 n of the QCD beta function (12. 101 ). Moreover the 
field strength in Eq. (12.141) is subject to operator mixings which, however, give either vanishing 
contributions to the on-shell matrix elements considered in Eq. (12.131 ) or, as in case of a quark 
mass term m q \j/\|/, vanish in the present limit of massless quarks. 

In order to determine the coefficient functions C^., to order aj, we thus 'only' need to evaluate 
the forward scattering amplitudes (12.131) to three loops. The corresponding Feynman diagrams 
have been generated automatically with Qgraf [36]. The resulting number of diagrams is shown 
in Tabled] Note that, as already in Refs. [23], the generation and counting of the diagrams are 
non-standard, as some tricks (e.g., using symmetries valid only at the relevant even values of N) 
have been employed to reduce the number of diagrams. 



process 


tree 


1-loop 


2-loop 


3-loop 


qty -> q§ 




1 


23 


696 


g§ -> gty 


1 


8 


218 


6378 


h§ — > h§ 




1 


33 


1184 


sum 


1 


10 


274 


8258 



Table 1 : The number of diagrams for the forward amplitudes employed in the calculation of the 
structure function F§ to three loops, h stands for the external ghost discussed in the text. 
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For all external partons states in Eq. (12.131) we need to project on their physical polarizations. 
To that end it has proven efficient to contract the external gluon lines only with the metric tensor 
_g«P instead of a physical projector 

p-q [p-qY 

matching the kinematics of Eq. (12.131) with p 2 = 0, q 2 = —Q 2 < 0. To compensate the additional 
unphysical degrees of freedom we have to consider an extra class of diagrams with external ghosts 
instead of external gluons, i.e., the process h§ ^ h§ listed in the third line of Table [T] In particular 
the resulting absence of higher tensor integrals from the numerator of Eq. (12.151) leads to a vital 
simplification of the calculations. For the same reason our all-TV computations have been performed 
in the Feynman gauge. As before the gauge independence has been verified for a few low values 
of N along the lines of Refs. [23] where also the projector (12.151) has been employed for checks. 

The calculation of the Mellin-/V projection of the forward- scattering integrals for all values 
of N has been discussed extensively, if not exhaustively, in previous publications [6,7,9,24,37]. 
Also the present calculation has been carried out using version 3 of the symbolic manipulation 
program Form [38] to which new features, such as an efficient database facility for huge numbers 
of intermediate integrals, were added in order to make such computations feasible [39]. We would 
like to stress again that the calculation has been set up in such a manner that checks of fixed low 
moments could be performed at all intermediate stages using the program Mincer [40,41] for 
massless three-loop self-energy integrals. Note also that the recent computations of fixed moments 
of heavy-quark structure functions [42] and of the three-loop quark and gluon form factors [43] 
provide strong independent checks of different aspects of our calculations. 

Having calculated the perturbative corrections to the individual parton contributions to 7^ in 
Eq. (12.121) it remains (first) to perform the renormalization of a s and the scalar current in Eq. (12.141) . 
and (second) to disentangle the coefficients of the Laurent expansion in terms of anomalous dimen- 
sions and coefficient functions. The second task also known as mass factorization is non-trivial, 
requires an explicit representation for the matrix Z% in Eq. (12.61) and introduces the dependence on 
a specific scheme, i.e., MS in our case. To that end, we set the factor exp(£{ln(47t) — y e }) = 1 where 
y e is the Euler-Mascheroni constant and choose, without loss of generality, the renormalization and 
mass factorization scales as /j 2 = Q 2 . 

A suitable normalization, already used in Eq. (12.111) above, at order a® implies 

r (0) _ c (0) _ r (0) _ (0) _ 2 6 

At the first order in oc s the respective forward amplitudes need to be calculated up to order £ 2 for 
our purposes, yielding 

T W - I Y (°) + c « + ea W , £ 2> (i) {1 1T) 

Vp ~~ £ ^ gP <t»,P 0,P ' ^ ' 

with p = q, g. Correspondingly the two-loop contributions are required up to order £. These 
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quantities are given by 

where 8i p is the Kronecker symbol and summation over repeated indices is understood. Finally we 
are ready to write down the third-order coefficients in Eq. (12.121) . 

1% = ^{vL° ) vL 0) v£ , -3PoiMp 0, +2Po 2 Yg ) } 

+ ^ + ia<: i )y«+^Y|p 0, + 5^ (yL '- W^Yg) , (2.19) 
from which the N LO coefficient functions ci ' are extracted. 

As mentioned above, the mass factorization in Eqs. (12.171) - (12.191 ) provides the splitting func- 
tions and coefficient functions in the MS scheme at the scale /j 2 = Q 2 . The full scale dependence, 
including the case of unequal renormalization and factorization scales, can be reconstructed from 
these results using, e.g., Eqs. (2.16) - (2.18) of Ref. [44]. 

For transformations to other schemes it is convenient to combine the present 'gluon' structure 
function Fa, recall Eq. (12.161 ), with a corresponding 'quark' observable such as (the flavour-singlet 

part of) F2 in photon-exchange DIS [9], with = 1 and = 0. Hence we define the two- 
dimensional vector of the moments of the singlet quark and gluon distributions, 



q = with q s = £fe + <7/) 



(2.20) 



and the corresponding TV-dependent quantities (recall Eq. (12.81) ) 

F2 in/ ^°qq ^°qg I / ^2,q C2 



c = cT c7 (2 - 21) 

at /u 2 = Q 2 . Any scheme transformation can now be performed by inserting a suitable (finite, no 
relation to Eq. (12.61) above) transformation matrix Z with 

oo 

Z ik = 5 lk +£a s z Z^ (2.22) 

1=1 

which results in transformed coefficient functions C' and parton distributions q' via 

F = C-q = CZ l Zq = C'-q' . (2.23) 
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The corresponding transformation of the splitting functions reads 

P' = P+ (|3^ + \Z,P\)Z- X (2.24) 

where P represents the right-hand-side of Eq. (12.101) . and [Z,P] the standard commutator of the 
matrices Z and P. To order (N 3 LO) the insertion of the perturbative expansions thus leads to 

P' = a s P 

+ a s 2 (p( 1 ) + [z( 1 ), J p(°)]-(3oZ( 1 )) 

+ « 8 3 (/>( 2 ) + [Z( 2 \P^] + [Z« P (1) ] - [ZP\FM\ZF> + Po({Z (1) } 2 -2Z^) - frzW) 

+ at(p^ + [Z( 3 \PW} + [Z<V (1) ] + [Z<V (2) ] - ([Z<V (0) ] + [Z«P (1) ]) Z(D 
+ [Z«P (0) ] ({Z (1) } 2 -2Z( 2 )) -po({Z«} 3 -Z^Z^ -2Z®ZM + 3Z( 3 )) 
+ p 1 ({z( 1 )} 2 -2Z( 2 ))-p 2 z( 1 )) + o(a s 5 ). (2.25) 

In order to transfer the momentum sum rule, q s + g = 1 at Af = 2 for all scales, from MS to the 
transformed parton densities q' in Eq. (12.231) . the coefficients in Eq. (12.221) need to satisfy the 
relations 

Z $ (N = 2) + Z$ (JV = 2) = , Z W (JV = 2) + zg (JV= 2) = . (2.26) 

Eqs. (12.221) - (12.251) can be employed to transform to the physical evolution equations for the 
system F = (F2,F§) by choosing 



7 (0 



z ( M 

Z, qg 








^gg / 


1 = 1 




c (/) 



(2.27) 



in Eq. (12.221 ). This leads to 



ding 2 " : ^ n s ^(0 ^(0 l F 4 

where the matrix elements of the physical evolution kernel K are given by Eq. (12.251) after inserting 
Eq. (12.271) . As far as we know this transformation has first been suggested (at NLO) in Ref. [16]. 
The above relations refer to the choice /j 2 = Q 2 of the renormalization scale - the mass factorization 
scheme and scale have now been eliminated from the problem. From these results the physical 
kernel at /j 2 ^ Q 2 (or in other renormalization schemes) can be reconstructed in the usual way. As 
we shall see below, Eq. (12.271) does not represent a 'normal' scheme transformation of the quark 
and gluon distributions, since the MS coefficient functions for F2 and F§ on the right-hand- side do 
not fulfill the momentum sum rule constraints (12.261) . 
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3 Coefficient functions for Higgs-exchange DIS 



In this section we present and discuss the x-space results, obtained from our above TV-space calcu- 
lations by a by now standard inverse Mellin transformation [24,45], for the previously unpublished 

(/) (i) 9 7 

MS-scheme coefficient functions c^ q and ci ' , I = 1, 2, 3, at the standard choice /u = Q of the 

renormalization and factorization scale. 

We express our results in terms of the harmonic polylogarithms H mir .. )/n)V (x) with mj = 0, ±1. 
Our notation for these functions follows Ref. [45] to which the reader is referred for a detailed 
discussion. For completeness we recall the basic definitions: The lowest-weight (w = 1) functions 
H m {x) are given by 

H (x) = \nx , H ±l (x) = =F ln(l ^jc) . (3.1) 
The higher- weight (w > 2) functions are recursively defined as 

-^\n w x, if ra 1 ,...,ra w = 0,...,0 

H mu ..., mw {x) = { rX W - (3.2) 



I dzf mi {z)H mi ^ mw {z) , else 

J 



with 



fo(x) = -, /±iW = t^. (3.3) 
For chains of indices zero we employ the abbreviated notation 

#0,...,0,±1,0,...,0,±1,...M = H ±(m+l),±(n+l),...( x ) (3-4) 
m n 

and usually suppress the argument x. 

The Z-th order coefficient functions involve harmonic polylogarithms up to weight 21— 1 . Hence 
only the results up to two loops, with w < 3, can be expressed in terms of standard polylogarithms. 
A complete list can be found in appendix A of Ref. [24]. A Fortran program for the harmonic 
polylogarithms including weight w = 4 has been published in Ref. [46]. Its extension to w = 5, 
required for the third-order coefficient functions, is also available [47]. 

In this notation the first-order coefficient functions are given by 

4!qW = C F (2p gq (x)(-H -H 1 )-3x- 1 +2x) , (3.5) 
c^(x) = C A (4p gg W(-Ho-H 1 )-ll/3((l-x)- 1 +x- 1 )+(67/9-4C 2 )5(l-x)) 

+ n f (l/3((l-xy 1 +x- 1 ) - 10/98(1-*)) . (3.6) 

Here and below C, n represents the Riemann zeta-function, and as above denotes the number 
of effectively massless flavours. Ca and Cp are the usual QCD colour factors specified below 
Eq. (12.101) . Finally we have employed the abbreviations 
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Pgq(x) = 2x —2 + x 

p gg (x) = (1 -x)~ l +x~ l -2+x-x 2 . (3.7) 

The corresponding two-loop (NLO) coefficient functions read 

C SW = C|(4pg q (x)(-6Hi,i -2Hi,i )0 -Hi,i,i -2Hi, 2 -2H 2 , -2H 2 ,i -2H 3 ) 
+ 8p gq (-x) (-H-1C2 -2H_i_i j0 +H_i A o -H0C2) - 16(1 +x)H_i, 
-4^(2 +8x^ + 3x) +2(2-x)(5Ho A o+2H 2l o+2H2,i + 4H 3 ) 
+ 4(5-9x _1 -6x)H 2 -4£ 2 (5 -9x _1 -2x) -8^ 2 (6 -*)H -(12 -x)H Oj0 
+ (15-10x~ 1 + 18x)H +2(16-18x~ 1 - 1 x) Hi, + (56 - 64x~ ! + 3x)Hi 
+ 1/2 (208 - 59x _1 - 30x) - 32H_ 2;0 - 4Hi,i ) 

+ C A C F (20 /? gq (x)(£ 2 Hi -Hi.o.o-Hi.i.o-Hi.i.i -H 1;2 ) 
+ 4p gq (-x)(-2H_iC 2 +3H_i A o+2H_ l52 )+16(l + 2x~ ] + 2x) (£ 2 H -H 3 ) 
+ 4(2-10x~ 1 -7x)H 2) i -4(10+4x~ ! + 7x) H ,o,o - 4^ 3 (16 - 18x _1 -llx) 
+ 4/3 (24 + 22x~ l +9x + 4x 2 )H_i, +4/3 (46 -53X" 1 +x + 4x 2 )Hi j0 
- 1/9(60 + 1362JT 1 + 513x + 352x 2 )H +2/3 (70 -84X" 1 + 13x + 8x 2 )H M 
-4/3^2 (88 -53X" 1 + 7x + 12x 2 ) + 4/3 (88 - 31x~ ! + 16x + 12x 2 )H 2 
+ l/9(242 -288x~ 1 - 127x - 176x 2 ) Hi +2/3 (268 - 106X" 1 + 13 x + 32x 2 )H ,o 
+ 1/54 (6172 - 13457X" 1 + 76x + 1216x 2 ) -8(4x" 1 + 3x) H 2 , - 16H_ 2)0 ) 

+ C F n f (2/9 p m (x) (24H ,o + 29Hi + 12Hi, + 6H U + 12H 2 ) + 8/3 & (2 - 2x _1 -x) 
- 2/3 (26 - 32x _1 - llx) H - 1/27 (332 - 737 x" 1 -28*)) (3.8) 



and 



cg(x) = C\ (p gg (x) (-2570/27 +68^3 + 176/3^2 +24H_ 2 , - 778/9 H +56H C 2 
-44H ,o -28H ,o,o -778/9 Hi +40Hi £ 2 - 176/3 Hi, -40Hi, ,o - 176/3 H u 
-56Hu,o -48Hi iU -56Hi )2 - 176/3 H 2 -48H 2;0 -56H 2 ,i -56H 3 ) 
+ 4p gg (-x) (7£ 3 +6H_ 2;0 -8H_i^ 2 -8H_i -1,0 +10H_i A o + 4H_i, 2 +2H C 2 
-3H ,o,o -2H 3 ) + 16(1 +x)(2^ 3 + 4H C 2 -5H 0A0 -2H 2 , -2H 2) i -4H 3 ) 
+ 8/3 (6 + 1 lx" 1 +6x+ 1 lx 2 )H_i,o + 1/27 (10 - 5659x _1 - 1916x + 3177x 2 ) 
+ 4/3 (14 -44x _1 -25x + 33x 2 ) (Hi )0 + H M ) -4/3^ 2 (45 -44x^ - 24x + 77x 2 ) 
+ 4/3 (45 -22x _1 - 12x + 77x 2 )H 2 +2/3 (171 -44x^ -33x + 220x 2 )H ,o 
-2/9(182 + 121x" 1 -58x)Hi - 1/9(1107 + 778x~ ! +699x + 536x 2 )H 
+ (30425/162 -242/3 £3 - 778/9 £ 2 + 101/5 £f) 8(1-*)) 
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+ C A n f (p gg (x) (224/9 - 32/3 ^ 2 + 56/3 H + 8 H , + 56/3 Hj + 32/3 (H 1;0 + H M ) 
+ 32/3H 2 )+4(l +x)H 0)0 +8/3(2 -x+x 2 )^ +Hi,i) + 8/3 (3 +x 2 )(H 2 - £ 2 ) 
+ 2/9(40+48x" 1 +7x-4x 2 )Hi +2/9(47 + 4Sx^ + l4x-4x 2 )U 

+ 2/27 (389 + 320;T 1 - 46x + 59x 2 ) - (4112/81 + 28/3 ^ 3 - 56/3 ^ 2 ) 8(1 -x)) 

+ C F n f (2p gg (x) +4(1 +x)(-2& -4H £ 2 +5H ,o )0 + 2H 2j0 +2H 2 ,i +4H 3 ) 
+ 4/3 (3 +4x- 1 -3x-4x 2 )(H h0 +H U ) +4/3 (9 +4x~ l +6x-4x 2 )(H 2 - £ 2 ) 
+ 2/3 (21 +8jc _1 + 15jc-8x 2 )H ,o +2/9(93 +32x _1 -57x-68x 2 )Hi 
+ 2/9(129 +32JC" 1 -3x-68x 2 )H + 1/27 (1671 +307JT 1 - 1212* -442x 2 ) 
-(63/2-24C 3 )5(l-x)) 

+ n} (s/21p gg (x)(-5 -3H -3Hi) + 8/27(2 -x + jc 2 )(-5 -3H -3Ui) 
+ (100/81 -8/9^ 2 ) 5(1 -xj) . (3.9) 



As mentioned in the introduction, the above results have also been derived, in a quite different and 
completely independent manner, in Ref. [30]. As their counterparts for photon- and W- exchange 
DIS presented in Refs. [9, 10], the full third-order expressions, not obtained by any other group 
so far, for the present Higgs-exchange coefficient functions are exceedingly long. Therefore these 
expressions are deferred to Appendix A. 

Before we illustrate the numerical size and perturbative stability of these coefficient functions, 
it is instructive to discuss their behaviour close to the endpoints x = 1 and x = 0. For this we will 
use to abbreviations 

Xj = 1— x, Lq = In*, L\ = \nx 1 , <D k = [x^ L*]+ . (3.10) 

As usual, the +-distributions are defined by 

■l /•! 



L 



f dxa(x)+f(x) = [ dxa(x){f(x)-f(l)} (3.11) 

JO 



for regular functions f(x). It is understood that all 1 /(l — x) poles in Eqs. (13.61) . (13.91 ) and (IA.2I) 
have to be read as -^-distributions. 

(n) 

The leading large-* contributions to the gluon coefficient functions ^(x) are given by soft- 
gluon emission contributions £>, with k = 0, . . . , 2n— 1. At the first order these and all other 



endpoint contributions can be read off directly from Eqs. (13.51) and (13.61) with the help of Eq. (1371 
The +-distribution coefficients of the two-loop gluon coefficient function read 



(2) 



(2) 



8Q 2 , (3.12) 
88 2 16 

-—Cl+ — C A n f , (3.13) 
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.(2) 



,(2) 



c 



778 



25, 



A V~9 
2 /2570 



56 

T 



Ci(^-32C 3 -f &) + C^ / (^-ffc) + 2C f „ / - 



224 32 



40 

27' 



and the corresponding contributions at the third-order are given by 



,(3) 



.(3) 



.(3) 



,(3) 



25, 



.(3) 



25, 



,(3) 



'Dr 



8C 



A > 



550 



■Cf +■ 



100 



■C 4 n 



Cj (340 -128£ 2 ) - 



256 



16 



C 



/9623 1892 



3 
292 



— C A rif + — Ca«/ , 
C 2 -256C 3 ) + C A 2 n / ( 



(3.14) 
(3.15) 

(3.16) 
(3.17) 
(3.18) 



3106 344 



^2 



+ 16C> C^ny 

192268 16268 



-C A n f +-n } 



(3.19) 



Ci 



81 



C2 



2816 



^3 + 



1316 2 



+ C A n / 2 ( 



/ 67730 4192 



128 



1 9 
2 / 6652 272 



/598 



^2- — CsJ -C^Ca^— -128£ 



Ci 



81 

1616486 



9 

169910 



20^ 2 
+ — C F ny 



80 , 
27 / 



(3.20) 



40454 „ /yjl „ 2 128U„ „ \ 
£2 - k + Cl + — C2 ^3 - 80 C 5 J 



+Q 2 «/( 

+ C F C A n f f 
-C A n}( 



729 

1234307 



81 

55388 



27 
1816 



1458 
7810 
27 



81 

32C2 



^3 + 



7931 
1538 



1280 



138493 584 



1624 
152 



32 

T 



Cpflj: 

350 80 



/200 16 



1458 



^^2 + ^^+n^— --C 2 



(3.21) 



Together with the 8( 1 — x) contributions arising from soft-gluon emissions and virtual correc- 
tions - which can be read off directly from Eqs. (13.61) , (13.91) and (IA.2I ) - these results have be 
employed in Ref. [25] to obtain the soft-gluon exponentiation of Ca „ to the next-to-next-to-next-to- 
leading logarithmic (N 3 LL) accuracy. Consequently the highest seven +-distributions are known 
at the fourth and all higher orders in a s , which the exception of the (almost certainly numerically 
irrelevant) contribution of the four-loop cusp anomalous dimension to the coefficient of a"£>2«-7 - 

As discussed in Refs. [26, 48], the computation of the forward amplitudes and 7^g in 
photon- and Higgs-exchange DIS to order a s 'e with / < 3, in the above soft+virtual limit, 
facilitates the determination of all 1/e pole terms of the three-loop y*qq and Hgg form factors 
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(the latter, of course, in the heavy-top limit). The corresponding results, including the additional 
8°ny contributions in the J*qq case, have been verified recently in direct calculations of these form 
factors [43], thus providing a check of Eqs. (13.161) - (13.211) and the 8(1— x) in Eq. (IA.2I) and the 
corresponding results in Refs. [6, 7, 9] including, in particular, the full results for the three-loop 
quark and gluon cusp anomalous dimensions. Those results, in turn, involve a considerable part of 
the three-loop forward-scattering diagrams and integrals entering the complete (all-x) calculations 
of the third-order splitting functions and coefficient functions. 



Returning to the large-x behaviour of the gluon coefficient functions cjj^Jx) we note that, as 
Eqs. (13.121) - (13.141) and (13.161) - (13.181) above, the coefficients of the highest three subleading 



ln*(l- 



(2) (2) 

-x) terms, k = 3, 2, 1 for (x) and k = 5, 4, 3 for (x), do not include contributions with 
the colour factor Cf, i.e., gluon emission from quarks. Hence they are guaranteed to originate from 
'non-singlet like' diagrams with a (modulo self-energy insertions) unbroken gluon line connecting 
the incoming gluon to the scalar §. Consequently non-singlet considerations and structures hold 
for these contributions, cf. Refs. [15, 19,29], and these subleading logarithms can be predicted to 
all orders n along the lines of Ref . [20] . 



We now turn to the quark coefficient functions cjj? Jx) where the integrable logarithms h\ = 
In (1 — x) with k = 1 , . . . , 2n — 1 form the leading large-* terms . Again the first-order coefficient 



is obvious from Eq. (13.51) . and the second-order contributions are given by 



.(2) 

<t>,q 



10 



Cf Ca + - C 



F ! 



(3.22) 



,(2) 



.(2) 



3 



■C F C A - 14C F +^C F n f , 



CfCa ( 



349 



20^ +5Q 



The corresponding three-loop results read 



58 r 

Y cpn f 



(3.23) 



(3.24) 



.(3) 



L 5 , 



10 2 2 3 

— CfC a + - C F , 



(3.25) 



J3) 

'<t>,q 



.(3) 



/ 4 



T2 



61 2 
~ Ta CfC a - 

c F cl 



442 

~7n 



C f Ca 



157 
~L8" 



47 



13 



Cf + 2jC F C A n f + — 



368 



1148 
~8T 



81 

■ Cf Ca n 



C 2 )+c F 2 c A ( 



5161 296 



81 



88 



/74 

( T + T ^2 



1018 



[f 



2 ( 8213 
CfC a y 



+ Cj( 



81 

295 104 
+ — 



V 3 



81 
296 
~9~ 

C2- 



CpHjr + 



16 
27 



Cfti 



'f ' 



^2- 

184 



84^) -C}Ca(^ 



/3236 
+ C F C A n f ^ — 



608 



248 



316 n 



+ Cin 



/6616 A y \ 104 2 ^ 
+ -^2) - —njC F 



27 



(3.26) 



(3.27) 



(3.28) 



13 



.(3) 



CpC A 
638 



2 /227065 17384 
V 243 2J~ 



590 



368 



2954 



27 



, 2674 „ 

^3+^C 



r 3 



5 

/2927 
V~12~ 



,2- ( 524873 



+ 170^2 



972 

328 „ 1804 ^ 



/ 75052 3260 



52 



/ 70747 812 



15 

272 



486 



27 



c) 



/ 544 40 



-C2 

27 9 S 



(3.29) 



The two-loop coefficients (13.221) - (13.241) contribute to the large-x behaviour of the NNLO physical 
kernel (x), recall Eqs. (12.271) and (12.281) above, which will be discussed in the next section. 

The third-order results (13.251) - (13.271) will enter our large- x predictions for the four-loop splitting 

(3) 

function Pgq (x) in Section 5. 

At small x both the coefficient functions ci Ax), p = q, g, show a double-logarithmic enhance- 
ment, i.e., terms x~ l \n k x contribute with k = 0, . .., 2n—\. The second-order coefficients are 



xc 



xc 



xc 



xc, 



(2) 
(2) 
(2) 

M 

:2) 



-CfCa , 

106^ ^ 16 _ 
Cf C a + — Cf ri f , 

/454 



— CfCa{^ 



— -32^) -Ci(lO-l6^ 2 )+^-C F n f , 
13457 ^ -72fe) -C*(f -36^2 + 32 C3 ) 



(3.30) 
(3.31) 
(3.32) 



54 3 
/737 16 
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and 
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110 
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(3.35) 
(3.36) 
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Their three-loop counterpart are found to be 
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This behaviour is different from that of the (spacelike) splitting functions and gauge-boson 
exchange coefficient functions which receive an only single- logarithmic ('BFKL') enhancement 
of the leading l/x contributions to all orders in oc s [49-51]. In fact, this seemingly surprising 
feature of the tyG^G^ probe of the hadronic system was already briefly stated in Ref. [50], and 
the leading double-logarithmic terms were explicitly calculated in Ref. [52] for the related case of 
Higgs production in proton-proton collisions in the heavy-top limit. 

However, the functions cl ~{x) are probed at very small values of x only at very high centre-of- 
mass (CM) energies y/s. Given the finite physical top mass, the large top-mass approximations for 
Higgs exchange/production breaks down in this region. Indeed, as explicitly shown recently for 
Higgs production in Ref. [53], the full finite-m t op coefficient functions receive only single logarith- 

(n) 

mic contributions at small x. Consequently the effective-theory coefficient functions '(x) cease 
to provide useful approximations for Higgs-exchange DIS in this limit. Where exactly in x this oc- 
curs cannot be established without the very non-trivial calculation of at least all NNLO l/x terms 
in the full theory. This point has not been reached yet even in the phenomenologically important 
case of Higgs production in proton-(anti-)proton collisions, for the present status see Refs. [54]. 
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We are now ready to discuss the numerical size of the coefficient functions (13.51) - ( 13.91 ), (IA.1I) 
and (IA.2I) . In Figs. 1 and 2 their cumulative effect is shown for «y = 4 quark flavours and an order- 
independent value 0C s (/j 2 = /j 2 = <2q) = 0.2 for the strong coupling constant corresponding to a 
scale of « 20 ... 50 GeV 2 . The same reference point was used in Refs. [8,9], hence the present 
results can be compared directly to the photon-exchange case discussed in those articles. Note 
that Fig. 2 does not provide the full information about the effect of the gluon coefficient function 
C&o(x) due to presence of +-distributions, recall Eq. (13.1 II) . and 8(1— x) contributions. 

As shown in the right parts of both figures, the perturbative expansions of C^ p (x) are rather 
stable in the mid-* region, but, unsurprisingly, large beyond-NLO corrections are found in the 
soft-gluon large-* and the high-energy small-* regions. While in the former the NNLO and N 3 LO 
are similar, the perturbative expansion appears to break down below x ~ 10 ~ 3 and x ~ 10 ~ 4 for 
Cb q(x) and C^Jx), respectively. This breakdown thus occurs in a region where this coefficient 
functions are definitely not expected to represent Higgs exchange anymore, as a CM energy of 
about m top is required to access x — 10 ~ 3 at the chosen scale Q 2 ~ 30 GeV 2 . 

Obviously the (perturbative) coefficient functions enter physical quantities such as the structure 
function F§ only through convolutions such as 



with the (partly non-perturbative) quark and gluon distributions which, due to the general shape of 
q s (x) and g(x), shift the onset of the small-* instability to considerably lower values of * than in 
Figs.l and 2. The resulting contributions of Eqs. (13.51) - ( 13.91 ), (IA.1I) and (IA.2I) to F§ are illustrated 
in Fig. 3 for the schematic, but sufficiently realistic order-independent input 



already employed in Refs. [7-9]. In reality both oc s and the parton distributions do depend on the 
perturbative order in a manner that reduces the relative higher-order corrections in particular at 
large *. Keeping this in mind, the medium-<2 2 perturbative stability in Fig. 3 is satisfactory over a 
wide range in * despite, as expected for a gluon-dominated quantity, the presence of considerably 
larger a s corrections than found for the standard-DIS structure function F2 . 

Finally we need to address the Mellin moments of the new coefficient functions c^(x) for 
n = 1, 2, 3 (p = q, g). The second moments for the quark case are given by 




(3.50) 



xq s {x,Qo) = 0.6jc-°- 3 (1-jc) 3 - 5 (1+5.0* - 8 ) , 
xg(x,Ql) = 1.6jc- a3 (l-*) 4 - 5 (l-0.6* - 3 ) 



(3.51) 




1.77778 , 



(3.52) 
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32^ 3 ) +cj 
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-32ty+ — C F n 



81 
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194.729 + 19.5556/j 



7> 



(3.53) 
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Figure 1: The perturbative expansion to order oc s 3 (N 3 LO) of the quark coefficient function C^Jx) 
for DIS via the exchange of a scalar (|) with a tyG^G/jy coupling to gluons. The results are shown 
for a standard medium-scale reference point and have been multiplied by x for display purposes. 




Figure 2: As Fig.Q], but for the gluon coefficient function C§ jg (x) at x < 1, multiplied by an addi- 
tional factor ( 1 —x) compensating the large x-dependence due to the +-distribution contributions. 
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Figure 3: The perturbative expansion of the quark (left) and gluon (right) contributions to the 
structure function F§, normalized to the lowest-order result F$(x, Q 2 ) = xg(x, Q 2 ), for the order- 
independent input (13.511) at our reference point corresponding to a scale Q 2 ~ 20 ... 50 GeV 2 . 
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= -12116.73 + 1902.16^-50.5213^ (3.54) 

where the second (approximate) equalities are obtained by inserting Ca = 3, Cp = 4/3 and the 
numerical values of the ^-function. The results for the gluon coefficient function read 

18 9 



c£l(N=2) = —C A - — n f = 33.8333- 1.11111^, (3.55) 



4>^ - ^ (w-«&Ww +12 & 



„ / 1726 nAt , \ 100 2 



1242.06-185.349^+1.23457^ , (3.56) 
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(3),„ ^ w 27878711 66970,, 880,, 
4>(N = 2) = Q 3 (^2 27" ^ 3 + — ^ 5 

b A f \ 2916 Si 15 S2 3 s - 
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Note the conspicuous absence of coefficient linear in ^2 from the above second- and third-order 
results, a feature that has been observed and discussed before for the corresponding moments of 
the gauge-boson exchange structure functions in Refs. [23]. 

For rij = 4 flavours Eqs. (13.521) - (13.571) yield the rather benign expansions 

C^ q {N=2) = - 0. 1415 oc s - 0.7378 oc s 2 - 2.6791 oc s 3 +... , (3.58) 
C +)g (JV=2) = 1 + 2.3387 a s + 3.2956a 2 + 5.3704a 3 + . . . (3.59) 

which nevertheless, as expected, are rather different from their counterparts for the photon-exchange 
structure function F2 given by [23] 

C 2) q(N = 2) = 1 + 0.0354a s - 0.0785a s 2 - 0.1986a 3 + ... , (3.60) 
C 2 , g (iV = 2) ^ -0.1592a s - 0.2259a 2 - 0.0274a 3 +... . (3.61) 

Consequently, as already indicated at the end of Section 2, the scheme-transformation relation 
(f2T26h is not fulfilled with the choice (I2T271) involving Eqs. (13381) - (I33TT) . 

The Af-dependences of C§, q and C^g are shown graphically, again for a s = 0.2, in Fig. 4. The 
corresponding results for F2 and F F in standard DIS have been presented in Figs. 11 and 12 of 
Ref. [9]. The present pattern is rather similar to that for F2 (with the quark and gluon coefficient 
functions interchanged), albeit with considerably larger a s corrections. In fact, the ,/V-dependent 
values of the strong coupling, 

for which the n-th order corrections are half as large as those of the previous order - recall that 
the coefficient functions are expressed in terms of the small expansion parameter a s = a s / (4ft) 
of Eq. (12.101) - shown in Fig. 5 are rather more similar to those for Ft than for F2 which can be 
found, respectively, in Figs. 13 and 14 of Ref. [9]. Note, however, that also in Fig. 5 the values for 
(//) increase, rather than decrease, with the perturbative order n. This behaviour indicates that 
also here we have not yet reached the expected asymptotic regime of the expansion in powers of 
a s with a factorial-type growth c'"' ~ n\ of the expansion coefficients. 
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5 10 15 5 10 15 

N N 

Figure 4: The perturbative expansion to order a 3 (N 3 LO) of the quark and gluon coefficient func- 
tions for the structure function in Mellin-./V space, shown for four effectively massless flavours 
and the order-independent value of oc s specified in the right panel. 




5 10 15 5 10 15 

N N 

Figure 5: The Af-dependent values (13.621) of a s at which the effect of the n-th order (N"LO) quark 
and gluon coefficient functions for F§ is half as large as that of the previous order. A NLO curve 
can only be shown for the gluon case (right) since only here the LO contribution does not vanish. 
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4 Physical evolution kernel for F 2 and 



We are now in a position to write down and discuss the physical kernels for the coupled flavour- 
singlet evolution of the structure functions F 2 and F§ as defined at the end of Section 2. To NNLO 
(1 = 2) these kernels, 

Kab = Y. a i +lR al with a,b = 2,$, (4.1) 

1=0 

are specified in terms of the three-loop splitting functions [6, 7] together with the second-order 
coefficient functions for F2, first calculated in Refs. [55], and F§ provided in the previous section. 
Also these results will by given in terms of harmonic polylogarithms [45] as summarized in Eqs. 
( 13.11) - (13.41 ) above. From / = 1 we will present K22 in the form 

4 , = < + 4 , , s . (4.2) 

i.e., we separately provide the scalar evolution kernel for the non-singlet structure function F% ns 
(addressed, but not written down in this form, before in Refs. [20, 56]), which is of direct phe- 
nomenological interest for, e.g., the determination of 0C S from scaling violations. As above, our 
results are given in the MS renormalization scheme for the scale /j = Q. 

At leading order (LO) the elements of the physical kernel are identical to the splitting functions, 
recall the first line of Eq. ( 12.241 ), i.e., 

= C F (2p qq (x) + 35(l-x)), (4.3) 

= 2n fPqg (x), (4.4) 

4? = 2C FPgq (x), (4.5) 

K$ = C A (4p gg (x) + y8(l-x))-^n / 8(l-x), (4.6) 
where we have used Eq. (13.71) and the additional abbreviations 

p qq (x) = 2(l-x)~ 1 -1-x 

p qg (x) = 1 - 2x + 2x 2 . (4.7) 

The next-to-leading order (NLO) contributions to these kernels are given by 

K 2%(*) = C f (Pqqto [-6 H + 8 Hj.o + 8 H 2 ] +S Pqq (-x) [-C 2 -2H_i |0 +H 0j0 ] 
+ 2 (-2 +2x +H - 3xH - 2H ,o - 2xH , ) + [3/2 - 12^ 2 + 24^ 3 ] 5(1— jc)) 

+ C A C F (p<n(x) [367/18 +2/3 (-6^2 +22H +6H ,o + 11 Hi)] 
+ 4p qq (-x)[C2+2H_i ) o-H 0) o] +1/6(13 -167x + [215 + 176^2-72^3] 8(l-x))) 

+ C F n / ('p qq (x)[-29/9-8/3Ho-4/3Hi] +1/3(1 +13x-[19 + 16C 2 ]5(l-x))) , (4.8) 
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K 22, P s( x ) = On / (2/3(32jr 1 -69 +45x-8x 2 ) + 4/3(4;T 1 -9x+Sx 2 )H Q 

-8(l+*)Ho,o) > (4-9) 

K$(x) = C A n / (8/7 q g(x)H 2 -8/7 qg (-x)H_i ; o +2/3(1 -50x + 50.* 2 )Hi 
-8(1 +2jc)H ,o - 8^2 (1 +2jc 2 )+2/3(8jT 1 -11 -98x + 66x 2 )H 
+ 2/9(100j»T 1 -322 + 155x+18x 2 )) + C F n f (4 Pqg (x) [2 H ,o -3 Hi + 2Hi, ] 
-4(1 -2x)H ,o+2(l +4jc-12x 2 )H +2(2-7x + 4x 2 ) + 16Hi) 
+ n / 2 (4/3p qg (x)[H +H 1 ] +4/9(-2x- 1 + ll -25x + 21x 2 )) , (4.10) 

^ W = C| ^8 p gq (x) H 2 + 8 ^2 (-2JT 1 + 2 - x) - 4 (2 -x) H ,o + 2 (2 + 3x) H 
+ (l(br 1 -18 -x) + l2 X - 1 nA + C A C F ^ gq (x)[4/3H 1 +8Hi, ] 

- 8 pgq(-x) H_i,o - 6 (2 - x) Hi + 8 (2 +x) H ,o + 1/3 (87JT 1 - 22 + 83x + 64x 2 ) 
+ 2/3 (22*- 1 - 138 + 15x- 16x 2 )H +8^2 (2x _1 +x)\ 

+ C F n / (p gq (x)[-4H -4/3Hi] -2/3 (29x~ l -26 + 13+)) , (4.11) 

K$(x) = Q 2 (p gg (x)[389/9 -8^ 2 +44/3Ho+8Ho ! o+44/3H 1 +16H 1 ,o + 16H 2 ] 

- 8 p gg ( -x) [ £ 2 + 2 H i , - H ,o ] + 32 (1 + x) H ,o - 4/3 (25 - 1 1 x + 44 x 2 ) H 

- 1/9 (268X" 1 - 728 + 607 x - 389 x 2 ) - (449/27 - 44/3 £ 2 - 12 £ 3 ) 8(1-*)) 

+ C A n f U/3p gg (x) [-7-2H -2Hi] -8/3(1 +x)Ho -4/9 (Or 1 + 13 -2x-2x 2 ) 
+ 4/27 (43- 18 £ 2 ) 8(1-*)) + C F n f (-8(1 +jc)H ,o -4/3 (4a:" 1 +6 -9jc)H 

- 2/9 (44x _ 1 - 27 - 45jc + 28jc 2 ) - 28(1 -x)) 

+ «/ (4/9pgg(jc) +4/9(2 -x+x 2 )- 20/27 8(1 -x)) . (4.12) 

The five corresponding third-order (NNLO) corrections fill about 10 pages (in total) and therefore 
are deferred to Appendix B. 

The leading large-x contributions to the evolution kernels (14.81) - (14.121) and (IB.ll) - (IB .51) are 
given by terms of the form CD k (defined in Eq. (13.101) ) for the diagonal entries and , and 

h\ = \n k (l— x) for their off-diagonal counterparts and . It turns out that in all cases the 
higher-order enhancement (as far as known, obviously) is only single-logarithmic, k = 0, . . . , n. 
This is in contrast to the behaviour of the off-diagonal MS splitting functions, cf. Eqs. (4.17) - 
(4. 19) of Ref. [7], and of all coefficient functions for F 2 and F§ as discussed in Ref. [9] and Section 3 
and Appendix C of the present article. Hence the key non-singlet observation of Refs. [19,20] 
generalizes to the present flavour-singlet case. We will utilize this fact in the next section. 



23 



The NLO and NNLO large-* coefficients for K22 (x) are given by 
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The corresponding terms for K^} (x) and K$ (x) read 
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Note that the leading large-* coefficients at all three orders are simply given by (— Po) w Ai )P , where 
Ai ;P denotes the appropriate one-loop quark or gluon cusp anomalous dimension, i.e., the coeffi- 
cient of (l-x)~ l in Eq. (H3T) or (l4T6b. 



The log-enhanced large-x contributions to (x) and (x) are 
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The corresponding terms for K^ix) are given by 
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Thus already the leading large-jc terms are of a more complicated form in the off-diagonal cases. 



At small-* these physical evolution kernels, unlike the MS splitting functions [49], include 
double-logarithmic 1 /^-contributions originating, via Eqs. (12.251) and (12.271) . in the corresponding 
behaviour for the coefficient functions c*f\ discussed in the previous section. Specifically, the 



quantities K^{x) and K^Mx) are found to include terms up to x Hn 2 " l x for n = 1, 2, 3. The 



-(«) 



corresponding leading small-x contributions to the N 2 LO and N 3 LO kernels (x) and K^j} (x) 



are of the form x l \n 2n 2 x. Since also these terms are artifacts of the large-m t0 p approximation 
discussed above, we will not write them down here. 
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Instead we directly proceed to a brief discussion of the numerical implications of Eqs. (14.81) - 
(14.121) and (IB.ll) - (IB.5I) . The resulting perturbative expansions (14.11) of the elements K a /, of this 
evolution kernel are illustrated in Figs. 6-9. Here the initial conditions F 2 (x, Qq) and FJx, Qq) 
are physical quantities and consequently, unlike the MS parton densities in Section 3, in principle 
independent of the perturbative order. Hence it is easier here to provide a realistic (and not only 
schematic, as in the previous section) account of the perturbative (in-) stability of the quantities 
under consideration. 

The only remaining order-dependent quantity is the MS coupling at our chosen reference scale. 
Here we employ the NLO and NNLO values obtained from a schematic analysis of non-singlet 
scaling violation in Refs. [9, 56], supplemented by a corresponding LO result, 

a s (£ 2 )LO = 0-255 
«s(2 2 )nlo = 0.208 

a.s(<2o)NNLO = O- 201 - (4-29) 

These values are consistent with a NLO — NNLO difference of about 0.002 to 0.003 at the mass of 
the Z-boson as found in recent fits of structure functions and related data [57]. The actual scale Qq 
in Eq. (14.291) depends, of course, on the value of a s (M|) with, at present, 20GeV 2 < Qq < 50GeV 2 . 

For the illustrations of the resulting scaling violations in Figs. 10 and 1 1 we use the sufficiently 
characteristic form (13.511) for the structure functions at the reference scale Qq, i.e., we employ 

F 2 {x,Ql) = 0.6x-°' 3 (l-x) 3 - 5 (l+5.0x - 8 ) , 

F*{x,Qo) = 1.6jc-°- 3 (1-jc) 4 - 5 (1-0.6jc ' 3 ) . (4.30) 

Note that, as implied by the identity of Eqs. (13.511 ) and (14.301 ), we have absorbed the average (e 2 ) 
of the quark charges into the definition of the singlet structure function F 2 , which is thus given at 
leading order by F 2 = xq s instead of the standard normalization F 2 = (e 2 )xq s . 

At medium to large x especially the diagonal entries K 22 and K^, respectively shown in Figs. 6 
and 7, are rather stable, with the relative NNLO corrections not exceeding 10% in the range 0.2 < 
x < 0.9 in both cases. Larger higher-order effects are found in this region for the (absolutely 
smaller) off-diagonal kernels displayed in Figs. 8 and 9. In all cases the perturbative expansion 
appears ill-behaved at small values of x, starting from about x ~ 10 1 for K§§ to x ~ 10~ 3 for 
K§ 2 . We expect that a (physically realistic) inclusion of the finite-m to p effects as discussed below 
Eq. (13.491) would fundamentally change these small-* results. 

Due to the ubiquitous Mellin convolutions (13.501) - and the diagonal distribution effects not 
visible in Figs. 6 and 7 - even the present evolution of the structure functions F 2 and F§, illustrated 
in Figs. 10 and 11, is perfectly stable down to about x = 10 and rather acceptable at x-values as 
low as x ~ 10~ 3 . At large x these physical scaling violations are, unsurprisingly, much stronger 
(by a factor between 2 and 3) for the gluonic observable F§ than for F 2 . Here the relative NNLO 
corrections exceed 5% only at x > 0.7 for the scale derivative of F 2 , and at x > 0.6 for that of F§. 
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Figure 6: The perturbative expansion to order a s 3 (NNLO) of the diagonal element A^C*) of the 
physical evolution-kernel matrix for the system (F2, Fa) of flavour-singlet DIS structure functions 
at x < I. The poles at x = and x = 1 have been suppressed by multiplication with x( 1 — x) . Also 
shown is the pure-singlet contribution K22, ps defined in Eq. (14.21) . 
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Figure 7: As Fig. [61 but for the second diagonal element K^(x). As in all figures in this section, 
the order-dependent values ( 14.291 ) are employed for the strong coupling a s . 
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Figure 8: The perturbative expansion to order oc s 3 (NNLO) of the off-diagonal element K2^(x) of 
the physical evolution-kernel matrix for the system (F2, F§) of singlet structure functions. 
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Figure 9: As Fig.[H but for the second off-diagonal element K^x). In both cases the curves have 
multiplied by a factor of x compensating the main x-dependence. 
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Figure 10: The perturbative expansion (12.281) to / = 2 (NNLO) of the normalized scale derivative 
d In F2 jd In Q 2 of the photon-exchange singlet structure function F2 for the initial conditions (14.301) 
at a mid-scale reference point with the coupling constants (14. 291 ). 
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Figure 1 1: As Fig.[l0l but for the scaling violations of the scalar-exchange structure function Fa . 
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5 Large-* behaviour and four-loop predictions 



We now return to the expansion of the physical kernels and the MS splitting functions in powers 
of L\ = ln(l— x). The leading large-* terms, (1— x)~ l L\ for K22 and K§§ and (1— x)°L^ for K2§ 
and K§2, have been given to order oc s 3 (NNLO) in Eqs. (14.131) - (14.281) . In fact, the only single- 
logarithmic enhancement - k < n at N' 7 LO - does not only hold for these leading contributions, 
but to all orders in (l—x) for all matrix entries. Thus the physical evolution kernels fulfill 



ab 



/»+*/(*) 



for k > 1 , (5.1) 



to (at least) n = 2, where f(x) denotes functions with a non-logarithmic expansion at x = 1 . 

Already the leading L"f(x) terms are not particularly simple for the off-diagonal evolution 
kernels, hence we refrain from writing them down here for brevity. The observation made below 
Eq. (14.221) for the diagonal entries, on the other hand, also generalizes to all orders in (1—*), i.e., 
the leading terms are those of the non-singlet kernels, and therefore given by [20] 

K${x) = ln(l-x)p qq (x)[-2C F p -8C|H ] + (m°(l-x)) , (5.2) 

K$(x) = \n 2 (\-x)p qq (x)[2C F $ 2 + 12C2(3 Ho + 32C^Ho,o] + o(ln(l-x)) (5.3) 

and 

K$(x) = ln(l-x)2p gg (x)[-2C A ^ -SC 2 H } + o(ln°(l-x)) , (5.4) 

K$( x ) = ln 2 (l-x)2p gg W[2C A p 2 + 12C A 2 (3oHo + 32C A 3 Ho,o] + o(\n(l-x)) . (5.5) 

Note that the additional factor of 2 in the latter two equations arises from a mismatch, in terms 
of the normalization of the leading (1— x)~ terms, of the definitions of p qq (x) and p gg (x) in 
Eqs. (13.71) and (14.71) . These results imply that even the L"f(x) contributions to the pure-singlet 
quark and gluon kernels (recall the discussion in the paragraph above Eq. (13.221) ) vanish to at least 
n = 2, 

= . (5.6) 



22,ps 



= K 

L»f(x) * 



C F L"f(x) 



The above behaviour of the physical evolution kernels is in striking contrast to that of the 
MS-scheme splitting functions, which do include double-logarithmic contributions, ln^( 1 — x) with 
k > n > 1 at N"LO, at all orders in the expansion at x = 1 . At NLO these contributions, and the 
single-logarithmic terms of the diagonal splitting functions, are given by 

P${x) = \n 2 (l-x)C AF n f {-4p qg (x)} + o(ln(l-x)) , (5.7) 

P g { q \x) = \n 2 (\-x)C AF C F {+4p gq (x)} + 0(ln(l-x)) (5.8) 
with Ca F = Ca~C f , and 

P&\x) = ln(l-jc)C|{-8p qq (x)}H + 0(ln(l-x)°) , (5.9) 

P g {l g \x) = ln(l-Jc)C A 2 {-16p qq W}H + 0(ln(l-jc)°) . (5.10) 
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The corresponding NNLO results read 

ln 4 (l-x)C| F n / 4/3p qg (x) 



and 



+ ln 3 (l-x) C A >n / {-(34+16/x-92x + 64x 2 )/9-8/3(l+4x)H } 

+ C AF C F n f {-40/9 - 16/(9x) - 28/9 x + 80/9 x 2 + (4 - 24x + 32/3 x 2 ) H } 



+ C AF n / 2 4/9 p qg (x) J + 0(ln 2 (l-x)) 

= ln 4 (l-x)C A 2 F CV4/3p gq (x) 
+ ln 3 (l- 



(5.11) 



cl F C F {(344/x- 316+ 170*- 16x 2 )/9-8/3 (4/x-8+x)H } 
+ C AF C% {128/(9*) -76/9- 16/9 (x + x 2 ) - (32/(3x) -24 + 4x)H } 
+ C AF C F nx (-20/9 p gq (x))l + o(ln 2 (l-x)) (5.12) 



p q (2) ( x ) = hi 3 (l-x) [c AF C F /^ {-32/(9x)- 8/3 + 8/3x + 32/9x 2 - 16/3(1 +x)H } 

+ ln 2 ( 1 -x) C| {2/3 ( 16/x + 7 - 7x - 16x 2 ) + 8 (3 + 4x + 4/3 x 2 ) H } 
+ C F C A n f {-S/(9x) - 6 + 6x + 8/9x 2 + 4/3 (4/x - 1 1 - 17x - 4x 2 )H } 
+ C F nj {16/(9x) +4/3 -4/3x- 16/9x 2 + 8/3 (1 +x)H } 



+ C AF C F n f {8 (l+x)(H 0) o-2Hi )0 )}J + Cj 32 /7 q q(x)H 0) o 
+ 0(ln(l-x)) , 



(5.13) 



i> g ( g 2) (x) = ln 3 (l-x)[c AF C F /^{32/(9x) +8/3 - 8/3x- 32/9x 2 + 16/3 (l+x)H } 

+ ln 2 (l-x)[c^/i / {-2/3(16/x + 7-7x-16x 2 )-8(3 + 4x + 4/3x 2 )H } 

+ C F C A n f {8/(9x) +6 -6x-8/9x 2 - 4/3 (4/x - 11 - 17x-4x 2 )H } 

+ C F ^{-16/(9x)-4/3 + 4/3x+16/9x 2 -8/3(l+x)H } 

+ C AF C F n / {-8(l+x)(H 0) o-2Hi j0 )}] + C 3 64p gg (x)H , 
+ o(ln(l-x)) . (5.14) 

The respective final terms of the last two equations represent the non-singlet contributions directly 
corresponding to the final terms in Eqs. (15.31) and (15.51) . 

Eqs. (15.71) - (15.141) show a couple of interesting features: The double-logarithmic contribu- 
tions, unlike the single-logarithmic terms, all vanish for the case C F = C A which forms a subset 
of the colour factor conditions leading to an j\£ = 1 supersymmetric theory. In fact, the coeffi- 
cients of the leading, next-to-leading etc double logarithms successively include terms of the form 
(C A — C F )' l ~ l , (C A — C F )' l ~ 2 etc. Moreover, to the accuracy of Eqs. (15.131) and (15.141) . the pure- 

(2) (2) 

singlet parts of Pqq and Pg g J turn out to be completely identical up to an overall sign. 



31 



Both the single logarithmic enhancement of the physical kernels also in the present flavour- 
singlet case, and the above features are non-trivial (and not obviously related) properties. In the 
non-singlet case the former property is known to hold also at N 3 LO [20]. Recall also that it is 
implied [56] for the leading (1 — *)~ by the soft-gluon exponentiation [21], and established at all 
powers of (1 — *) in the large-n^- limit by Refs. [58, 59]. 

Hence it appears very natural to conjecture that also Eqs. (15.11) and Eqs. (15.61) are valid at 
N 3 LO as well. As the three-loop coefficient functions are known from Refs. [9] and Section 3, this 
implies that the known double logarithms originating from the matrix (12.271) in Eq. (12.251) at order 

have to be compensated by the (hitherto unknown) corresponding contributions to the four- 
loop splitting functions P^ 3 K These terms can therefore be predicted from the above conjecture, 
and checked against the colour-factor pattern expected from Eqs. (15.71) - (15.141) . 

For the off-diagonal splitting functions these predictions are found to be 



= ln 5 (l-*) CXf n/{32/(27x) +46/27 + 52/27 *-4* 2 + 16/9 (1+4*) H } 
+ Cl F C F n f {32/(27jc) + 58/27 + 100/27 x - 68/9 x 2 + 8/9(1 + 10*) H } 
-C 2 F nj {4/27 p qg (x)} 

+ ln 4 (l-x) \cl F n f {152/(27*)- 106/9+ 1178/27*- 719/27 x 2 

- 2/9 (8/(3*) — 1 + 16jc— 11 14/3 x 1 ) H - 80/9 ( 1 + 2x + 2x 2 ) H_i, 

- (104/9 + 160/9*) H ,o + 8/9 (1 1 + 38* + 2* 2 ) H 1)0 - 80/9 ^ 2 P qg (*)} 
+ Cl F C F n f {-176/(81*) +685/162- 1426/81*+ 1166/27 * 2 

- 1/27 (32/*-71+742*-2156* 2 )H -80/9(l+2* + 2* 2 )H_i !0 
+ 8/9 (31 - 2* + 52* 2 ) Hi, + 16/9 (7*+ 10* 2 ) H ,o - 8 ^ 2 P qg (*)} 

+ cl F nj {8/(27*) -98/81 + 184/81 *- 226/81 * 2 + 4/9 (1 +4*)H } 
+ C AF Cjn f {-632/(81*) -877/81 + 1349/81 * + 49/27* 2 

-l/27(16/* + 211+28* + 552* 2 )H -(10 - 220/3 * + 80/3 * 2 ) H 0)0 } 
+ C AF C F nj {-40/(81*) +719/81 - 1060/81 * + 398/81 * 2 

+ 2/27 (29 + 50* - 52* 2 ) H + 8/9 (2 - 29*) H , } 



Cafti} {4/81 p qg (x)} + 0(ln 3 (l-*)) 



(5.15) 



and 



Pgq j (*) = ln 5 (l-*) Cj F C F {(292/*-236 + 46*-32* 2 )/27 + 16/9(4 + *)H } 
+ Cl F C} {(148/*-44-86*-32* 2 )/27 + 8/9(10+*)Ho> 
-Cl F C F n f {4/27 p gq (*)} 
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+ In 4 ( 1 -*) Cl F O {3029/ (8 1 *) - 428/8 1 + 295/8 1 * + 128/27 x 2 

- ( 1168/jc - 2300 - 68jc + 176* 2 )/27 H + (80/(3 x) - 400/9 + 16*) H 0j0 
+ 80/9 (2/* + 2 + *) H-1.0 + 8/9 (26/ x+ 14 + 23*) Hi, + 16/9 ^ 2 PgqM} 

+ C^Q? {7142/(81*) - 1418/27 - 505/54 x + 736/8 lx 2 - (1472/(27 x) 

- 1 022/9 - 89/9 x + 64/9 x 2 ) H + ( 1 60/ (3*) - 1 28 + 280/9 x) H ,o 
+ 80/9 (2/* + 2 + *) H_ 1;0 + 8/9 (52/ x - 2 + 3lx) Hi, - 8 ^ 2 P gq (*)} 

+ Cf F C F n y { -2290/ ( 8 1 x) + 656/27 - 70/9 x + 1 84/8 1 x 2 
+ (8/(3x) - 440/27 - 56/27 jc) H } 

+ C AF C^ {47/(3 jc) - 1855/81 + 185/81 x + 352/81 x 2 

- (304/x - 424 + 17* + 16jc 2 )/27 H + (80/(3*) - 220/3 + 10*) H ,o} 
+ C AF C%n f {-(1474/*- 1460 + 217*-248* 2 )/81 

+ (8/* - 532 + 62*)/27 H + 16/9 (2 - *) H ,o} 
+ C AF C F n} {32/81 p gq (*)}] + 0(ln 3 (l-*)) . (5.16) 

Here, and in Eqs. (15.191) and (15.201) below, we have employed the modified basis of Ref. [20] for 
the harmonic polylogarithms, in which all logarithms and ^-functions are made explicit in the 
expansion at * = 1 to all orders. The corresponding functions entering the above expression are 

H_i,o(*) = H_i, (*) + C2/2 , Hi,o(*) = Hi )0 (*) + (5.17) 

and, e.g., Ho,i(*) has been expressed by 

H j(*) = -H 1)0 (*) - ln(l-*)H (*) - C2 • (5.18) 

The corresponding predictions for the four-loop pure-singlet splitting functions (which, as ex- 
pected, are suppressed by another power of ( 1 — *) for * — > 1) read 

fqq,psM = ln 4 (l-*) \cl F C F n f {-12%/(%\x) -602/27 + 122/27 * + 2168/81 * 2 
+ 4/27 (40/* - 1 2 1 - 24 1 * - 40* 2 ) H + 80/9 ( 1 + *) H ,o} 

+ C AF C} n f { -248/ (8 1 *) - 62/27 + 62/27 * + 248/8 1 * 2 
+ 4/27 (40/* - 1 - 6 1 * - 40* 2 ) H + 1 60/9 ( 1 + *) H , } 

+ C AF C F nj {176/(81*) +44/27 (1 - jc) - 176/8 1* 2 + 88/27 (1 +*)H } 

+ In 3 ( 1 -*) \cj n f { 1 072/ (27 *) + 6 1 0/9 - 2 1 82/27 * - 80/3 * 2 

- (320/(27*) + 80/9 (1 -*) - 320/27 * 2 ) £2 + (2252/27 - 884/27* 

- 2560/27 x 2 - 160/9 C2 ( 1 + x) ) H - 8/3 (5 + 25/3 * + 1 52/9 * 2 ) H , 
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+ (1216/(27*) - 112/3 - 1264/9*- 1216/27 * 2 )Hi, 

- 320/9 ( 1 +x) Hi ,o,o + 32 ( 1 + *) H , i ,o + 560/9 ( 1 + *) H ,o,o} 
+ C}C A n f {-24232/(243*) + 8954/81 - 16730/8 lx + 47560/243 * 2 

+ (64/(27 x) + 16/9(1 -x) - 64/27 x 2 ) & ~ (4864/(81 jc) +2584/27 

- 4952/27 x - 6944/8 1 x 2 - 32/9 (1 +x) £ 2 ) H - 16/27 (142 
+ 30 1 x) H ,o + (640/(27 x) - 160/9 ( 1 +x) + 640/27 * 2 ) H_ 1 , 

- 16/27 (88/* - 205 - 457* - 88* 2 ) Hi, - 320/9 ( 1 -x) H ,-i,o 

- 128/3 (1 +x)Ho,i,o- (176/3 - 208/9*) H ,o,o} 
+ C F cln f {1580/(27 x) - 13028/81 +23804/81 x- 1724/9* 2 

+ (256/ (27 x) +64/9 ( 1 — jc) — 256/21 x 2 ) ^ 2 + (2720/(81*) +632/27 

- 3656/27* + 3872/81 x 2 + 128/9 (1 +jc) C2) H - (832/(27*) - 1376/27 

- 4064/27 * - 1 28/9 x 2 ) H ,o - 1 60/9 (4/ (3*) - 1 - * + 4/3 x 2 ) H_ 1 , 
+ (64/(9jc) -2272/27 - 3520/27* -64/9* 2 )Hi, + 32/3 (1 +x)H 0) i,o 
+ 320/9 (1 -*) H -i,o + 320/9 (1 +*) Hi A o - (32/9 + 256/3*) H ,o,o} 

+ C}n 2 f {-(1504/*+ 1820- 860*- 2464* 2 )/81 - (368/9 + 160/3* 

+ 704/8 1 x 2 ) H + 1 6/27 ( 1 + *) H ,o - 352/27 ( 1 + *) H x , } 
+ C F C A nj {3280/(243*) + 1472/81 -512/81*- 61 60/243 * 2 

- (704/ (8 1 *) - 832/27 - 448/9 * - 704/8 1 * 2 ) H 

- 32/27 (13 + 10*) H , + 352/27 ( 1 +*) Hi, } 

+ C F nj {-128/(81*) - 32/27 (1 -*) + 128/81 * 2 - 64/27(1+*) H } 
+ 0(ln 2 (l-*)) (5.19) 



and, in perfect agreement with the NNLO observation below Eq. (15.141) . 

-P$M + 0(ln 2 (l-*)/(x)) . (5.20) 



C F 

Eq. (15.191 ) includes three more modified basis functions, 

Ho,-i,oW = H ,-i,o(*) +H (*)W2 + 3^3/2, 
H ,i,oW = H ,i,o(*) + H (*) ^2 + 2C 3 , 

Hi,o,oto = Hi )0) oW - • (5.21) 

Both the full quark-quark and gluon-gluon splitting functions differ from the above results 
by a term with one unknown coefficient, respectively given by ^ p (3)C^ln 3 (l— *) p qq (x)Ho,o,oW 
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and 4p(3) C\ In 3 ( 1 — x) 2 p g% (jc) Ho,o,o W , which directly enters the large-x behaviour of the diagonal 
physical kernels, given by [20] 



K 



(3) 
22 



p qq W [-2C F p 3 - 44/3C2p 2 Ho-64C^poH 0; o + ^ 3 C;Ho,o,o] , (5.22) 



and correspondingly 



K 



(3) 



2/> gg (x) [-2Qp 3 " 44/3C A 2 (3 2 Ho-64C A 3 poHo,o + ^ 3 ClHo,o,o] ■ (5.23) 



Obviously the N 3 LO predictions (15.151) and (15.161) . and the double-logarithmic part ~ ln 4 ( 1 — x) 
of Eq. (15.191) , exactly follow the colour-factor pattern expected from the NLO and NNLO results 
(as discussed below Eq. (15.141) ). In our view this represents, for the time being, an acceptably 
strong check of these predictions and thus of the underlying conjecture about the four-loop physical 
kernel matrix for the structure functions (F2,F§). Furthermore we expect that a formal proof of this 
property, presumably first for subleading the (1—x) In (1— x) terms corresponding to iV _1 lniV 
contributions in Mellin-/V space, can be obtained using, for instance, the recent path-integral ap- 
proach of Ref. [60] or soft-collinear effective theory (SCET) [61]. For other recent research into 
N~ l InN contributions to hard processes see Refs. [62]. 

In fact, the (1— x)°ln k (l — x) contributions to the four-loop off-diagonal splitting function form 
the part of the above predictions which can be expected to become phenomenologically relevant 
first. These terms can be obtained from Eqs. (15.151) and (15.161) by removing all H-functions (the 
expansions of which at x = 1 generally start with (1— x) m with m > 1) and setting x = 1 in the 
remaining terms. For the convenience of the reader we finally write down the resulting expressions, 



d(3), 

qg 



x ) = hi 6 (l-jc) ■ 



+ ln 5 (l 



22 „ 3 
27 



27 



AF 



+ ln 4 (l-. 
81 

+ (ln 3 (l-x)) 



/293 


80 


V 27 


~ ~9~ 




116 




81 



27 
/4477 



AF n f 



CAFH f + \l62 
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c 
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AF 



CfH f 



81 CApn f 



(5.24) 



\n 6 (l-x) -0 



c, . 70 Q 14 9 9 4 9 



+ In 3 ( 1 -x) — Ci F C F - — C l AF C L F - — Ci F C F n f 

16 „ \ , „ /637 



+ ln* (1-x) 



27 
3280 



49^ 3 256^ 2 ^ 17 ^ ^ 2 32 ^ 2 
— CafC p — —C AF Cffij- + —CAFC F rif + —CafCfv-j 



+ (In 3 (1-x)) 



(5.25) 
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Figure 12: Large-A^ expansions of the NNLO and N 3 LO gluon-quark splitting functions Pq g (N) 

(3) 

and Pqg (N) , in the former case compared with the exact result of Ref . [7] . The curves ' 1 ' , '2' , etc 
are successively including the leading logarithms, ./V -1 In 4 /V (left) and iV -1 In 6 N (right), next-to- 
leading logarithms, A^ _1 In 3 N (left) and TV -1 In 5 N (right), etc. All log-enhanced terms are shown 
at NNLO. Only the double-logarithmic contributions obtained in this article are available at N 3 LO. 



0.2 



P gq } ( N ) (x 1/2000) 



0.03 




-0.2 



Figure 13: As Fig. [[21 but for the NNLO (left) and N 3 LO (right) quark-gluon splitting functions. 
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The only surprise in these results and Eqs. (15.151) and (15.161) is the absence, indicated explicitly in 
Eqs. (15.241 and (15.251) . of the expected leading ln 6 (l— x) terms. In general the - related, recall 
Eqs. (5.7/8) and Eqs. (5.11/12) - coefficients of the leading \n 2n (\-x) terms of P$ and 
receive contributions with both signs from Eq. (12.251) . which happen to cancel at n = 3. Such 
'accidental' cancellations of leading coefficients at some orders do occur, for a famous example see 
Refs. [49], and do not indicate any problem. We expect a non-vanishing coefficient of ln 8 (l— x) 
at five loops. Analogous comments apply to the (presumably related) absence of the expected 
ln 5 (l— x) contribution to the diagonal fourth-order splitting functions in Eqs. (15.191 ) and (15.201) . 

The Mellin moments of Eqs. (15.241) and (15.251) and their NNLO counterparts in Eqs. (4.18) 
and (4.19) of Ref. [7] are numerically illustrated in Figs. 12 and 13. Unlike analogous results 
for the non-singlet coefficient functions in Refs. [19,20], even all \\\N enhanced terms do not 

(2) 

provide a meaningful approximation at relevant values of N for Pg q , and the double-logarithmic 
contributions (blue curves in the electronic version) are nowhere near the full results in either case. 
Hence our predicted double-log terms are not phenomenologically useful without further four-loop 
information, e.g., in the form of an extension of Refs. [23] to the next order. Indeed, given the fact 
that a first fourth-order calculation of a splitting-function moment has already been presented [63], 
such an extension should become feasible in the foreseeable future. 



6 Discussion and outlook 

It is an old idea to replace the factorization- scheme dependent NLO (and now NNLO) quark and 
gluon distributions by suitable reference observables. One may expect to obtain smaller higher- 
order corrections to the coefficient functions for other observables in such a scheme, having ab- 
sorbed some universal corrections entering all coefficient functions in, e.g., the usual MS scheme. 
The predictions of (massless) perturbative QCD then depend on only one unphysical scale and 
scheme, for the renormalized strong coupling oc s . Moreover physical constraints on the chosen 
observables, for instance on their positivity or endpoint behaviour, would directly apply to the 
unavoidable non-perturbative initial distributions for the scale evolution in this case. 

An obvious quantity to replace the overall (flavour-singlet) quark density is the (singlet part of 
the) standard structure function F2 in (say, charged- current) deep-inelastic scattering, i.e., one can 
define physical quark densities q s by F2 = q s to NLO and NNLO. Any such scheme is usually 
called a DIS scheme [64]. Of course, this condition does not specify the lower row of the scheme- 
transformation matrix Z discussed in Section 2. The standard choice adopted in this respect, used 
in parton-distribution analyses mainly in the 1990s, was to arbitrarily use the coefficient functions 
for F2 also for this lower row. Consequently the gluon distribution (which, in practical analyses, 
would profit most from being defined in terms of an observable) remains entirely unphysical in 
this version of the DIS scheme, which is therefore mostly useful for quark-dominated observables 
such as Drell-Yan lepton-pair production in /^-collisions not too far from the threshold. 
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Various alternatives to the above arbitrary choice - which, given that gluon-dominated quanti- 
ties usually receive larger higher-order corrections than F2, 'under-transforms' the gluon density - 
have also been discussed already long ago, e.g., in Ref. [16]. The cleanest possibility is obviously 
the addition of a second observable which exactly mirrors the behaviour of F2, i.e., is directly 
sensitive only to the gluon distribution at LO. An example already mentioned in Ref. [16] is deep- 
inelastic scattering through the exchange of a scalar (]) interacting with hadrons only via a (]) G^G^ 
coupling to the gluon field- strength tensor. This interaction is realized in the Standard Model by 
the Higgs boson, for a very heavy top quark and negligible couplings to all other flavours [22]. 

In fact, as discussed in Section 2, the corresponding coefficient functions had to be calculated 
to order a s 2 in the approach of Refs. [23] to the determination of the complete set of third-order 
(NNLO) splitting functions, carried out for all-V/all-x in Ref. [7]. The numerous three-loop inte- 
grals for the latter were evaluated to order £° in dimensional regularization in order to obtain also 
the third-order coefficient functions for F2, Fi and F3 [8-10]. Hence also the three-loop coefficient 
functions for the ^-exchange structure function F§ could be calculated with a relatively minor extra 
effort. In Section 3 of this article we have presented and discussed these coefficient functions. 

While, of course, Higgs-exchange DIS is of no experimental relevance at any past or fore- 
seeable lepton-nucleon facility, these results have already proven useful in theoretical research as 
discussed above. Both the quark and gluon coefficient functions for F§ show, as expected, consider- 
ably larger higher-order corrections than their counterparts for F2. Their most prominent feature is 
the double-logarithmic enhancement of the 1 /^-contributions dominating their small-* behaviour, 
i.e., the occurrence of terms up to x~ 1 In 2 ' 1 x in the n-loop coefficient functions c^ n ] (x) . These terms 
arise from using the large-m t0 p effective interaction outside its region of validity, see also Ref. [53] 
for a recent discussion in the context of Higgs production in proton collisions. 

The above double-logarithmic contributions also render the physical evolution kernels for the 
system (F2,F&), presented and discussed in Section 4, unstable within the physically interesting 
region of x > 10~ 4 relevant to crucial high-scale processes at the LHC. Consequently this scheme 
is unattractive for practical analyses of hard processes. There are various possibilities, theoretical 
and phenomenological, to improve upon this situation which we, however, leave to future research. 
Already the present heavy-top results for the coefficient functions, however, can provide useful 
checks of, e.g., the order-dependence of fitted MS gluon densities by comparing the resulting 
shapes of the (in principle order-independent) quantity F§ outside the small-* region. 

The above flavour- singlet evolution-kernel matrices share a theoretically very interesting fea- 
ture with the (scalar) physical kernels for a rather wide range of non-singlet observables [19,20]: 
They show a single-logarithmic enhancement at large values of the scaling variable x, i.e., terms 
of the form (1 — x) a \n k (1—x) with k > n do not contribute to the N"LO physical kernels for any 
value of a. This structure is in contrast to that of all coefficient functions and (pure-) singlet splitting 
functions. We consider this feature, which is already established to N 3 LO for non-singlet structure 
functions, as a general property of the physical evolution kernels, and consequently conjecture that 
it holds to N 3 LO also for the present (F2,F$) singlet case. 
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Hence the double-logarithmic contributions from (a) the four-loop splitting functions and (b) 
the scheme transformation including lower orders and the three-loop coefficient functions have to 
cancel, i.e., the so far unknown corresponding terms in (a) can be inferred from our now com- 
plete results for (b). We have provided and discussed these predictions in Section 5. In particular, 
we note that the known double logarithms in the splitting functions to NNLO exhibit a particular 
colour structure in powers of (Ca—Cf) which is non-trivially fulfilled also by our N 3 LO predic- 
tions, a fact that we consider a significant check of their correctness. 

(3) 

As also discussed in Section 5, the predicted contributions to the splitting functions P-i (x) 
alone are clearly insufficient for use in phenomenological applications. However they will become 
very useful in connection with future fourth-order calculations, such as the extension of Ref. [28] to 
higher moments and the flavour- singlet sector. Already the corresponding three-loop predictions, 
had they been known at the time, would have led to more accurate x-space approximation on the 
basis of the moments calculated in Refs. [23] than those actually obtained in Refs. [44,65]. In fact, 
this situation clearly illustrates the general importance of calculations to at least NNLO: It usually 
requires two non-trivial orders to reliably identify most general structures - for another example 
recall the 'curious incident' of the coefficient of ln(l— x) in the three-loop diagonal splitting func- 
tions [6,7] which triggered Ref. [15] which in turn greatly assisted further developments including 
the determination of the NNLO splitting functions for the final-state fragmentation case [29,66]. 
We hope, accordingly, that also our present results will assist further developments in perturbative 
QCD, e.g., in the approach introduced in Ref. [60]. 

Form files of our main formulae, and Fortran subroutines of the ^-exchange coefficient 
functions and the resulting physical evolution kernels can be obtained from the preprint server 
|HTTP: //arXiv.orgl by downloading the source of this article. Furthermore they are available 
from the authors upon request. 
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Appendix A: The third-order coefficient functions for 



In this appendix we present the very lengthy exact expressions for the third-order scalar-exchange 

(3) (3) 

coefficient functions cl ^(x) and (x). The notation is the same as used for the first- and second- 
order results given in Eqs. (13.51) - ( 13.9b above. The result for the quark coefficient function reads 

c®(x) = C F CI ^ gq (x)(-1652/15H^ 2 2 +240H 1 ,_ 3 ,o- 160H lj _ 2 ^ 2 + 256H 1 ,_ 2 ,o,o 

- 704/3 Hi _ 2) -i,o + 128/3 Hi _ 2)2 + 712/3 H^Cs + 680/3 H^ok + 32H M) _ 2)0 

- 376/3 Hi.0,0,0,0 + 80/3 H M C 3 + 168H M , ^ - 248/3 Hi,i A o,o + 736/3 H u ^ 2 
-896/3H 1 j, l! o,o-1208/3H l!U j,o-400Hi i i,i iU -1208/3H 1 j, l!2 -392H lil , 2i o 

- 1232/3H 1 j, 2 a-288H lil)3 +776/3H l!2 C2-784/3H 1:2A o-384H li2)li o 

- 1216/3 Hi.2,1,1 - 1192/3 Hi A2 - 1024/3 Hi, 3 ,o- 1144/3 Hi )3 ,i -272H M ) 

+ Pgq(-x) (-988/15 H_iC 2 2 - 992/3 H_ 1; _3,o + 1712/3 H_ 1; _ 2 C 2 - 432H_ l5 _ 2 ,2 
+ 832/3 H_i _2,-i,o - 1688/3 H_i _ 2A o + 2240/3 H_i + 272H_i _i _ 2 ,o 
-2120/3H-i_i_iC2-944/3H_i_i_i_i,o + 2128/3H_i -1-1,0,0 
+ 1648/3 H-i-i -1,2 + 2464/3^1-1,0^2 -1832/3^1-1,0,0,0-27^-1 -1,2,0 

- 992/3 H_i _i )2) i - 2104/3 H_i _i, 3 - 536H_i, C 3 - 1 160/3 H_i A0 £ 2 

+ 232^1,0,0,0,0 + 64/3 H_ 1)2 C 2 + 304/3 H_i )2)0> + 224/3 H_i A i,o + 96H_ 1;2 , 2 
+ 256/3 H-1,2,1,1 + 632/3 H_i )3 ,o + 784/3 H_i )3 ,i + 1024/3 H_i, 4 ) 

- 35864/27 ( 1 - 8/4483/x 2 + 46675/53796/x + 7755/ 17932x 

+ 3883/13449 x 1 )H_i |0 + 32/3 (1 - 75/jc - 56x)H 2 ,i,i,i + 80/3 (1 - 162/5/jc 

- 267/10x)H_2,o A + 208/3 (1 - 48/13/x - 279/26x)H 5 + 112/3(1 - 24/7/x 

- 29/14x)H_ 3j2 - 208/3 (1 - 36/13/jc- 127/13 x)Ho A o£ 2 

+ 2368/15 (1 - 349/148/jc - 96/37 %)H C 2 - 6112/9 (1 - 841/382/* + 85/764 x 
+ 54/ 191 x 2 )H 2 ^ 2 - 2620/9 ( 1 - 1356/655/x + 1 87/262* + 44/655 x 2 )Hi£ 3 
+ 1 0040/8 1 ( 1 - 1 6057 / 1 0040/x + 24 1 9 1 / 1 0040 jc-865/251jc 2 )H u 

- 72928/27 ( 1 - 5560 1 /36464/x - 946 1 /36464 x - 1 697/9 1 1 6 x 2 ) ^ 3 

+ 59026/81 (1 -86857/59026/x + 24793/59026x-17892/29513jc 2 )Hi )0 
+ 7976/9 (1 -131 2/997 fx - 229/ 1994 x + 60/997 x 2 )Hi )0 ,o,o 

- 9376/9 ( 1 - 382/293/x + 55/586 x + 32/293 x 2 )Hi j C2 
-2648/3 (1 - 1289/993/jc- 1/331 x + 80/993 jc 2 )^^ 

+ 11420/9(1 -3704/2855/x + 257/2855x + 336/2855x 2 )H MA o 

+ 12092/9 ( 1 - 3614/3023/jc + 107/3023 x + 3 12/3023 x 2 )H M , 2 

+ 11200/9(1 -208/175/x-73/1400x + 3/35x 2 )Hi, 3 + 11240/9 (1 - 1664/ 1405 /* 
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+ 161/2810*+ 148/1405 * 2 )H M , M + 12076/9 (1 - 3562/3019/* + 103/3019* 

+ 312/3019* 2 )Hi !2J + 1336(1 - 195/167 /* + 6/167* + 52/501* 2 )H M)li0 

+ 1 3252/9 ( 1 - 3852/33 1 3 fx - 53/33 1 3 x + 320/33 1 3 x 2 )H lj2 ,o 

+ 1 289237 /243 ( 1 - 1478082 / 1 289237 fx - 205504 / 1289237 x 

+ 167284/1289237 * 2 )Hi- 68764/27 (1 -38843/34382/* -4969/34382* 

+ 432/17191 x 2 )Hi£ 2 + 1328/9 (1 - 91/83/jc — 89/166* + 8/83x 2 )Hi,_ 2 ,o 

+ 72760/27 (1 - 1 9543/1 8 1 90 fx - 419/1819 x + 9/9095 x 2 )Hi )0 ,o 

+ 1 1480/9 ( 1 - 1527/ 1435/* + 421/2870 x + 268/1435 x 2 )H 2 ,o J o 

+ 54620 /27 (1 - 2785 1 /273 10 fx -2458/1 3655 x - 146 /273 1 x 2 )Hi , 2 

+ 48764/27 (l-24351/24382/*-4361/24382*-782/12191* 2 )H Mil 

+ 56864/27(1 - 28103/28432/* -365/1777*- 383/7108 * 2 )H U;0 

+ 94 1 2/9 ( 1 - 2040/2353/jc + 1046/2353 jc + 800/2353 x 2 )H 4 

+ 1 3468/9 ( 1 - 2866/3367 fx + 557/3367 jc + 648/3367 x 2 )H 2)2 

+ 13736/9(1 - 1431/1717/* + 623/3434* + 320/1717* 2 )H 2!M 

+ 4904/3 (l-507/613/* + 247/1839* + 316/1839* 2 )H 2 ,i, 

- 9412/9 ( 1 - 1932/2353/* + 956/2353 x + 800/2353 x 2 )H ^ 2 

+ 5200/3 (1 -523/650/jc- 597/1300 *)^ 5 + 12472/9(1 - 1128/1559/x 

+ 51 l/1559x + 368/1559x 2 )H 3! o + 15968/9 (1 - 1095/1996/*+ 1303/3992* 

+ 1 02/499 x 2 )H 3 ,i +85748/27(1 - 10598/21437/*- 3388/21437* 

+ 70/1261 x 2 )H 2 , - 3896(1 - 1 888/4383/jc + 200/4383 x + 580/4383 x 2 )H C, 3 

+ 122288/27 (1 - 3219/7643/*- 6949/61 144* + 1455/15286* 2 )H 3 

- 122288/27 (1 - 12395/30572/*- 3535/61 144 x+ 1455/15286* 2 )H ^ 2 

- 160(1 -2/5 fx- 1/10jc)H 2 ,-2.o + 27832/9 (1 - 8245/20874/jc- 739/6958 x 
+ 39 1 / 10437 x 2 )H 2) i + 1 600/3 ( 1 - 9/25/* - 6 1 /200 *)H_ 3 A o 

+ 3047527/243 (1 - 1062210/3047527/jc - 139147/6095054 x 
+ 522652/3047527 * 2 )H + 170344/27 (1 - 13635/42586/* -4134/21293* 
+ 3017/21293x 2 )H , ,o + 4904/3(1 -166/613/x-129/613jc + 72/613x 2 )H_ 2 , 
+ 7796/3(1 - 270/1949/*+ 121/1949* + 216/1949* 2 )H_ 2 J; 2 + 480(1 -2/15/ 

- 1 1 / 10 jc)H ,o,oao - 3 1780/9 ( 1 - 206/ 1589/jc - 54/ 1 135 x + 928/7945 x 2 )H_ 2 
-381824/81 (1 - 180471/1527296/*- 2124809/4581888* + 250877/572736* 
+ 2384/3 (1 - 12/149/* -25/298*)H_ 4 , - 5344/3 (1 - 13/167/*+ 125/668* 
+ 18/167jc 2 )H_ 2! 9 +74912/27 (1 -481/18728/*- 1707/18728* 

+ 1 305/4682 x 2 )H_ 2 , - 25544/9 ( 1 + 54/3 193/* - 45/3 193 * + 400/3 1 93 * 2 )H 
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+ 832(1 + 1/6/* + 69/104 *)H ,o^ 3 +9464/9(1+ 512/1 183/* + 509/1 183* 
+ 40/1183* 2 )H_ 1 _ 1 , 2 -4060/3 (1 +226/435/*+ 13/35* + 136/3045 ^)H_i 

- 1 112/3(1 + 76/139/* + 347/278*)^ 3 + 2800/3 (1 + 4/7/jc+ 1 17/350x)H_ 2 -1,2 
+ 4732/3 (l+748/1183/x + 789/2366x + 76/1183x 2 )H_iC 3 

-3488/3(1 +72/109/* + 85/218*)H_ 2 ,_i£ 2 -7528/9 (1+648/941/* 

- 167/1882* - 268/941 * 2 )H , ,o,o - 16780/9 (1 + 3232/4195/* 

+ 1274/4195* + 336/4195 x 2 )H_ u + 4648/3 (1 +65/83/* + 92/581* 
+ 8/83* 2 )H_i,_i, ,o + 3176/3 (1 +314/397 /* + 389/794 *)H_ 2 £ 3 

- 7984/9 ( 1 + 396/499/* + 3 1 9/998 * + 44/499 * 2 )H_ i , 2 ,i 

+ 20036/9(1 +3982/5009/*+ 1414/5009* + 416/5009* 2 )H_i,oC2 

+ 186364/81 (l+37486/46591/*-6167/186364* + 22118/46591* 2 )C 2 

-5432/9(1 +558/679/*+ 113/679* + 8/97* 2 )H_i ; _i_i, 

- 7240/9 ( 1 + 752/905/* + 283/905 * + 72/905 * 2 )H_ i , 2 , 

+ 576(1 +71/81/*+ 17/216* + 8/81* 2 )H_i,_ 2 , + 144(1 +8/9/* 

+ 37/54*)H_3_i ; o- 19328/9 (l+601/604/* + 751/4832* + 63/604 * 2 )H_i A o,o 

-424(1 + 160/159/* + 427/318*)H 3i0 ,o- 1376/3 (1 +44/43/* 

+ 53/86x)H_2-i-i ) o- 293866/81 (1 + 151 137/146933/*: + 8738/146933*: 

+ 73960/146933 * 2 )H ,o + 59324/27 ( 1 + 15328/1483 1 /* + 2398/1483 1 * 

+ 4676/14831* 2 )H_ 1A0 + 24536/27(l +3303/3067 /* + 385/6134* 

+ 1 178/3067 * 2 )H_i, 2 - 704/3 (1 + 25/22/* + 43/44 *)H 2 , ,o,o 

-896/3 (1 +8/7/*+ ll/14jc)H_ 2l 2,i -2024/3 (1 + 292/253/* + 37/46 *)H_ 2 , 3 

-38680/27 (1+6051/4835/* + 823/19340*+ 1759/4835* 2 )H_i^ 2 

-656/3 (1 + 52/41/* + 75/82*)H_ 2 , 2 , - 186364/81 (1 + 121647/93182/* 

+ 40363/186364* + 22118/46591* 2 )H 2 + 2144/3(l +88/67/*+ 125/134*)H_ 2 , ^ 2 

-1336/3(1 +224/167 /* + 415/334*)H 2 j,o !0 + 240(l + 8/5/*+ ll/10*)H_ 2 _ 2j0 

+ 1912/3 (1 + 348/239/x + 485/478 Jc)H_ 2 -i,o,o + 408(l +224/153/* 

+ 509/306*)H 3 ^2-28288/27(l+687/442/* + 53/7072* + 581/1768* 2 )H_i,_i, 

+ 512/3 (1 +53/32/* + 93/64 *)H 2 C 3 - 1168/3 (1 + 124/73/x+271/146x)H 3 ,i, 

+ 880/3 (1 + 102/55/*+ 171/1 10*)H 2 ,i£ 2 -336(1 + 124/63/*+ 13/6*)H 3 , 2 

-632/3(1 + 168/79/* + 547/158*)H 4 , - 896/3(1 + 18/7/* + 37/14*)H 3 ,u 

-616/3 (l+204/77/* + 629/154*)H 4 ,! +368/3(1 +82/23/*+ 119/46*)H 2 ,oC 2 
-488/3 (1 +260/61/* + 421 /122*)H 2 , 2 ,o+ 104/3 (1 +72/13/* + 95/26*)H_ 3 £ 2 

- 9986/45 ( 1 + 27874/4993/* + 1 529/4993 * - 328/4993 * 2 ) £f 
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- 248/3 ( 1 + 224/3 1 fx + 331/62x)H 2 , 3 -64(1+ 73/6/x + 28/3 x)H 2 ,i,2 

- 160/3 (l + 73/5/x+56/5x)H 2 ,2,i - 128/3 (1 +73/4/*+ 14x)H 2 ,i,i,o) 

+ C|C A (p gq W(-1856/15HiC|- 176Hi,_3,o + 264Hi i _ 2 ^2 + 1 12Hi,_ 2 ,_i, 

- 640/3 Hi _2,o,o - 208Hi _2,2 - 32H lj0 ^ + 640/3 Hi, ,o£2 - 544/3 Hi, ,o,o,o 
+ 632/3 H M C 3 - 464/3 Hi.i _ 2| o + 1400/3 Hi,i,oC 2 - 928/3 Hi,i A o,o 

+ 784/3^,1,1^2-320/3^,1,1,0,0-40^,1,1,1,0-104^,1,1,2-120^,1,2,0 

- 304/3 Hi,i )2) i - 432H M)3 + 832/3 H^ - 592/3 Hi, 2j0 , - 368/3 U l>2 ,i,o 

- 80Hi )2 ,i,i - 488/3 Hi, 2 ,2 - 176Hi )3 , - 440/3 Hi )3) i - 616/3 Hi, 4 ) 

+ p gq (-x)(1016/15H_iC| + 192H_i _ 3) o - 344/3 H_i _ 2 £ 2 - 688/3 H_i _ 2 -1,0 
+ 488/3 H_i -2,0,0 - 136H_i ,_iC 3 - 224H_i _i _2,o + 296/3 H_i _i _iC 2 
+ 944/3 H_i,_i,_i,_i,o- 736/3 H_i,_i,_i, ,o + 176/3 H_i,_i,_i, 2 

- 112/3 H_i _i,oC2 + 80H_i _i A o,o + 1 12/3 H_i _i |2 ,o + 160/3 H_i,_i, 2 ,i 

+ 8/3 H_i ; _i, 3 - 448/3 H_i,oC 3 - 320/3 H_ 1 ,o,oC2 + 96H_ 1 , ,o,o,o- 232/3 H_i, 2 C2 
+ 368/3 H_i j2 , ,o + 128/3 H_i, 2 ,i,o + 32/3 H_i, 2 ,i,i + 64/3 H_ 1;2 , 2 + 184/3 H_i i3 ,o 
+ 64H_uj + 344/3 H_i, 4 ) + 1280(1 - 7/1080/x 2 + 139/640/x + 15683/17280* 
-7/240 ;c 2 )H_ 1 ,o + 742/27 (1 - 13896/371/jc + 59/2 jc - 144/53x 2 )^ 3 
+ 904/9(1 -543/1 13/x + 853/226x- 64/1 13x 2 )H ^ 3 - 128/3 (1 - 15/8/x 

- l/4x)H_ 2 ,o,o,o - 27938/27 (1 - 22578/ 13969/jc - 35941/27938x 

- 720/ 1 3969 x 2 )H ^ 2 + 27938/27 ( 1 - 2 1 1 20/ 1 3969/x - 2778 1 /27938 x 
-720/13969x 2 )H 3 + 5032/9(l-936/629/x-1715/1258x+12/629x 2 )H 2 ,i,i 
+ 176(1 -16/ll/x-107/66x)H 3j i, + 1120/3 (l-46/35/x-31/35x)H 2j i A0 
+ 34460/27 (1 - 10674/8615/x- 1888/8615x-72/8615x 2 )H 2 , 

+ 784/3 (1 -8/7/x- 13/14x)H 2i2! o + 41896/27(l -5901/5237/x 

- 4243/20948x + 288/5237 x 2 )H 2 ,i + 688/3 (1 - 48/43/x - 125/86x)H 3i2 

+ 43340/27 ( 1 - 4745/4334/x - 3007/10835 x + 24/2 1 67 x 2 )Hi ,2 

+ 38416/27 (l-20597/19208/x-961/4802x + 30/2401x 2 )H Mi0 

+ 43012/27(1 -23041/21506/x-2357/10753x + 544/10753x 2 )H Mil 

+ 6128/9(1 -807/766/x-605/766x + 30/383 x 2 )H 2 j, 

+ 78766/27 ( 1 - 8207 1 /78766/x - 2 1 843/78766 x + 20/39383 x 2 )Hi ,0 

+ 7304/9(1 -927/913/x- 1483/1826x + 60/913x 2 )H 2 , 2 + 736/3 (1 - 1/x 

-35/46x)H 2 ,i, 1,0 + 8728/9 (1 - 1073/1091/x- 436/1091 x + 24/1091x 2 )Hi i2 ,o 

+ 274384/81 (1 -33685/34298/x-20807/68596x-331/34298x 2 )H M 

+ 10688/9(1 -653/668/x-503/1336x + 3/167x 2 )H U! o, + 768(l-23/24/x 
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- 143/288*+ l/72* 2 )H M) i )0 + 1923029/486 (1 - 1747009/ 1923029/* 
-795761/3846058x-40576/1923029x 2 )H 1 +9224/9(1 -1003/1153/* 

- 1105/2306 jc+ 12/1 153 jc 2 )Hi j2j1 -286819/81 (1 - 19056/22063/* 

- 181409/286819* + 360/286819* 2 )^ 2 + 1262/27(1 -540/631/* 

- 17713/1262jc - 1600/631 * 2 )H ,o,o + 7528/9 (1 - 796/941/jc - 583/941 * 

- 4/941 x 2 )Hij,ij +8552/9(1 - 898/1069/*- 1039/2138* + 12/1069 * 2 )H M , 2 
+ 8792/9 (1 - 894/1099/*- 1945/2198* + 132/1099 * 2 )H 3 ,i 

+ 848/3 (1 -42/53/*- 97/106*)H 3jM +286819/81 (1 -225210/286819/* 

- 873 1 1 /2868 19 x + 360/2868 19 x 2 )H 2 - 52040/27 ( 1 - 4073/5204/* 
-4477/13010* + 36/6505 * 2 )Hi£ 2 -480(1 -7/9/x- 31/60 *)H 2 ^ 3 
-5800/9(1 -112/145/*+98/145*- 16/725 * 2 )H_i _i,o 

+ 1023301/486(1 -789822/1023301/* -203321/1023301* 

- 6016/1023301 jc 2 )H - 13160/9(1 - 181/235/*- 1627/3290* 
+ 24/1645* 2 )H M C 2 + 13972/9(l -2686/3493/*- 1709/6986* 
- 32/3493 * 2 )Hi A o+ 1088/3 (1 - 13/17/jc-43/68jc)H 2 , 2 ,i 

+ 1216/3 (1 -29/38/*- ll/19*)H 2 j, 2 - H56(l - 1916/2601/*- 281/578* 
+ 80/2601 * 2 )Hi^ 3 - 2368/3 (1 - 47/74/* - 65/148*)H 2 ,^ 2 

+ 3008/3 (1 -467/752/*- 293/752*+ 1/47 * 2 )Hi,o,o,o 

+ 12904/9(1 -995/1613/*- 1553/3226* + 24/1613* 2 )H 1;3 

+ 4472/9 (l-1719/2795/*-4153/11180* + 348/2795* 2 )^ 2 2 

-10856/9(1 -825/1357/*- 1729/2714* + 96/1357* 2 )H 2 £ 2 

+ 640/3 (1-3/5/* -2/5*^-2-1,2+ 1472/3(1 -55/92/*- 13/23 *)H 2 ,o, ,o 

-13864/9(1 -1028/1733/*- 1667/3466* + 36/1733* 2 )H lj0 C2 

+ 10120/9(1 -732/1265/*- 1469/2530*+ 144/1265 * 2 )H 3)0 

-528(1 -6/ll/*-68/99*)H 3 ^ 2 + 304(l- 10/19/*- 17/38jc)H 2 ,i,i,i 

+ 1088/3 (1 - 71/136/* - l/2*)Hi _2,o + 8968/9 (1 - 567/1 121/* - 989/2242* 

+ 72/1 121 * 2 )H 2 ,o,o+ 1360/3 (1 -42/85/*- 117/170*)H 3) o,o 

+ 23 1 5033/324 (1 -1407 1 5/3307 19/* - 306555/23 1 5033 * - 7 104/23 15033 * 2 ) 

- 2560/3 (1 - 17/40/* - 31/80*)H 2 ^2 + 2432/3 (1 - 31/76/* - 29/76*)H 2 , 3 

- 14632/9 ( 1 - 24/59/* - 2 1 59/3658 * + 240/ 1 829 * 2 )H ^ 2 

+ 1 184/3 (1 - 15/37/* - 19/74 *)H 2 _ 2 ,o + 1472/3 (1 - 37/92/* - 75/184 *)^ 3 

+ 14632/9(1 -690/1829/*- 1919/3658* + 240/1829* 2 )H 4 

+ 800/3 (1 - 9/25/*- 4/25 *)H_ 2 ,oC2+ 127235/81 (1 -44982/127235/* 
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- 301739/254470 x - 3352/127235 x 2 )H ,o + 304/3 (1 - 6/19/x + 26/ 19x)H , ^3 

- 8084/9 (1-51 4/202 1/jc+l 900/202 1 x - 48/202 1 x 2 ) H x & 

-2368/3 (1 - 17/74/x + 41/148x + 3/37jc 2 )H_ 2 -1,0 + 2848/15 (1 — 39/178/jc 

- 407/7 1 2 x)H £ 2 - 928/3 ( 1 - 6/29/x - 2/29 x)H- 2 ,3 + 7556/9 (1 -3 1 8/ 1 889/* 
+ 1750/1889x- 32/1889 jc 2 )H_i,o,o- 1120(1 - 4/105 /jc-11/420x 

+ 3/35x 2 )H_2^2 + 2704/3(l-7/338/x+12/169x + 20/169x 2 )H_ 2 ,o,o 
+ 576(1 +23/648/x + 235/216x-2/81x 2 )H_ li2 + 2176/3 (1 +9/136/* 

- 13/68x + 3/34x 2 )H_ 2 ,2 + 464(l +2/29/jc+ 13/174 jc)H_ 3;0 ,o 

- 1312/3 (1 +3/41/*+ 19/164 x)H_3^ 2 + 1904/3 (1 +9/119/x + 20/119x 
+ 12/1 19jc 2 )H_ 3 ,o- 1792/3 (1 +3/28/jc+ 11/1 12x)H_ 3 ,_i, 

- 1 120/3 ( 1 + 6/35/x + 33/35 x)H 5 + 1 120/3 ( 1 + 6/35/x + 34/35 x)H , ,o^2 

+ 4576/9(1 +243/1 144/jc + 85/143jc-4/143jc 2 )H_ 2 ,o + 16(1 + 1/3/jc-13/6jc 

+ 4/3x 2 )H_ 1 ,2,i -824(1 + 109/309/x + 235/618jc + 8/103 jc 2 )H_!_ 1A o 
-576(1 +31/72/x + 53/108x+ 1/27 x 2 )H_i,_2,o + 1024/3 (1 + 111/256/x 

+ 97/256jc + 3/16jc 2 )H_i,oao- 480(1 +7/15/x + 3/10x)H_2,-2,o-64H_ 2 ,2,i 
-432(1 +38/81/x + 31/162x + 4/81x 2 )H_ 1 ,_ l! 2-312(l +56/117/*+ l/6x)£ 5 
+ 768(1 + l/2/x+7/24x)H_ 2 ,-i,-i,o + 944(1 + 199/354/* + 43/1 18 x 
+ 2/59x 2 )H_i_iC 2 - 836(1 + 122/209/x+ 149/418^ + 32/627 x 2 )H_i^ 3 

-464(1 + 18/29/x + 37/174x+ 14/87 x 2 )H_^ ^2 -240(1 +83/90x)H OA oao 
+ 1024(1 + 41 /64/jc + 49/96 x+ 1/48 ^)H_i _i _i )0 + 1072/3 (1 + 97/134/* 
+ 9/67x+ 12/67 x 2 )H_i, 3 + 272/3 (1 + 13/17/x + 2/17x + 4/17x 2 )H_i, 2 , 

- 1024/3 (1 + 15/16/x+ 35/64x)H_2,-i,o,o - 320/3 ( 1 + 6/5/x + 59/20 jc)H 4 ,o 

- 160(1 +9/5/x + 21/20jc)H_ 2 C3 +512/3 (1 + 15/8/x + 37/32jc)H_ 2 ,-iC2 
-64(1+3/x+67/12x)H 4; i + 13952/9 (1 - 337/872*+ 13/109jt 2 )H ,o,o,o 

- 160/3 (1 - l/10x)H_ 2) 2,o + 1952/3 (1 + 1/61 +)H_ 4j0 + 416/3 (1 + 2/13x)H_ 3 , 2 ) 

+ C^p gq W(1288/15H 1 C2 2 +640/3Hi ! _3,o-304Hi i _2^2-160Hi,_2,-i,o 
+ 256H 1) _ 2 ,o,o + 224H l5 _ 2 ,2 + 400/3 H^Cs - 96H 1A o^2 + 64H 1 , ,o,o,o 

- 352/3 H M C 3 + 416/3 H M ,_ 2 ,o - 944/3 H M) <£2 + 96Hi,i, ,o,o - 368/3 H Ui ^ 2 
-608/3^,1,1,0,0-400/3^,1,1,1,0-80^,1,1,1,1-208/3^,1,1,2-128^1,2,0 

- 64Hi,i, 2 ,i + 208Hi,i, 3 - 96Hi, 2 C2 - 272/3 Hi, 2)0 ,o - 400/3 Hi, 2 ,i,o 

- 272/3 Hi, 2 ,i,i - 80Hi,2,2 - 64Hi, 3 ,o - 80H Ui i + 64Hi, 4 ) 

+ ^(-^(736/15 H_^ 2 2 -832/3 H_ 1 ,_ 3> o + 400H_ 1 ,_ 2 ^ 2 + 352H_ 1 ,-2,-1,0 
- 1504/3 H_i_2,o,o-224H_i,_2,2 + 464H_i,_iC3 + 352H_i,_i,_ 2 ,o 
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-448H_i_i _i^2 -384H_i_i_i_i ) o + 704H_i_i_i,o,o + 256H_i _i_i,2 
+ 1472/3 H_i _i,oC2- 1552/3 H_i _i A o,o -128H_i -1,2,0 -128H_i _i A i 

- 384H_i,_i )3 - H2H_i,oC3 - 240H_ 1,0,0^2 + 448/3 H_i,o,o,o,o + 80H_i, 2 £ 2 

- 96H_i,2,o,o + 64H_i i3j0 + 64H_i, 3 ,i + 208H_ M ) - 6752/3 (1 - 2/633/x 2 
+ 1 83/844/x + 4349/5064 x)H_ 1)0 - 88/3 (1 - 210/1 1/x - 21/2 x) £ 5 
+424/15 (l-519/53/x+4/53x)£f +284/3 (l-581/284/x+887/284x)Hi,i 

+ 1088/3 (1 - 171/136/x-65/136x)Hi, 2 ,i + 920/3 (1 -57/46/x-67/115x)H 2 ,i 
+ 976/3 (1 -75/61/x- 101/244x)Hi,i,i,i + 1216/3 (1 -93/76/x 
-71/152x^1,1,2+ 1340/3(1- 358/335/x- 88/335 x)Hi,i,i 
+ 1856/3 (1 - 30/29/x - 27/58x)Hi,i,i, + 1 136/3 (1 - 72/71/x - 44/71 x)H 2;2 
+ 1436/3 (1 -364/359/x- 23/7 18 x)Hi, 2 + 2036/3 (1 -504/509/x 

- 193/1018jc)H m , + 1424/3 (1 - 84/89/x- 83/178x)Hi, 2 ,o 

+ 1600/3(1-9/10/*- 131/200x)H 2 ,i, - 1664/3 (1 -93/104/x 

- 25/52x)Hi ,_ 2 ,o + 1472/3 (1 - 81/92/x - 187/368x)H 2 ,i,i 

+ 916/3 (1 - 180/229/x- 167/458x)H 2)0 - 3265/6 (1 -493/653/x 
-4527/6530x)H! +896/3 (1-3/4/x- l/2x)H 2 , 2) i +784/3 (1 -36/49/x 

- l/2x)H 2 ,i,i,i +800/3 (1 - 18/25/x-27/50x)H 2 ,i,2+ 1232/3 (1 - 54/77/x 

- 34/77 x)H 3 ,i + 1216/3(1 - 12/19/x- 9/19x)H 2 ,i,i,o + 240(1 -3/5/x 

- 3/5x)H 3 , + 1024/3 (1 - 9/16/x- 15/32x)H 2 , 2 , - 70927/36 (1 - 2098/3733/x 

- 11754/70927 x) + 1088/3 (1 -9/17/x- l/2x)H 3 , 2 - 1088/3 (1 -9/17/x 

- ll/34x)H 2 -2,o + 368(1 - 12/23/x- l/2x)H 3 ,i,i +424(1 -55/106/x 

+ 91/212x)Hi )0 + 1312/3 (1 - 18/41/x- l/2x)H 3 ,i, + 1232/3 (1 -24/77/x 

- 13/22x)Hi,iC 2 +448(1 - 2/7/x- 3/7x)H 4 ,i - 596(1 - 36/149/x 
+271/447x)HoC3 + 512/3(l-15/64/x-77/128x)Hi,i, 0) o + 832/3(l-3/13/x 
- 3 /26x)H 2 , 1^2 — 1000/3 (1 -27/125/x- 189/250x)H l!0 ,o + 568(1 - 15/71/x 

- 121/213x)H^ 3 + 336(1 -4/21/x- 17/42x)H 4 , + 880/3 (1 -3/22/x 

+ 689/220 x)H 3 - 1592/3 (1 - 24/199/x- 55/398x)H ,oC3 + 976/3 (1 - 6/61 /x 
-33/122x)H 2 ,i, 0) o + 2060/3(l-32/515/x-861/1030x)HiC2-288(l-l/18/x 

- 19/36x)Hi A o,o + 4352/3 (1 + 3/68/x+ 15/272x)H_ 3 ^ 2 - 4192/3 (1 + 6/131/x 
+ 21/262x)H_ 3 , , - 1304(1 + 19/326/x+ 1607/978x)^ 3 - 3584/3 (1 + 9/112/* 
+ 33/224 x)H_ 3 ,o- 1056(1 + l/ll/x + 5/33x)H_ 2 ,2 + 3712/3 (1+3/29/x 

+ 15/1 16x)H_ 3 ,_i,o - 7120/3 (1 + 99/890/x+ 121/445 x)H_ 2 ,o + 845(1 + 8/65/x 
+ 39127/15210x)H ,o + 1424/3 (1 + 12/89/x- 9/178x)H 2 ,o^ 2 
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+ 5248/3 (1 +27/164/* + 167/656*)H_ 2 £2 -736(1 + 4/23/* + 9/46 *)H_ 2 , 3 
-2720/3 (1 + 3/17/* + 41/170*)H_ 2 ,o,o,o + 176(1 + 2/ll/*-5/22*)H 2 C 3 
-5392/3 (1 +63/337/* + 217/674*)H_ 2 ,o,o + 960(l + l/5/* + 2/9*)H_ 2 ,o^2 
+ 1648/3 (1 + 24/103/* - 273/412jc)Hi,oC2 - 1264(1 + 19/79/* + 62/79 *)H_i, 2 
+ 4160/3 (1 + 18/65/*+ 107/260x)H_ 2 -i,o + 168(1 + 1/3/*+ 1079/504*)H 2 
+ 4576/3 (1 +48/143/* + 81/286*)H_ 2 ,_i,o,o + 7288/3 (1+312/911/* 
+ 593/911*)H_^ 2 -2528(1 +85/237/*+ 151/237*)H_ 1A o+ 1152(1 +7/18/* 
+ 7/24 *)H_ 2 ^3 +2432/3 (1 + 15/38/x+6/19x)H_ 2 ,- 2 ,o + 640(1 +2/5/* 
+ 3/10x)H_ 2) _ 1)2 -224(l + 3/7/*+4/7x)(H_i ) 2,o + H_i )2) i) 

- 1184(1 + 16/37/* + 23/74*)H_ 2 ,-i^2 + 6992/3 (1 + 198/437/* 
+ 221/437 jc)H_i _i,o- 1088(1 + 8/17/*+ 11/34 jc)H_ 2 ,-i -1,0 

- 1088(1 +65/136/* + 55/102*)H_i i0 , ,o + 3472/3 (1 + 108/217/* 

+ 249/434x)H_i, C 2 + 256(1 + 1/2/*+ 19/48*)H 3 £ 2 - 864(1 + 1/2/* 
+ 7/12*)H_i i3 - 752/3 (1 + 24/47/*+ 17/94 *)H 2 , 3 + 368/3 (1 + 12/23/* 

- 5/46 *) ££3 + 944/3 ( 1 + 36/59/* - 19/118 *)H 2 £ 2 + 5984/3 (1+21 /34/x 
+ 395/748jc)H_i_i A o + 928(1 + 18/29/* + 31/58*)H_i _ 1)2 

+ 1640(1 + 138/205 /* + 221/410*)H_i^ 3 + 3424/3 (1 + 147/214/* 
+55/107x)H_i -2,0- 1744(1 +75/109/*+58/109x)H_i -1^2 -832H_ 3 , 2 

- 880/3 (1 +42/55/*+ 1173/220 *)H ^ 2 - 1632(1 + 13/17/* 

+ 9/17x)H_i ,-1,-1,0- 256(1 + 17/16/x-51/64x)Hi )3 - 236/3(1 +72/59/x 
+ 365/59 *)H ,o^2 -208/3 (1 + 18/13/*+ 105/52*)H 2 , , - 272/5 (1+40/17/* 
+ 2x)H £! - 168(1 +68/21/* + 20633/1512*)^ 2 + 2437/18 (1 +7972/2437/* 
+ 5058/2437 *)H - 160/3 (1 - 1 47/20 *)H , ,o,o + 328(1 -21/82*)H 5 
-416/3 (l-5/26x)H 2 ,o,o,o-328(l-13/82x)Ho,o,oC2-2816/3(l+3/88x)H_4, 

- 192(1 + l/6*)(H_2,2,o + H_ 2 ,2,i) + 120(1 +43/90 *)H ,o, ,o,o 

+ 1 172/3 ( 1 + 950/293 *)H ,o,o + 236/3 ( 1 + 233/59 *)H 4 - 64( 1 /* + 3/2*)H 3 A o) 

+ C f C A n / (p gq (x)(100/9HiC3 + 160/9^^2,0- 184/3 Hi i0 C2 + 688/9Hi, ,o,o 
-560/9^,^2 + 544/9^,1,0,0 + 128/3^,1,1,0 + 376/9^,1,1,1+448/9^,1,2 
+ 512/9Hi,2,o+416/9Hi,2,i+560/9Hi, 3 )+Pg q (-x)(-20H_iC3-16H_i,_2,o 
+ 80/3 H_i,_iC 2 + 224/9 H_i,_i,_i,o-32H_i,_i,o,o- 128/9 H_i,_i,2 
+ 1 12/9 H_i,oC 2 - 256/9 H_i, ,o,o- 128/9 H_i, 2 ,o- 80/9 H-1,2,1- 104/9 H_i, 3 ) 
-5744/27 (l+2/359/* 2 + 1751/4308/*+ 1049/1436*+ 113/1077 * 2 )H_i, 

- 64/9 (1 - 29/2/*- 23/2*)H 2 , 2 - 176/9 (1 - 63/1 1/*- 103/22*)H 3jl 
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- 176/9 (1 - 5/x-7/2x)H 2 ,i,i - 15527/27 (1 - 398483/139743 fx 

- 171 14/139743x- 512/46581 x 2 ) - 128/3 (1 - ll/4/x-5/3x)H 2 ,i, 
+ 448/9 (1 - 5/2/x - 7/4x)H 2 C2 - 5 12/9 (1 - 5/2/x - 53/32x)H 2 ,o, 
-202498/243(1 - 183399/101249/x- 84161/202498* - 17064/101249 x 2 )H 
+ 80/3 (1 -5/3/jc + 7/15x)H ^ 3 + 7376/27(1 - 1325/922/jc- 1517/1844x 

+ 68/46 1 x 2 )H(^2 - 7376/27 ( 1 - 627/46 1 jx - 1 27 1 / 1 844* + 68/46 1 x 2 )H 3 

- 1 10708/243 (1 - 133181/1 10708/x - 26707/1 10708 x - 646%/ 21 611 x 2 )Hi 
+ 11128/27 (1 -3301/2782/jc- 1556/1391 x- 10/1391 x 2 )^ 

- 12256/27(1 - 1761/1532/x-2001/3064x + 59/383x 2 )H , ,o 
-2560/9(1 -533/480/jc- 2/5 jc + 7/80jc 2 )H 2 ,i -45940/81 (1 - 24377/22970/jc 

- 5353 /22970x - 36/ 1 1485 x 2 )Hi )0 - 9 1 52/27 ( 1 - 2407 /2288 fx - 661 /2288 x 
+ l/22x 2 )Hi ;2 - 9728/27 (1 - 1261/1216/x - 467/ 1216x + 29/304 x 2 )H 2 ,o 
-34412/81 (1 - 1273/1229/jc- 1474/8603 jc+ 164/8603 x 2 )H M 
-86756/81 (l-22113/21689/x-7885/21689x+1064/21689x 2 )H ,o 

- 13 120/27 ( 1 - 833/820/x - 577/ 1640x + 7/205 x 2 )Hi )0 ,o 
+ 10624/27 (1 - 84/83/x- 899/2656* + 29/664 x 2 )HiC 2 

- 8264/27 (1 - 2085/2066/x- 643/2066* + 44/ 1033 * 2 )H M ,i - 64/3 (1-1/* 
-*)H- 2 ,2 -2704/45 (1 - 168/169/jc— 185/338x)£f - 9776/27 (1 - 2391/2444/* 

- 823/2444* + 22/61 1 x 2 )Hi,i, + 59024/81 ( 1 - 3207/3689/* - 8753/29512* 
+ 303/3689 x 2 )£ 2 - 59024/81 (1 - 11077/14756/jc-2459/29512jc 

+ 303/3689 * 2 )H 2 - 640/9 (1 - 3/10/*- 3/10*)H_ 3j0 - 320/9 (1 - 1/5/x 
-3/5x)H_2,o,o + 1312/27(l-l/8/x + 83/82x + 9/82x 2 )H_i^2 

- 1 664/27 ( 1 - 7/208/x + 79/ 104* + 1 /4* 2 )H_ i , ,o - 3392/27 (1+71 /424/x 
+ 123/424* + 13/106 * 2 )H_ 2 ,o + 256/3 (1 + l/2/x+ l/8x)H_ 2 ,-i,o 

+ 1072/9 ( 1 + 42/67 jx + 143/ 134x)H , ,o,o + 2944/27 ( 1 + 28 1 /368/x + 29/46* 
+ 3 /92 x 2 )H_ i _ i ,o + 1 60/27 (1 + 161 /20/x - 5/2 x - 3 /5 * 2 ) H i , 2 

- 80/9 (1 +54/5/x+ 1 13/10x)H ,oC2 + 80/9 (1 + 66/5/x+ 137/10x)H 4 
+ 32/9 ( 1 + 33/x + 55/2*)H 3;0 + 64(1 - l/4x)H_ 2 £ 2 ) 

+ C 2 n / (p gq (x)(-128/3HiC 3 -64/3Hi^ -8/9Hi^ 2 -16/3Hi A o,o 
+ 104/9^,1^2 + 88/9^,1,0,0 + 56/3^,1,1,0+104/9^,1,1,1 + 88/9^,1,2 
+ 272/9Hi,2,o + 184/9Hi,2,i + 104/9Hi, 3 )+Pg q (-x)(112/3H_iC3 
+ 128/3 H_i _ 2 ,o - 128/3 H_i,_^ 2 - 128/3 H_i _i _i,o + 224/3 H_i _i >0)0 
+ 64/3 H_i _i )2 + 64/3 H_i )0 £2 - 32H_i, ,o,o - 32/3 H_i, 3 ) 
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+ (2 -x)(-56/5 HoCl + 152/3 n ,o& + 88H OA oC2 - 120H 0) o,o,o,o + 16H 3 £ 2 

- 16H 3 , ,o - 16H 3il! o - 16H 3 ,i,i - 16H 3 , 2 - 48H 4 , - 48H 4 ,i - 88H 5 ) 

+ (2 - x) (24^ 5 - 32C2C3 ) - 1 1 2/9 ( 1 - 25/1 /x - 2/7 x) H 2 ,o,o + 64/3 (1 - 3 1 / 12 fx 
+ 13/4x- l/3x 2 )H_i, 2 + 28030/81 (1 - 16948/14015/x + 2234/14015x 

- 1488/14015jc 2 )^ 2 -544/9(1 -39/34/x-7/68x)H 2 ,i,i 
-28030/81 (1 -3168/2803 /*- 9646/14015*- 1488/14015x 2 )H 2 
-608/9(1 -21/19/x- ll/76x)H 2 j,o-704/9(l - 12/11/x- 17/88x)H 2 , 2 

+ 1072/9 (1 -72/67/x- 53/268 x)H ,oC2 - 848/9 (1 - 54/53/jc+ 19/106jc)H 3 j 

- 880/9 (1 -54/55/*+ 17/ 110x)H 3i0 + 4640/27 (1 - 142/145/x- 101/232* 

+ 9/145x 2 )Hi^ 2 - 1072/9 (1 - 60/67 /jc+ 19/268x)H 4 -7712/27 (1 - 823/964/jc 

- 1087/1928* + 6/241 x 2 )Hi j2 - 1624/9(1 - 139/174/jc- 95/203* 

+ 32/609jc 2 )Hi A0 -6448/27(1 - 1279/1612/jc-445/806x+ 12/403 jc 2 )Hi iM 
-7552/27 (l-185/236/jc- 1007/1888x + 3/118x 2 )Hij, 
-9148/27 (1 - 1790/2287/x- 1323/2287*- 192/2287 * 2 )Hi,o 

- 2465 17/243 (1 -1 85803/2465 17 fx - 1 8 1 6 1 5/493034 x - 5280/2465 17 x 2 )Hi 

- 28300/81 (1 - 5272/7075/*- 907/1415*- 576/7075* 2 )H M 

+ 2528/9(1 -43/79/*- 247/632*)H(£ 3 - 19352/81 (1 -2481/4838/* 
+ 19385/19352*- 392/2419* 2 )H ,o - 585697/486 (1 - 842416/1757091/* 

- 129025/585697* - 9872/251013* 2 ) - 12872/27 (1-705/1609/* 
-llll/3218* + 72/1609* 2 )H 2 ,o- 12832/27 (1 -351/802/*- 139/401* 
+ 18/401 x 2 )H 2 ,i - 440/9 ( 1 - 24/55/* + 89/1 10*)H ,o,o,o 
-254765/243 (1-14591/36395/*- 103486/254765 *- 2592/254765 * 2 )H 
+ 18752/27 (1 -439/1172/*- 839/9376*+ 12/293 * 2 )H(£2 

- 18752/27 (1 - 189/586/*- 1343/9376* + 12/293 * 2 )H 3 
+ 12176/27 (1 -239/761/* + 3553/6088* + 14/761 x 2 )^ 

+ 320/3 (l-l/5/x-3/10x)H_ 2 , 2 + 2528/9 (1 -61/474/x+21/158x 

- 10/237 x 2 )H- 2 fi + 832/3 ( 1 - 277/2808/* + 55/52 x + 5/702 x 2 )H_i, 

+ 320(1 -l/15/x-l/10x)H_ 3j o- 18896/27(1 - 15/1181/x- 151/4724 x 
+ 52/ 1 1 8 1 x 2 )H ,o,o + 832/3(1 + 1 /39/jc - 3/26 *)H_ 2 ,o,o - 608/3 ( 1 + 1 / 1 9/jc 
-5/38jc)H_ 2 ^ 2 + 1504/9 (1 + 19/282/* + 55/47 x- 10/141 x 2 )H- L0fi 

- 1216/9 (1 +23/152/* + 175/152*- 3/38* 2 )H_^ 2 - 192(1 + 1/3/x 

+ l/18x)H_ 2 ,-i,o + 5096/45 (l + 246/637/x-l/14x)C 2 -2048/9 (1 + 85/128/x 
+ 97/128*- l/32x 2 )H_i,_i,o- 160/9 (1 +3/x+ ll/20x)H 2 ^ 2 ) 
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+ C F nj ( Pgq (x) (8/3 £ 3 + 64/9 Hok- 16H OA o + 40/9H^ 2 - 64/9H 1)0 , 
- 40/9 Hi,i j0 - 32/9 Hi,!,! - 40/9 H12 — 64/9 H 2 , - 40/9 H 2 , 1 - 64/9 H 3 ) 
+ 13726/243 (1 -25291/13726/x- 5807/13726*) +5440/81 (1 -959/680/x 
-59/136jc)H + 1712/27 (1 -241/214/x- 101/214jc)H ,o 
+ 1036/27 (1 - 551/518/jc- 239/518x)H 1 - 848/27 (1 - 1/x- 25/53 x)(C 2 -Hi, 
-H 2 ) +536/27 (1 — 125/134/jc — 61/134jc)Hi j i) . (A.l) 

The corresponding expression for the gluon coefficient function is given by 

c^( x ) = cl (-44228/27(1 - 8/1 1057/x 2 + 8977/1 1057/x+ 11031/1 1057x 

+ 8959/1 1057x 2 )H_i, -7262/27(1 - 19412/3631/x- 842/3631jc+ 16705/3631jc 2 )^ 3 

+ 12724/27(1 - 48356/9543 fx- 16641/6362x + 67957/19086x 2 )H u 

-608/3(1 -275/57/x- 153/76* + 869/228x 2 )Hi 5 oC2 +7672/9(1 - 59035/ 17262/x 

- 739/41 1jc+ 12371/4932x 2 )Hi i0 + 544(1 - 77/34/jc - 1 13/34x + 66/17jc 2 )H 2 j,i 
+ 408(1 - 110/51/x-383/306x + 583/306x 2 )H li o,o,o + 5968/9(l -759/373/x 
-1955/1492jc + 2563/1492jc 2 )Hi, 3 -6200/9(1 - 1573/775/x-941/775jc 

+ 1397/775x 2 )H M ^ 2 - 4616/3(1 - 20489/10386/x - 1483/1 154* 

+ 1831 1/10386jc 2 )Hi£ 2 + 8128/9(1 - 913/508/jc - 623/508jc + 803/508x 2 )Hi j ; i i0 

+ 13448/9(1 - 8963/5043 fx - 4355/3362* + 7874/5043jc 2 )Hi i2 

+ 6544/9( 1 - 726/409/x - 4373/ 1636x + 5357/ 1636x 2 )H 2! 2 + 800( 1 - 44/25/x 

-61/50jc + 77/50jc 2 )H MiM + 8536/9(1 - 169/97 /jc- 2597/2134* 

+ 295/1 94x 2 ) (Hi, 1 ,2 + Hi )2 ,i ) + 1 3960/9 ( 1 - 2943 / 1 745/x - 87 1 /69Sx 

+ 5039/3490jc 2 )Hi,i,i +9832/9(1 -2013/1229/jc- 1426/1229jc+ 1826/1229jt 2 )Hi, 2 ,o 

+ 137744/27(1 - 13910/8609/x- 492137/413232* + 568297/413232x 2 )Hi 

+ 15400/9( 1 - 9326/5775/x - 4593/3850* + 7874/5775jc 2 )H u , 

+ 10904/9(1 -2189/1363/x- 1529/1363x + 2013/1363jc 2 )H MA0 

+ 1 6040/9 ( 1 - 304 1 /2005 jx - 462 1 /40 1 Ox + 5477 /40 1 Ox 2 ) H j >0)0 

+ 7304/9(1 - 123/83/x- 1916/913jc + 248/83jc 2 )H 2j i,o + 7544/9(1 -33/23/x 

- 1931/943x + 2541/943x 2 )H 2! o, + 7784/9(l-880/973/x-4183/1946x 

+ 1067/278x 2 )H 3il +2668/3(1 -583/667 /jc-3995/2001x + 7139/2001x 2 )H 3 ^ 

+ 68572/27(1 - 13324/17143/x- 599/553x + 28492/17143x 2 )H 2 ,o 

+ 62572/27(1 - 10990/15643 fx- 1421 3/ 15643x + 27766/ 15643x 2 )H 2 ,i 

+ 788(1 -11 1 1/1773/jc - 4202/1773* + 8558/1773jc 2 )H 4 + 9506/3( 1 - 23900/42777/jc 

- 13907/14259x + 25964/14259x 2 )H 3 -788(1 -935/1773/x-4160/1773x 
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+ 8558/1773jc 2 )H ,o^2 - 9506/3( 1 - 20684/42777/jc - 4463/61 1 Ix 
+ 25964/ 14259x 2 )H ^ 2 + 3 1 84/9(1 - 88/ 199/jc - 301 /39Sx + 253/398x 2 )H 1 _ 2)0 
+ 959401/81(1 — 306237/959401 fx - 230274/95940 \x + 32945 1 5/2878203x 2 )H 
+ 146410/27(1 - 16284/73205/x- 57243/73205* + 101818/73205jc 2 )H ,o,o 
-22004/9(1 -957/5501 fx- 1450/5501jc + 9372/5501x 2 )H ^ 3 

- 16168/9(1 -341/2021/x + 619/2021x+1661/2021x 2 )H_ 2 ,2 
+ 23116/9(1 - 891 /5779/x + 998/5779* + 4752/5779x 2 )H_ 2 ^2 
- 10360/3(1 -517/3885/jc- 13/2590jc + 7337/7770x 2 )H_ 2 ,o,o 

+ 1544(1 - 209/1737/x - 80/579*+ 1430/1737x 2 )H_ 2 ,-i,o + 63584/27(1 - 201/1987/x 
-655/1987x+1959/1987x 2 )H_ 2 , - 15722/3(1 -31615/1910223/x-477635/848988jc 

+ 3411235/7640892x 2 ) -7568/3(1 +4/129/x + 7/946x + 89/86x 2 )H_ 3 ,o 

+ 1344(1 + 143/378/jc + 263/252jc + 253/756x 2 )H_i,_i, 2 - 3464/9(1 + 198/433/jc 

+ 2302/433x-2970/433jc 2 )H , ,o,o- 14432/9(1 +39/82/x + 3871 /3608x 

+ 133/328x 2 )H_i_i^ 2 + 16532/9(1 +2321/4133/x + 4363/4133x + 2101/4133x 2 )H_i^ 3 
-2072(1 + 1441/2331/jc + 788/777jc+ 1408/2331x 2 )H_i, 3 + 22792/9(1 + 169/259/jc 
+ 5765/5698x + 331/518x 2 )H_ l!0 ^ 2 -2960/3(1 + 121/185/x + x+121/185x 2 )H_ l!2J 
-2384/3(1 +99/149/x + jc + 99/149jc 2 )H_i i2 ,o + 5552/3(1 +275/347 /jc+ 1497/1388* 
+ 1001/1388jr)H_i_i A o- 302995/81(1 + 244674/302995/* + 3 13876/302995* 
+ 19086/43285x 2 )H , - 6776/3(1 +72/77/x+157/154x+141/154x 2 )H_ 1A0 ,0 
+ 4364/9(1 + 1045/1091/x- 532/1091* -506/1091jc 2 )Hi£ 3 -4672/9(1 + 143/146/jc 
+ 715/584x + 451/584x 2 )H_i _i _i, + 5104/9(1 + 30/29/* + 368/3 19* 

+ 26/29* 2 )H_i,_ 2 , + 1848(1 +7127/6237/* + 2074/2079* + 7127/6237* 2 )H_i A0 
+ 1 72045/8 1 ( 1 + 243046/ 1 72045 fx + 207847/ 1 72045* - 26002/ 1 72045* 2 ) £ 2 
+ 5888/9(1 + 3557/2208/*+* + 3557/2208* 2 )H_i, 2 - 172045/81(1 +70154/34409/* 
+ 3402 1 9/ 1 72045* - 26002/ 1 72045* 2 ) H 2 - 2096/3 ( 1 + 9677/47 1 6/jc + 26 1 /262* 
+ 9677/47 16x 2 )H_ ^2- 17752/45(1 + 7172/2219/* + 6750/2219*- 6061/4438* 2 )£f 

- 800/9(1 + 2563/300/* + 47/50jc + 2563/300* 2 )H_i _i >0 + 404/9(1 + 3 1 13/ 101 fx 
+ 4613/101x-3927/101x 2 )H 2 C 2 -2752/3(1 - 13/43x)H_ 3 -1,0 

+ 6464/3(1 - 16/101x)H_ 3)0 , + 5440/3(1 -7/85*)H_ 4 ,o- 4384/3(1 - ll/137x)H_ 3 ^ 2 
+ 1856(1 - 3/58*)£ 5 + 3008/3(1 + l/47*)H_ 3i2 - 1760(1 + 131/165*)C 2 £ 3 
+ 6656/3(1 + ll/13x)H ,oC3 - 1472(1 + 85/69x)H 3jl)0 - 1344(1 + 9/7x)H 3 , 1;1 

- 4288/3(1 + 87/67jc)H 3:2 - 3616/3(1 + 157/1 13x)H 3A o - 5152/3(1 + 35/23*)H 4i i 

- 4384/3 ( 1 + 22 1 / 1 37*)H 4j0 + 29 1 2/3 ( 1 + 23/ 1 3*)H 3 C 2 + 1344( 1 + 2x)H ,o,o^2 
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- 1344(1 + 19/9jc)H 5 - 2336/15(1 +452/73x)H C| - 192(kH ,o,o,o,o 

+ pgg(-*) (504^5 + 1876/3^3 + 6464/8 1^ - 1888/3^3 + 5038/45^ + 1376/3H_ 4 , 

- 1 184H_ 3 ^2 - 832H_3 _i )0 + 3872/9H_ 3 , + 1 120H_ 3 , , + 768H_ 3 , 2 - 1712H_ 2 £3 

- 616H_ 2 ^2 - 1984/3H_ 2 ,-2,o + 1712H_ 2 ,-i^2 + 800H_ 2 ,-i,-i,o - 5104/9H_ 2 -i,o 

- 5408/3H-2 -1,0,0 - 1312H_ 2 _i )2 + 536H_ 2)0 - 1808H_2,o£2 + 704H_ 2 ,o,o 

+ 4144/3H_2,o,o,o + 2992/9H_ 2 ,2 + 1376/3H_ 2 ,2,o + 544H_ 2 ,2,i + 4576/3H_ 2 , 3 

-7436/9H_iC3-2144/3H_iC2 + 1528/15H_iC2-832H_i,_3,o + 800H_i_2,_i,o 

+ 1712H_i -2^2 - 5192/9H-1 -2,0 - 5408/3H_i,_ 2 ,o,o - 1312H_i,_ 2 , 2 + 6448/3H_i,_i£ 3 

+7568/9H_i_iC2 + 800H_i_i_2,o-2112H_i_i_iC2-768H_i_i_i_i >0 

+ 6688/9H_i_i _i,o + 2336H_i_i _i,o,o + 1728H_i_i _i,2 -2144/3H-1 _i,o 

+ 7648/3H_i,_i,o^2 - 8624/9H_i -1,0.0 - 2096H_i,_i, ,o,o - 1408/3H_i _i )2 

- 2176/3H_i -1,2,0 - 2432/3H_i,_i,2,i - 6496/3H_i ,_i, 3 + 12928/81H_i, 
-4784/3^1^3-1936/3^1,0^2 + 2680/3^1,0,0-4400/3^1,0,0^2 

+ 5500/9H_ 1 , ,o,o + 2864/3H_i A oAO+1072/3H_i,2 + 48H_ 1 ,2^2 + 704/3H_ 1 ,2,o 

+ 256H-i,2,o,o + 704/3H-i,2,i +704/3H_i, 2 ,i,o+ 192H_i, 2 ,i,i +704/3H_i, 2 , 2 

+ 528H_i,3 + 672H_i, 3 ,o + 2336/3H_i, 3 ,i + 1328H_i, 4 + 120H + 5360/27H ^ 2 

+ 32 12/9Ho^ 3 +904/1 5HoC 2 2 - 6464/81H , + 1712/3H ,oC3 + 1804/9H ,oC2 - 268H ,o,o 

+ 432H ,o,o^2 - 1364/9H ,o,o,o - 208H , ,o,o,o - 496/3H 2 ^ 3 - 836/9H 2 ^ - 1 12H 2 , C 2 

- 176/9H 2A o + 64/3H 2 , ,o,o - 304/3H 2 ,iC 2 + 88/9H 2 ,i,o + 64/3H 2 ,i, ,o - 88/9H 2 , 2 
+ 64/3H 2 , 3 - 536/3H 3 - 608/3H 3 £ 2 - 308/3H 3 , - 352/3H 3 , ,o - 352/3H 3 ,i 

- 352/3H 3 ,i, - 96H 3 ,i,i - 352/3H 3 , 2 - 1628/9H 4 - 784/3H 4 ,o - 976/3H 4 ,i - 400H 5 ) 
+ p gg (x)(- 1616486/729 - 480£ 5 + 21994/9£ 3 + 170666/81^2 - 2896/3^2^3 

- 25102/45^ + 1 184/3H_ 4 ,o - 960H_ 3 £ 2 - 640H_ 3 ,_i, + 3520/9H_ 3 , + 896H_ 3 , , 
+ 640H_ 3 ,2 - 2144/3H_2^ 3 - 572/3H_ 2 C 2 - 544/3H_ 2 _ 2 , + 752H_ 2 ,_i£ 2 

+ 608/3H_ 2 ,_i,_i,o - 616/3H_ 2 ,_i,o - 2048/3H_ 2 ,-i, 0| o - 1952/3H_ 2 ,_i, 2 
+ 13400/27H_ 2 ,o - 2320/3H_ 2 ,oC 2 + 2552/9H_ 2 , ,o + 1312/3H_ 2 , ,o,o + 88H_ 2 , 2 
+ 544/3H_ 2 , 2 , + 608/3H_ 2 , 2 ,i + 2048/3H_ 2 , 3 - 191858/81H + 1 1792/9H C3 
+ 6656/3H ^ 2 - 2408/5H Cf - 1305 14/8 1H , + 2224/3H ,oC3 + 1 1264/9H ,o^ 2 

- 10952/9H OA o + 688H , ,oC2 - 5060/9H ,o,o,o - 272H , ,o,o,o - 192268/8 lHj 
+ 9196/9H^ 3 + 17236/9Hi^ 2 - 296H 1 ^ 2 2 + 1664/3Hi _ 3>0 - 400H 1 ,_ 2 ^ 2 

- 1696/3Hi _ 2 ,_i,o + 3080/9Hi _2, + 1792/3Hj _ 2 , ,o + 352/3Hj _ 2 , 2 - 176374/81H 1 , 
+ 2032/3Hi,oC3 + 12760/9Hi,oC2 - 47888/27Hi, ,o + 688Hi, ,o£ 2 - 9196/9Hi i0 . ,o 

- 1456/3Hi, ,o,o,o - 19246/9Hi,i + 240Hi,i£ 3 + 10736/9^,^2 + 32H 1 , 1 ,_ 2 ,o 
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- 19484/9Hi i1j0 + 640H Mi oC2 - 10384/9H l5 i A o - 592H 1 , 1A o,o - 2040H Uil + 768Hl U ^ 2 
-14696/9H 1 , 1)l! o-1216H 1)1;1A o-4400/3Hi ;1) i ;1 -1152H 1)1)1)1)0 -960H 1;1) i ;1)1 

- 1 152Hi, u , 2 - 14080/9H M , 2 - 1280Hi,i, 2 , - 1152Hi,i, 2 ,i - 1024H M , 3 - 20452/9Hi, 2 
+ 880Hi, 2 C2 - 14432/9Hi, 2 ,o - 3296/3Hi, 2 , ,o - 1584Hi A i - 1280Hi, 2 ,i,o 

- 1152Hi,2,i,i - 1280Hi )2 , 2 - 13816/9Hi, 3 - 3488/3Hi, 3> o - 1216Hi )3 ,i - 2480/3Hi, 4 

- 177130/81H 2 + 1568/3H 2 ^ 3 + 3520/3^2 + 96H 2 - 2 ,o - 16424/9H 2 , 

+ 752H 2 ,o^2 - 8360/9H 2 , ,o - 1616/3H 2 , OA o - 6656/3H 2)1 + 880H 2)1 £ 2 - 14432/9H 2)1 , 
-3296/3H 2 ,i A o- 14344/9H 2 ,i,i - 1280H 2 ,i,i, - 1152H 2 ,i,i,i - 1280H 2 ,i, 2 - 13112/9H 2 , 2 

- 1 184H 2 , 2 ,o - 1280H 2 , 2 ,i - 3104/3H 2 , 3 - 60440/27H 3 + 2432/3H 3 ^ 2 - 3784/3H 3 , 

- 832H 3 , ,o - 14344/9H 3 ,i - 3680/3H 3 ,i 5 o - 1280H 3 ,i,i - 3680/3H 3 , 2 - 1 1440/9H 4 
-2704/3H 4 ,o-3568/3H 4 ,i -720H 5 ) + (1 -x)(-1504/3H_ 2 C 3 - 1088/3H_ 2 ,_ 2 , 

+ 384H_ 2 ,_i£ 2 + 896/3H_ 2 ,_i,_i,o - 2656/3H_ 2 ,_i A o -704/3H_ 2 ,_i, 2 - 2560/3H_ 2 , C 2 
+ 3584/3H_ 2 ,o,o,o + 704/3H_ 2 ,2,o + 256H_ 2 ,2,i + 1952/3H_ 2j3 ) + (1 +x) (896/3H 2 ^ 3 

- 128/3H 2 ,-2,o + 512H 2 ,o^2 - 1600/3H 2 , ,o,o + 2240/3H 2 ,iC2 - 3296/3H 2) i,o,o 

- 896H 2 ,i,i, - 768H 2 ,i,i,i - 896H 2 ,i,2 - 992H 2 , 2 ,o - 896H 2 , 2 ,i - 2144/3H 2 , 3 ) 

+ 8(1 14160613/2916 + 4048/9^5 -90632/27^ 3 +464/3^ 3 2 - 196718/81^2 
+ 836^2^3 + 22501/27^| - 1 5644/105^1 } ) 



+ C 2 n / (- 15128/27(1 -4/1891/x 2 + 375/1891/x+ 1893/1891x + 381/189U 2 )H_ li0 
-631/27(1 - 264040/5679/jc - 45427/ 1 893* + 1 32 1 60/5679jc 2 )H 1 
-428/9(1 -848/107 /x-481/107x + 534/107x 2 )Hi, 2 -416/9(l - 102/13/x -215/52* 
+ 127/26x 2 )H U j-92(1-2692/621/x-491/207jc+1786/621x 2 )H U! o 
+ 26362/8 1 ( 1 - 34234/ 1 3 1 8 1 fx - 62089/ 1 3 1 8 he + 3 1748/ 1 3 1 8 lx 2 ) & 
-3284/27(1 -2048/821/x- 15 13/821* + 2594/821x 2 )H 2;1 - 9536/27(1 - 5465/2384/x 

- 9027/4768jc + 857/596jc 2 )Hi, - 26362/81(1 - 29734/13 181/jc- 39373/13181jc 
+ 3 1748/ 1 3 1 8 lx 2 )H 2 + 1228/9(1 — 668/307 jx - 68 1 /307x + 360/307x 2 )Hi ^ 2 

- 460/3 ( 1 - 2 1 28/ 1035/jc - 3 1 / 1 5* + 944/345jc 2 )H 2 , - 452/3 ( 1 - 656/339/x - 57 1 /339jc 
+ 442/339jc 2 )H 1 ,o, - 10538/27(1 -30380/15807 /jc-9705/5269jc+18482/15807x 2 )Hi,i 

- 13184/27(1 -599/412/x-22057/6592x + 3555/1648x 2 )H , + 1600/9(1 -9/10/x+x 

- 87/100x 2 )H_ 1 ,_ 1 ,o + 64/3(1 -2/3/x- 1/8x + 37/24jc 2 )H 1A o,o - 248/3(1 - 16/31/jc 

- 17/31x + 30/31jc 2 )Hi, ^ 2 -218/45(1 -48/109/x- 1782/109*- 804/109x 2 )^f 

+ 3472/27(1 - 185/434/x-53/7jc + 85/31x 2 )^ 3 - 1 168/3(1 - 226/657/jc + 21 8/219jc 
_443/1314x 2 )H_i A0 + 776/9(1 - 28/97 /x-109/194x + 71/97x 2 )Hi, 3 
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- 1672/9(1 -3/ll/* + 206/209*)H ^ 3 +3968/9(1 -61/248/* + *- 119/496* 2 )H_i^ 2 
-2312/9(1 -64/289/x + 266/289x-42/289x 2 )H_ 1 _ 1 _ 1)0 - 120(1 - 28/135/*+ 14/15* 

- 19/135* 2 )H_ 1AO ,o + 976/3(1 - 12/61/* + 225/244*- 15/122* 2 )H_! _ 1>0)0 

+ 1664/9(1 - 5/26/* + 191/208* - 3/26* 2 )H_i _ 2 ,o - 2524/9( 1 - 1 12/63 l/x + 584/631* 

- 66/63 1* 2 )H_ + 2440/9( 1 - 43/305 fx + 569/6 10* - 23/305* 2 )H_ i ^ 3 
+ 152(1 -8/57/* + 53/57*-4/57* 2 )H_i_i j2 -700/3(l-2/15/* + 761/525* 
+ 4/35* 2 )H 0! o^2 + 904/3(1 - 34/339/* + 38/113*- 20/1 13* 2 )H_ 2 ,-i,o 

- 776/9(1 - 9/97 fx + 241 /91x + 7 1 /97* 2 )H 2 £ 2 + 1268/9(1 - 28/3 17/* + 304/3 17jc 

- 14/317* 2 )H_!^ 2 - 352(1 - 8/99/*+*- 8/99* 2 )H_ 1;2 + 700/3(1 - 38/525/* 

+ 977/525* + 4/35* 2 )H 4 + 1520/9(1 -6/95/* + 40/19*)H ,o,o,o-968/3(l - 16/363/jc 

- 36/121*- 70/363* 2 )H_ 3 , - 4808/27( 1 - 26/601/* + 740/601* - 630/601* 2 )H_ 2 , 
+ 2056/9(1 - 10/257/*+ 116/257* + 32/257* 2 )H 2! i, + 712/3(l -3/89/* + 677/534* 

+ 31/267* 2 )H 3! o + 548/9(1 -4/137/*-52/137jc + 90/137jc 2 )Hi j i, 2 + 2368/9(1 -1/37/jc 

+ 817/592* + 53/296* 2 )H 3J - 1028/3(1 +4/257/* + 35/257*- 158/771* 2 )H_ 2A0 

+ 740/9(1 +4/185/*- 23/37* + 64/185* 2 )Hi A0 + 644/9(1 +4/161/*- 76/161* 

+ 82/161* 2 )H li2;1 +2132/9(1 + 16/533/* + 368/533* + 2/13* 2 )H 2)2 

-412/3(1 +4/103/* + 99/103* + 8/103* 2 )H_i, 3 + 3164/9(1 +6/113/* + 4/113* 

- 152/791* 2 )H_ 2 ^ 2 -836/9(1 + 12/209/*- 271/209*- 4/1 l* 2 )H U! o, 

+ 952/9(1 +2/17/* + 571/238* + 6/17* 2 )H 2A o- 1808/9(1 + 19/1 13/* -43/226* 

-23/113* 2 )H_ 2i2 - 160/3(1 + 1/3/*+*+ 1/3* 2 )H_i j2 j -40(1 +4/9/*+* 

+ 4/9* 2 )H_i i2 ,o + 608/9(l + 17/38/*- 107/76*- 33/38* 2 )HijC 2 

+ 496/9(1 + 77/62/* - 56/3 1* - 63/3 l* 2 )Hi £ 3 + 208*H ,oC3 - 32*H , ,o,o,o 

+ 16306/81(1 +52501/8153/* + 78053/16306*-22866/8153* 2 )H 

+ 544/27(1 +309/34/*+ 1667/34*- 1903/68* 2 )H , , 

+ 28211 /243 ( 1 + 1 582888/84633/* + 47572 1 /56422* - 1 63987 1 / 169266* 2 ) 

- 88/27(1 + 804/1 1/* + 6137/22* - 1624/1 1* 2 )H ^ 2 + 88/27(1 + 830/1 1/* 
+ 4657/22* - 1624/ 1 1* 2 )H 3 - 1600/9(1 - 23/25* + 7/ 100* 2 )Hi _ 2)0 

+ 240(1 - 11/15*)H_ 3 _i j0 + 824/9(l -62/103* + 40/103* 2 )Hi,i i i i o 

+ 64(1 - 1/2*:+ l/2x 2 )Hi,i,i,i -704/3(1 - 17/44jc)H_ 3 ,o,0 + 728/3(1 - 19/91*)H_ 3 £ 2 

- 656/3(1 - l/4Lc)H_4 )0 - 368/3(1 + 7/23*)H_ 3;2 - 80(1 +31/75*)^^ 

+ 224(1 +4/7*+ 1/7* 2 )H 2 jj +376/3(1 +31/47*)H 3 ^ 2 - 232/3(1 +21/29*)H 3)0 , 
-48(1 + 11/9*)H 5 +48(1 +4/3*)H OA oC2 - 232/3(1 +65/29*)^ 3 -56/3(1 +27/7*)^ 
+ p gg (-*)(-280/3C 3 - 896/81^2 - 916/45£f - 704/9H_ 3 .o + 1 12H_ 2 ^ 2 + 928/9H_ 2! _i, 
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- 80H_ 2 ,o - 128H_ 2 ,o,o - 544/9H_ 2 ,2 + 1352/9H_i^ 3 + 320/3H_i^ 2 + 944/9H_i _ 2 ,o 
-1376/9H_i_iC2-1216/9H_i_i_i )0 + 320/3H_i_i )0 + 1568/9H_i _i )0) o 

+ 256/3H_i _i, 2 - 1792/81H_ lj0 + 352/3H_ lj( £ 2 - 400/3H_ 1A o - 1000/9H_ 1A o,o 

- 160/3H_i, 2 - 128/3H_i j2j o - 128/3H_i, 2 ,i - 96H i, 3 - 48H - 584/9H ^ 3 

- 800/27H ^ 2 + 896/8 lH ,o - 328/9H ,oC2 + 40H , ,o + 248/9H , ,o,o + 152/9H 2 ^ 2 
+ 32/9H 2 , ,o - 16/9H 2 ,i,o + 16/9H 2j2 + 80/3H 3 + 56/3H 3j0 + 64/3H 3 ,i + 296/9H 4 ) 
+ p gg (x) (1234307/1458 - 8608/27^ 3 - 55712/81^ + 4852/45^ - 640/9H_ 3 , 

+ 104/3H_2^ 2 + 112/3H_ 2 ,_i,o - 2000/27H_ 2 ,o - 464/9H_ 2 ,o,o - 16H_ 2 ,2 

+ 68830/8 1H - 992/9H ^ 3 - 4672/9H C2 + 44684/8 1H , - 2048/9H ,oC2 

+ 2776/9Ho,o,o + 920/9H 0) o,o,o + 67730/81Hi-520/9Hi^3-4544/9Hi^2-560/9Hi_2,o 

+ 56284/81Hi,o - 2320/9Hi, C2 + 12128/27Hi A o + 1672/9Hi A o,o + 6212/9H M 

- 1952/9Hi,iC2 + 4672/9Hi,i, + 1888/9Hi,i, 0) o + 512Hi,i,i + 2672/9Hi,i,i,o 

+ 800/3H UiU + 2560/9Hi,i, 2 + 5024/9Hi, 2 + 2624/9Hi j2 , + 288Hi j2 ,i + 25 12/9H 1;3 
+ 56608/81H 2 - 640/3H 2 ^ 2 + 1360/3H 2 , + 1520/9H 2)0 ,o + 1616/3H 2J + 2624/9H 2 ,i i0 
+ 2608/9H 2 ,i,i + 2384/9H 2;2 + 14096/27H 3 + 688/3H 3j0 + 2608/9H 3 ,i + 2080/9H 4 ) 
+ (1 -x)(488/3H_2^ 3 + 368/3H_ 2! _ 2 ,o - 544/3H_ 2 ,_i£ 2 - 544/3H_ 2 ,-i,-i,o 
+ 216H_ 2 ,_i,o,o + 272/3H_ 2 ,_i, 2 + 224/3H_ 2 ,oC2-80H- 2 ,o, 0) o+16/3H- 2 , 2 ,o 

- 152/3H_ 2 , 3 + 16/3H 3 ,i,o - 16/3H 3)2 ) + (1 +x) (464/3H 2 C 3 - 224/3H 2 ,_ 2 , 

+ 176/3H 2 ,oC 2 - 32H 2 , ,o,o + 304/3H 2 ^ 2 - 232/3H 2 ,i, , + 32/3H 2 ,i, 1 , - 32/3H 2 , li2 
+ 40/3H 2 , 2 ,o-104/3H 2 , 3 -24H4,o-40H4,i) + 8(l -jc){ - 1988899/972 
+ 1064/9^5 + 6794/27^3 + 70774/81^2-24^3 -28526/135^|}) 



+ C F C A n f ( 11528/9(1 -4/1441/x 2 + 122/1441/* + 1384/1441* + 69/1441x 2 )H_l 

- 56/9(1 - 309/7/jc - 78/7* + 50jc 2 )H m , - 340/27(1 - 1338/85/x - 643/17jc 
+ 3304/85jc 2 )H 2;0 - 68/3(1 - 1450/153/jc - 41/17jc + 98/9jc 2 )Hi jU 

+ 2362/27 (1 - 2384/1 181/x + 3209/ 1 1 8 lx + 76 1 6/ 1 1 8 lx 2 )H ^ 2 

- 2362/27 ( 1 - 2256/ 1 1 8 1 fx - 1 843 / 1 1 8 Ix + 76 1 6/ 1 1 8 Ijc 2 ) H 3 + 800/3 ( 1 - 67/75/x + x 

- 67/75x 2 )H_i,_iC 2 - 169601/243(1 - 121601/169601/x- 89143/16960 Ix 

+ 2769 1 2/ 1 6960 lx 2 ) H -3952/9(1 - 141/247/x + 265/247* - 123/247x 2 )H_i _ lj0 
-904/3(1 -58/113/x+x-58/113jc 2 )H_^ 3 +6128/9(1 -466/1149/x + 847/766x 
- 689/2298jc 2 )H_i A o - 1400/3(1 - 212/525/x+x- 212/525x 2 )H_i _i A0 
-6712/9(1 - 259/839/jc + 91 1 /839x - 187/839x 2 )H_^ 2 + 544(1 -46/153/jc + jc 

- 46/153x 2 )H_i _! _ lj0 - 1 136/3(1 - 20/11/x + x- 20/71x 2 )H_i ,_ 2 , 
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- 23 150/27 ( 1 - 2664/ 1 1 515 /x - 5048/1 1 575jc + 9784/ 1 1575x 2 )H , ,o 

+ 4736/9(1 -59/296/* + 323/296* -4/37x 2 )H_ li2 - 1712/3(1 -20/107 /x + 66/\01x 

- 68/321x 2 )H_2,-i,o - 664( 1 - 44/249/x - 20/83* - 100/747jc 2 )H_ 2 ^2 

+ 1136/3(1 -12/71/x-63/71x-16/213x 2 )H_2,2 + 656(1 -44/369/*+ 1/246* 

- 26/369x 2 )H_ 2 ,o,o + 1904/3(1 - 8/1 \9/x - 39/1 19* - 12/1 19jc 2 )H_ 3 , 

+ 4384/9( 1 - 8/4 1 1 jx + 42 1 /548x - 475/822x 2 )H_ 2 ,o + 1 1 84/3 (l + 4/1 1 l/x - x 

-4/11 lx 2 )Hi _2, + 13612/45(1 + 352/3403 jx - 2282/3403* - 776/3403* 2 )^ 2 2 

+ 22804/27(1 +4697/34206/jc- 9101/1 1402jc-4828/17103jc 2 )H u 

+ 20020/27(1 + 159/910/*- 7283/10010* - 156/385* 2 )Hi )0 

+ 384329/243(1 + 455021/2305974/*- 19253/69878jc- 486427/1 152987* 2 ) 

+ 3140/9(1 + 192/785/* + 1919/785*)H ,o,o,o- 78058/81(1 + 10992/39029/* 

- 68023/39029jc - 9282/39029* 2 ) £ 2 + 2960/9( 1 + 1 2/37/* +104/1 85* - 54/ 1 85* 2 )H 3 ,i 
_ 4496/9( 1 + 92/28 l/x - 83/562jc - 76/28 1* 2 )H 2 ^ 2 + 3976/9( 1 + 172/497 /* 

- 13/497*- 152/497x 2 )H 2 , ,o + 2752/9(1 + 15/43 /* + 197/344*- 12/43* 2 )H 3;0 

+ 3100/9(1 +288/775/jc+ 1339/775*- 168/775jc 2 )H 4 + 78058/81(1 + 15384/39029/jc 

- 18199/39029*-9282/39029* 2 )H 2 - 1112/3(1 +68/139/*-*- 68/139* 2 )H M C 2 
-3100/9(1 +384/775/jc+ 1807/775*- 168/775* 2 )H ,oC2 + 1928/9(1 + 124/241/jc 
+ 115/241* - 84/241* 2 )H 2 , 2 + 344/3(1 + 24/43/* + x + 24/43* 2 )H_i A o,o 

+ 2024/9(1 + 148/253/* + 97/253*- 100/253* 2 )H 2i i 5 o + 664/3(1 + 148/249/*-* 

- 148/249* 2 )Hi, 3 + 1648/9(1 +64/103/* + 34/103*-40/103* 2 )H 2;U 

+ 280( 1 + 40/63/* - x - 40/63x 2 )Hi i 1>0|0 - 640/3(1 + 5/6/x -x- 5/6* 2 )Hi j0 C2 
-760/3(1 +50/57/*-*-50/57* 2 )H^ 3 + 152/3(1 +52/57/*+* + 52/57* 2 )H_ u 
+ 4468/27(1 + 1062/1117/* + 2947/1117*-2918/1117* 2 )H 2i1 + 296/3(1 + 112/111/* 
-x - 1 12/1 1 l* 2 )H Mj2 - 2888/9(1 + 1 1 14/1083/* - 306/361* - 1369/1083* 2 )Hi£ 2 
+ 224/3(1 + 22/21/* + * + 22/21* 2 )H_ li2jl + 31550/81(1 + 16644/1 5775 /* 

- 945 1 /3 1 55* - 2784/ 1 5775* 2 ) H , +136(1+ 200/ 153/* - x - 200/ 1 53* 2 )H i ,0,0,0 

+ 8(1 + 4/3/*-*-4/3* 2 )(22/3H M!M +7Hi )2) i) + 160/3(1 + 4/3/* + * + 4/3* 2 )H_i, 2 , 

- 164(1 +592/369/* + 455/123*-64/123* 2 )H C3 + 1792/9(1 + 11/6/* 

- 143/224jc- 1475/672x 2 )Hi A o + 152/3(1 + 112/57/jc-x-112/57x 2 )H 1)1i i i o 

+ 48(1 +2/jc-x- 2jc 2 )Hi j2i0 + 304/3(1 +691/342/jc- 41/1 14jc-500/171jc 2 )H!, 2 
-4262/27(1 +6900/213 l/x- 16385/2131*- 1952/21 3 1* 2 )^ 3 
+ 15112/81(1+ 33767/7556/* - 4061/7556* - 29945/7556* 2 )H! + 16/3(1 + 88/3/* 
+* + 88/3x 2 )H_i _i, 2 -256/3(1 - 155/16*)H 5 + 256/3(1 - 149/16*)H , ,oC2 
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_480(1 -47/30x)Ho,o,o,o,o- 544(1 - 14/17x)H_ 3 _i )0 + 304(1 -9/19x)H_ 3 ,o,o 

- 400(1 -8/25x)H_ 3 ^2 + 1024/3(1 +3/32x)H_ 4; o -448(1 +31/168x)^ 5 

+ 128(1 +3/4jc)H_ 3j2 + 560/3(1 +437/350x)H C| - 1280/3(1 + 107/80x)H 3 C 2 
+ 384(1 +25/l&c)H 3i0 ,o + 1232/3(1 + 138/77^)^3 + 560/3(1 +71/35jc)H 3i i, 
+ 464/3(1 +65/29jc)H 3;2 + 416/3(1 + 31/13x)H 3) i,i -496/3(1 + 163/31jc)H ,oC3 
+ 112(1+ 39/7x)H 4 ,i + 208/3(1 + 109/13x)H 4 ,o + P sg (-x) (-224/5^ - 128/3H_ 3 , 

- 24H - 128/3Ho^ 3 - 64/3Ho,oC 2 + 64/3Ho,o,o,o) + P gg (x) (7810/27 -2152/9^ 

- 100/3^2 + 32/15^1 - 256/3H_2^2 + 128/3H_ 2 ,o,o + 256/3H_ 2 ,2 + 2020/9H 

- 64Ho^ 3 - 176/9H C 2 + 24H , + 64H ,o^2 - 64/3H OA o,o + 598/3HJ - 128H^ 3 
+ 32Hi ;0 + 32Hi,i + 100/3H 2 - 176/9H 2)0 + 176/9H 3 - 128/3H 4 ) 

+ (1 -x) (-1 120/3H_ 2 ^3 - 992/3H_ 2 ,-2,o + 1408/3H_ 2 ,-i^2 + 1472/3H_ 2 ,-i,-i,o 

- 496H_ 2) -i )0 ,o -224H_ 2 ,-i,2 - 128/3H_ 2 ,o£ 2 - 16/3H_ 2 , , ,o - 64H_ 2 ,2,o 
-64H_2,2,i-16H_2 J 3) + (l+x)(-128/3H_ lj oC2-416H 2 C3 + 640/3H 2) _2 ) o 

- 992/3H 2 , ^2 + 880/3H 2 , ,o,o - 1232/3H 2 ,iC2 + 1072/3H 2) i,o,o +400/3H 2 ,i,i,o 
+ 352/3H 2 ,i,i,i + 496/3H 2 ,i )2 + 128H 2 , 2 ,o + H2H 2)2) i + 272H 2 , 3 ) 

+ 8(1 - x) { - 41 1 137/324 + 200^ 5 + 6664/9^ 3 + 2128/9^2 - 128^2^3 + 176/45^}) 



+ Cjn f (^-12016/9(1 - 2/2253/* 2 + 5/ 1502/x + 4505/4506* + 3/75 U 2 )H_ 1>0 

- 16/3(1 - 13/x-5x+ 17x 2 )H 1A0 + 44/3(1 -32/1 l/x+70/llx- 128/33x 2 )H 0) o,o,o 
-256/3(1 - 19/12/x+x- 19/12x 2 )H_i AO ,o + 592/3(1 - 172/111/x+x 

- 172/1 1 lx 2 )^! _i )0)0 - 296/3(1 - 56/37/x - 94/37x - 1 12/1 1 1* 2 )H C3 

+ 452/3(1 -280/339/x- 137/113x + 352/339jc 2 )H^ 2 + 704/3(1 - 20/33/jc+jc 

- 20/33jc 2 )H_i _ 2 ,o - 320(1 - 8/15/x+x- 8/15x 2 )H_i,_i,_i, - 192(1 -4/9/x- 13/6x 

-2/3x 2 )H_ 2 ,2 + 352(1 -4/11/x- 10/11x-20/33x 2 )H_ 2 ^2 -592/3(1 - 12/37/x 

+ 19/74jc-14/37jc 2 )H ,o^2 + 320(1 -4/15/x + 3/5jc-8/15x 2 )H_2-i,o 

+ 1424/3(1 - 70/267 /x + 79/89x- 100/267x 2 )H_i _i ;0 - 1216/3(1 -4/19/x-21/76x 

-37/57jc 2 )H_ 2 ,o,o + 712(1 - 122/801/jc + 257/267jc- 152/801x 2 )H_^ 2 

+ 592/3(1 -4/37/x + 51/74x- 14/37x 2 )H 4 - 1424/3(1 -26/267 /x + x-26/267jc 2 )H_i, 2 

- 1216/3(1 -2/19/x-4/19jc-26/57x 2 )H_ 3 , -9104/9(1 -48/569/* + 554/569* 
-63/569x 2 )H_i,o,o- 136(1 -4/51/x-x + 4/51x 2 )Hi !3 +2141/9(l+264/2141/x 
+ 42253/6423x-9776/6423x 2 )H ,o + 680/3(l + 16/85/x + 77/85x-28/85x 2 )H 3a 

+ 48631/162(1 +951 3/4863 1/jc- 1 270/4863 lx- 58494/4863 \x 2 ) +200(1 + 16/75/jc 
+ 24/25* - 8/25* 2 )H 3i0 + 7442/27 ( 1 + 836/372 1 jx + 773/7442* - 1 8254/ 1 1 1 63x 2 )H 
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+ 136(1 + 4/17/jc-x-4/17x 2 )H m ^ 2 +488/3(1 + 16/61/x + 49/61x-52/183x 2 )H 2 ,i, 
+ 200(1 +4/15/x + 58/75x-56/225x 2 )H 2 ,ij +216(1 +8/27 /x + 7/9x-20/81x 2 )H 2 ,2 
-56(1 + 8/21/*+ 33/7jc + 4/7x 2 )H 2 ^ 2 + 376/3(1 + 68/141/x -x - 68/ 141jc 2 )Hi i0 ^ 2 
+ 64/3(1 + 1/2/x + 39/8x- 1/2jc 2 )H 2A o + 7141/27(1 + 11065/21423/x+1340/7141jc 

- 36508/21423jc 2 )Hi + 220(1 + 284/495/x + 17/33jc - 472/495x 2 )H 2 ,o 

+ 78(1 + 272/351/x + 911/117x- 1312/351x 2 )H 3 -78(1 + 272/35 l/x + 445/39x 

- 1312/351x 2 )H C2 + 416/3(1 + 71/78/x + 29/52x-22/13x 2 )H 2 j + 152/3(1 +52/57/x 
-x-52/57x 2 )H ljl) i i0 - 112(1 +20/21/x + jc + 20/21x 2 )(2H_i_i !2 -H_i,3) 

+ 64(l + l/x-x-x 2 )H l!2 ,o + 856/9(l+374/321/x+1379/428x- 1321/321x 2 )H 2 
-856/9(1 + 389/32 1/x + 221 57/ 1284jc- 1321/321jc 2 )^ 2 + 224/3(1 + 10/7x)£ 5 
+ 8(1 + 4/3/x-x-4/3x 2 )(H MA0 + 14/3Hi,i,i,i +7Hi )2 ,i) -368/3(1 + 100/69/x+x 
+ 100/69x 2 )H_i ; oC2 + 820/9(1 + 3 16/205 fx - 49/205x - 472/205x 2 )H M , - 576H_ 4 ,o 

- 128/3(1 + U/6/x-x- ll/6x 2 )H 1AO ,o + 812/9(1 +384/203 /x -59/203* 
-528/203x 2 )H U j +260/3(1 +28/13/x-21/65x- 184/65x 2 )Hi, 2 + 24(1 +20/9/x 
-x - 20/9jc 2 )Hi ,i )2 + 686/9(1 + 3208/1029/x + 241/343* - 4960/1029x 2 )Hi )0 

+ 64(1 + 14/3/x + x+14/3x 2 )H_ w C 2 + 188/15(1 + 808/141/x + 577/47x + 44/47x 2 )£f 
+ 40/3(1 + 28/3/x -x- 28/3jc 2 )Hi ^ 3 + 50/9(1 + 88/5/x - 4609/25* + 1 144/25x 2 )^ 3 
+ 70/3(1 +424/45/x + 667/105x-5284/315x 2 )H M - 8(1 +92/3/x + x + 92/3x 2 )H_^ 3 
+ 32/3(1 - llx) ^3 -304/3(1 - 1 4/ 1 9x) HoCf + 1600/3(1 - 16/25x)H_ 3 _ 1)0 
-608(1 - 13/57x)H_ 3 ,o,o + 1568/3(1 - 10/49x)H_ 3 ^2 -712/3(1 - 15/89x)H ,o^3 
+ 320/3(1 + l/10x)H 3 C 2 - 160/3(1 + l/5x)H 3)0 ,o -256(1 + l/4x)H_ 3 ,2 
-7120/9(1 + 159/445x - 20/89x 2 )H_ 2 ,o + 946/9(1 + 4570/473x - 912/473x 2 )H 0A0 
+ pgg{-x) (224/5^ 2 + 128/3H_ 3 ,o + 128/3Ho£ 3 + 64/3H ,oC2 - 64/3H ,o,o,o) 
+ Pgg (x) (-1 - 128/15CI + 256/3H_2^ 2 - 128/3H_ 2 , , - 256/3H_ 2 , 2 - 64H ,oC 2 

- 64H<£ 3 + 64/3H OAO ,o + 128/3H 4 ) + (1 -x) (336H_ 2 £ 3 + 896/3H_ 2i _ 2 , - 384H_ 2; -iC2 

- 384H_ 2 ,_i,_i,o + 1504/3H_ 2 ,_ 1 ,o,o + 192H_ 2) -i,2+544/3H_ 2 ,oC2 - 640/3H_2,o,o,o 

- 96H_ 2 , 3 - 256Hi _2,o) + (1 +x)(-616/3H ,o,oC2 + 120H OA o A o + 560/3H 2 £ 3 + 616/3H 5 

- 128H 2 ,-2,o + 144H 2 ,o^2 - 128H 2 ,o A o + 368/3H 2) i£ 2 - 16/3H 2 ,i A o + 80H 2 , M ,o 
+ 224/3H 2 ,i,i,i + 208/3H 2 ,i, 2 + 320/3H 2 , 2 ,o + H2H 2;2) i - 176/3H 2 , 3 + 128H 3 ,i,o 

+ 448/3H 3 ,i,i + 160H 3j2 + 176H 4j0 + 640/3H 4 ,i) + 8(1 -x) {751/9 - 400^ 5 + 740/3^3 }) 



+ C A nj (4498/8 1 (1 + 11 28/2249/x +131 /2249x + 524/2249x 2 ) (£ 2 - H 2 ) 

- 1286/27(1 + 1024/ 1929/x-99/643x + 524/ 1929x 2 )H u - 1264/27(1 +269/474/x 
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- 47/316* + 59/237* 2 )Hi ;0 - 3524/81(1 + 642/881/* + 428/881* + 56/881* 2 )Ho,o 
-217/3(1 + 18169/17577/*+ 104/837* -262/17577* 2 )H! -2170/27(1 + 1208/1085/* 
+ 214/651*+ 1/9765* 2 )H - 1112/27(1 - 68/139* + 69/1 39* 2 )^ 3 

- 2275 1 /243 ( 1 + 343 12/2275 1 /* + 17705/45502* - 108 17/45502* 2 ) 

+ 16/3(1 - 1/2*+ 1/2* 2 )(16/3H^ 2 - 10/3Hi A o -4H M , - 14/3Hi,i,i - 16/3Hi, 2 ) 

-704/27(1 - 19/44* + 21 /44* 2 )H 2 ,i -544/27(1 - 1 1/34*+ 15/34* 2 )H 2 ,o 

+ 568/27(1 - 19/71* + 30/71* 2 ) (H ^2 -H 3 ) + p gg (*)(- 138493/1458 - 296/27^3 

+ 584/9^2 - 6740/81H + 256/9HoC2 - 56H ,o - 176/9H ,o,o - 6652/81Hi + 304/9HiC 2 

- 584/9Hi i0 - 256/9Hi )0)0 - 584/9H M - 272/9H M) o - 32H M) i - 304/9H ii2 

- 584/9H 2 - 256/9H 2 ,o - 32H 2 ,i - 256/9H 3 ) - 200/27H , ,o(l +*) 
+ 8(1 - x) { 174769/972 + 404/9^3 - 6740/8 1 ^ 2 + 608/45^1 } ) 

+ C F nj ( - 1 84/45 (1 +x) X£ + 2684/27 ( 1 + 96/67 1 /* + 5 1 5/67 1* - 1 6/67 1* 2 ) (H ^ 2 - H 3 ) 
+ 3092/27(1 + 136/773/* + 71/773* + 244/773* 2 )^ 3 -3508/27(1 + 168/877/* 
+ 826/877* - 56/877* 2 )H , ,o - 1406/3(1 + 38524/170829/* - 3974/6327* 
-53617/170829* 2 )- 1240/27(1 +36/155/*+ 134/155*- 16/155* 2 )H 2 ,i 
-68290/243(1 + 1941/6829/* + 4931/34145*- 11746/34145* 2 )H 

- 14240/81(1 +57/178/*+ 1177/1780*- 21 1/890* 2 )H , - 1016/27(1 +48/127/* 
+ 166/127* - 32/127* 2 )H 2)0 + 56/9(1 + 8/21/* - 9/7* - 20/21* 2 ) (H^ 2 - Hi )2 ) 

+ 8020/8 1 ( 1 + 948/2005/* + 230/401* - 160/401* 2 ) (£ 2 - H 2 ) 

-3316/27(1 +9397/14922/*- 887/1658*- 6143/7461* 2 )H! -32(1 + 157/162/* 

-7/9*- 21 l/162* 2 )Hi j0 - 8/9(1+ 4/3/* -*-4/3* 2 )(7H 1 ,o,o + 4H li i il ) 

-24(1 + 346/243/* - 7/9* - 400/243* 2 )H u -8/3(1+ 32/9/* - 1/3* - 20/9* 2 )H U;0 

+ pg g (*) (-350/9 + 1 12/3^3 - 100/9H + 32/9H ^ 2 - 20/3HJ + 32/9H 2 ,o - 32/9H 3 ) 

+ (1 +*) (40/3H ^ 3 + 328/9H ,oC2 - 536/9H , ,o,o + 80/9H 2 £ 2 - 1 12/9H 2 , ,o 

- 80/9H 2)li0 - 64/9H 2 ,i,i - 80/9H 2 , 2 - 208/9H 3;0 - 176/9H 3)] - 328/9H 4 ) 
+ 8(1 -*)|28945/162 - 1 144/9^3 - 100/9^ 2 - 32/45^|}) 

+ nj (32/9(1 - 1 /2* + 1 /2* 2 ) (25/1 8 - L, 2 + 5/3H + H , + 5/3Hi + H 1)0 + H u + H 2 ) 
+ 16/9p gg (*) (25/18 - £ 2 + 5/3H + U Qfi + 5/3Hj + Hi, + H M + H 2 ) 
+ 8(1 -*){- 1000/729 + 80/27^}) • (A.2) 

It may be worthwhile to note that, unlike the third-order coefficient functions for gauge-boson 
exchange DIS [9, 10], Eqs. (IA.1I) and (IA.2I) do not involve additional special functions including 
terms of the form (l±*p"H v -j(*) with n > 1. 
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Appendix B: The third-order physical kernels for (F 2 , F§) 

In this second appendix we write down the exact expressions, again using the notations introduced 
in Section 3, for the matrix elements of the NNLO contribution to the physical evolution kernel for 
the system (F2,F§) of singlet structure functions. The coefficients in the expansion in powers of 
of the strong coupling are normalized according to Eq. (I4.ll ). 

As the corresponding NLO quantity before, we present K 21 via the decomposition (14.21) into 
its non- singlet and pure- singlet components. The former quantity, which is of direct interest for 
analyses of data on, e.g., the difference of the neutron and proton structure functions, is given by 

<w = 

ClC F (p qq (x) (50689/162 - 20/3 ^ - 1504/9 C2 + 192/5 C 2 2 + 16H_ 3 ,o + 64H_ 2 ,_ l!0 
+ 64H_ 2 ,o - 16H_ 2i0 ,o + 32H_ 2 ,2 + 10585/27 H + 64H C 3 - 212/3 H ^ 2 
+ 1988/9H ,o + 256/3H ,o,o- 16H ,o,o,o + 4649/27 Hi + 144H^ 3 -88H^ 2 
+ 96H 1; _ 2! o + 484/9H 1;0 - I6HK0C2 + 44/3 H 1A o - 48H 1A o,o + 484/9H ia 

- 88/3 Hi,i (0 - 64Hi,i,o,o + 88/3 Hi, 2 + 64Hi ;3 + 968/9 H 2 - 32H 2 ,o,o + 44H 3 + I6H4) 
+ p qq (-jc) (176/3 ^ 3 + 536/9 £ 2 - 16^| - 16H_ 3 ,o + 128H_ 2 ^ 2 + 352/3 H_ 2 , 

- 64H_2, 0) o - 128H_ 2 , 2 + 192H_iC 3 - 176/3 H_^ 2 - 256H_i - 352/3 H_i _i, 
+ 128H_i _i, ,o + 256H_ l5 _ li2 + 1072/9 H_ 1)0 + 176H_ li0 C2 + 176H_i )0)0 

- 48H_! A o,o - 128H_i, 3 - 124/3 H - 64H C 3 + 80/3 H0C2 - 536/9 H ,o - 48H ,o^2 

- 256/3 H ,o,o + 16H ,o,o,o + 32H 2 ^ 2 + 32H 4 ) + (1 +x)(-96H_ lj2 + 1 12H 2 - 8H 2 £ 2 

- 4H 2i0 ,o + 48H 3 + I2H4) + (1 -x)(8H_ 3 ,o - 8H_ 2 C 2 - 16H_ 2 ,-i,o + 8H_ 2 , , ) 
+ 16/3(1- 43x)H_ i,-i,o + 8/3 (1 + 43x) (H^ - H 1)0)0 ) - 4^2 (3 + 5x)H , 

+ 2^|(5 + 3jc) - 8^2(6 + 1 1jc)H + 8/3 (13 - 15x)H_ 2 , - 8/3 (19 - 25x)H_i )0 ,o 

- 4/3 ^ 3 (27 - 91jc) + 8/3 £ 2 (37 - 7jc)H_i - 8/15 £ 2 (210 - 145*+ 198jc 3 ) 

+ 11/9 (229 - 467x)Hi + 8/15 (575 + 22j*T 2 + 355* - 198x 3 )H_i, + 1/45 (2619 

- 52&T 1 - 3324 Ix + 4752x 2 ) H + 1 /270 (64993 + 3 1 6Sx 1 - 283643* + 285 1 2x 2 ) 

- 8/15 (145jc - 198x 3 )H ,o +jc(-16H<£ 3 + 40H ,o,o + 8H ,o,o,o) 

+ (20/27 (827 + 54Cs) - 1544/3 ^3 + 22286/27 ^ 2 - 1592/15 Q) 8(1 - jc)) 

+ C A Cj [p qq (x) ( 1 1 /2 - 664/3 ^3 - 276/5 ^ 2 2 - 1 6H_ 3 ,o - 48H_ 2 ^ 2 - 224H_ 2 ,_ 1 , 

- 128H_ 2j0 + 80H_ 2)0) o - 64H_ 2)2 - 3 17/2 H - 136H ^ 3 - 188/3 H ^ 2 - 74H , 

- 64H ,o^2 - 4/3 H ,o,o + 80H , ,o,o - 384H^ 3 + 176/3 H^ 2 - 256Hi,_ 2)0 

+ 1468/9H lj0 - 32H lj0 £ 2 + 584/3 H lj0 ,o + 176H lj0 ,o,o + 352/3H Mj0 + 128H MA0 
+ 176/3 H 1)2 - 128Hi, 3 + 1072/9 H 2 - 32H 2 ^ 2 + 440/3 H 2 , + 80H 2A0 + 88H 2 ,i 
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+ 128H 3 + 16H 3 , + I6H4) + Pqq(-x) (-244/3 ^ - 1072/9 £ 2 + 4^ 2 + 80H_ 3 ,o 

- 5 1 2H_ 2 ^2 - 64H 2 , 1 ,0 - 776/3 H_ 2 ,o + 336H_ 2 ,o,o + 480H_ 2 ,2 - 672H { ^ 
+ 208/3 H_i^ 2 - 64H_i _ 2 ,o + 896H_i_iC 2 + 704/3 H_i -1,0 - 576H_i _i A o 

- 896H_i _i )2 - 2144/9 H_i )0 - 672H_ li0 ^2 - 376H_i )0)0 + 272H_i A o,o + 48H_i, 2 
+ 32H_ lj2 ,o + 5 12H_i,3 + 284/3 H + 232Ho£ 3 - 124/3 H ^ 2 + 1072/9 H , 

+ 208H ,oC2 + 620/3 H , ,o - 80H , ,o,o - 1 12H2C2 - 24H 3 - 16H 3 , - I6OH4) 
+ (1 +x)(384H_ 1;2 + I6H2C2 - 44/3 H 2 , + 32H 2 , , - 48H 4 ) + (1 -x) (16H_ 2 ^2 
+ 32H_ 2 ,-i,o - 48H_2,o,o + 224/3 H 1)0 ) + 16/3(1+ 43x)H 1A0 - 32(2 - x)H_ 3;0 

- 16^3(3 - 2jc)H + 16^(3 + 5x)H , - 64/3 (5 - 17jc)H_i _i |0 - 32/3 ^ 2 (5 + 17jc)Hi 

- 8C|(7 + 6x) - 44/3 (7 + 25x)H ,o,o - 8/3 (1 1 - 75x)H_ 2 ,o - 16(22 + 35x)H 2 
+ 32/3 (29 + 7x)H_i A o - 32/3 £ 2 (41 + 19x)H_i - 44/3 (49 - 72x)Hi 

- 8/3 (65 + 101x)H 3 + 8/3 ^ 2 (65 + 176x)H + 4/3 £ 3 (99 + 61x) - 32/15 (280 

+ llx~ 2 + 170x-99x 3 )H_ 1; o + 16/15^ 2 (330 + 185x+ 198x 3 ) -2/45 (635 + 7115x 
+ 4752x 3 )H ,o + 1/90 (799 + 21 Ux' 1 + 89679x - 19008x 2 )H 

- 1/90 (10817 + 21 12JT 1 - 29627x+ 19008x 2 ) - 64xH , ,o,o + (2/3 (47 + I8OC5) 
+ 2296/3 ^ 3 - 1235/3 £ 2 + 16^3 - 856/15 8(1 -xj) 

+ Cj (/> qq (x) (72/5 Q - 32H_ 3i0 + 96H_ 2 £ 2 + 192H 2 , 1,0 - 96H_ 2A0 - 3H 
+ I6H0C3 - 24H C 2 + 26H , - 32H ,o,o,o + 192HiC 3 + 128Hi,_ 2 , - 96H 1)0 , 

- 64Hi, ,o,o + 64Hl2,o + 128H 1;3 - 48H 2 , + 32H 2 , ,o + 64H 2 ,i,o + 64H 2)2 + 64H 3)0 
+ 64H 3 ,i + 32H 4 ) + Pq q (-x) (-72^ + 56^ 2 - 96H_ 3 ,o + 512H_ 2 ^2 + 128H_ 2 ,-i,o 
+ 48H_ 2 ,o - 416H_ 2l o,o - 448H_ 2j2 + 576H_iC 3 + 96H_i£ 2 + 128H_i,_ 2i0 
-768H_i_iC2 + 640H_i_i )0) o + 768H_i_i ) 2 + 640H_i ) oC2 + 48H_i i o,o 

- 352H_ 1,0,0,0 - 96H_i )2 - 64H_i )2)0 - 512H_i )3 - 24H - 208H ^ 3 - 24H C 2 

- 224H ,oC 2 - 72H , ,o + 96H , ,o,o + 96H 2 ^ 2 + 48H 3 + 32H 3 , + 192H 4 ) 

+ ( 1 + x) (296/5 C 2 2 + 480H 1 £ 2 + 1 92H 1 _ lj0 - 32H 1 , - 4 1 6H_ 1 ,0,0 - 3 84H 1;2 

- 48H 2 , ,o - 32H 3 , - I6H4) + (1 -x)(-62 + 64H_ 2 ,o,o + 16H ,o^2 - 48H , ,o,o 

+ 560Hi +96Hi^ 2 - 16Hi, ) -80(1 +3x)H_ 2 ,o + 8(3 -5x)H 2 , + 32(3 -x)H_ 3 , 
+ 16(3 + 20x)H ,o,o - (3 + 539x)H - 24£ 3 (5 + 27x) - 2(9 - 191x)H ,o 
+ 16(9 + 13x)H 3 - 16^(9 + 28x)H + 4(77 + 85x)H 2 - 4^(77 + 93x) +96^ 3 H 
+ (1/2 (29 - 480£ 5 ) + 68C3 + I8C2- 32^ 3 + 288/5 C 2 2 ) 5(1 -x)) 

+ C A C F n f (p qq (x) (-753 1/81 - 64/3 ^ 3 + 48^ 2 - 16H_ 2 , - 3536/27 H + 32/3 H C 2 

- 608/9 H ,o - 40/3 H ,o,o - 1552/27 Hi + 16H^ 2 - 176/9 Hi, - 8/3 Hi, ,o 
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- 176/9 Hi,i + 16/3 Hi,i,o-16/3Hi )2 - 352/9 H 2 -8H 3 ) 

+ p qq (-jc) (-32/3 ^3 - 80/9 ^2 - 64/3 H_ 2 ,o + 32/3 H_^ 2 + 64/3 H_i _i, 

- 1 60/9 H_ i j0 - 32H_ ! ,o,o + 1 6/3 Ho - 8/3 H ^ 2 + 80/9 H ,o + 40/3 H 0)0 ,o) 
-8(1 + x)U 2 + 16/3 (1 -5x)(HiC 2 -Hi A o) + 16/3 (1 + 5x)(H_ 1 ^ 2 + 2H_ 1 ,_ 1 , 
-H_i A0 ) - 16/45 (3 - 6x 1 -497jc + 54x 2 )H + 8/ 15^ 2 (15 -55jc + 36jc 3 ) 

- 16/15 (55 + 2jr 2 + 35x-18x 3 )H_i, -4/9 (73 -269jc)Hi- 1/135 (7303 + 288A+ 1 

- 43013* + 2592x 2 ) +8/15 (55x - 36x 3 )H ,o - 160/3 - 32/3 H_ 2 , - (5516/27 

- 224/3 £ 3 +7216/27 £ 2 - 296/15 £f ) 8(1 -x)\ 

+ Cjn f (p m (x) (-67/3 + 208/3 ^ + 32H_ 2 , + 19H + 80/3 H C 2 + 12H , 

- 32/3 H OA o - 4Hi - 32/3 H^ 2 - 232/9 H lj0 - 80/3 Hi |0 ,o - 64/3 H u , - 32/3 Hi j2 

- 160/9 H 2 - 80/3 H 2 , - 16H 2J - 32H 3 ) + p qq (-x) (64/3 ^ + 160/9 C 2 

+ 128/3 H_ 2 ,o - 64/3 H_^ 2 - 128/3 H_i -i,o + 320/9 H_i )0 + 64H_i, ,o - 32/3 H 
+ 1 6/3 H ^ 2 - 1 60/9 H ,o - 80/3 H , ,o) + 8/3(1+*) (-2H ^ 2 + 7H , ,o + H 2 , 
+ 2H 3 ) - 32/3 ( 1 - x)Hi |0 - 32/3 ( 1 - 5x) (H1C2 - Hi j0 ,o) + 16( 1 + 3x)H 2 

- 32/3 (1 + 5jc)(H_i£2+2H_i _i )0 - H_i, ,o) - 16/15 £ 2 (15 -25x + 36x 3 ) 
+ 8/3 (19 - 42jc)Hi + 32/15 (55 + 2x~ 2 + 35jc - 18x 3 )H_i, + 4/45 (85 - 95x 

+ 432x 3 )H , + 8/45 (134 + 24.X+ 1 -269x + 216;c 2 ) - 1/45 (209+ 192* -1 + 4929* 

- 1728x 2 )H + 320/3 + 64/3 H_ 2 , - (239/6 + 400/3 ^3 - 146/3 ^ 2 

- 208/15 Cf) 5(1 -x)) 

+ C F nj (p qq (x) (470/81 - 32/9 £ 2 + 268/27 H + 16/3 H ,o + 1 16/27 Hi + 16/9 Hi, 
+ 16/9 H M + 32/9 H 2 ) -4/9(1 + 13x)Hi - 4/9 (3 + 23jc)H + 2/27 (29 - 295jc) 
+ (406/27 + 536/27^)5(1 -*)) . (B.l) 

The corresponding pure- singlet component reads 

C A C F n f {%{\ +x)(4H 2 ^ 2 - 7H 2A0 - 4H 3 , - H 4 ) + 8(1 -x)(4H_ 2 ^ 2 + 8H_ 2) _ 1)0 
+ 4H_ 2 ,o,o + H ,oC 2 ) + 32^(2 +x)H + 64/3 (3 - x~ l + 3x - x 2 )U^ h0 
+ 32/3^ 2 (3+^ 1 -3x-x 2 )Hi-32/3(3 + 2^ 1 + 3x + 2x 2 )H_i ! 2-16/3(3 + 8^ 1 
+ 3x + 8x 2 )H_! ,o,o + 16(3 -x)H_ 3 ,o + 16(4 - 3*)Ho,o,o,o + 8/3 (5 - 13x + 8x 2 )H 2 , 
+ 32/3 ^2(6 +jc~ 1 + 6jc + jc 2 )H_i - 8/3 £ 3 (11 - 32a:" 1 + S9x + 28jc 2 ) 
+ 4/5 CfC 11 +31jc) -4/3 £ 2 (15 - 16x _1 + 15x+ 16jc 2 )H + 4/3 (15 + 16a:" 1 Six 
+ 16x 2 )H 3 - 8(15 + 32x)H ,o,o - 4/9 £ 2 (21 - 368x _1 - 363* + 160x 2 ) 
+ 4/9 (21 + 104jc _1 +21x+ 160x 2 )H 2 - 8/3 (27 - 16x _1 + 33x + 24x 2 )H_ 2 , 
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- 4/3 (39 + 16.XT 1 - 39jc - 16x 2 )Hi A o - 4/3 (49 - 32jT 1 -x- 16jc 2 )H 1j0 

- 4/9 (57 - \21x~ 1 - A%x + UZx 2 )^ + 32/9 (84 + 59;T 1 + 48* + 23x 2 )H_i j0 

+ 4/9 (743 + 2&c _ 1 - 874x + 744x 2 )H ,o + 1 /27 (8 1 1 5 - 1 597 lx~ 1 + 8688x - 832* 2 ) 

- 2/27 ( 1 1373 + 820-xT 1 + 1 860k - 1396x 2 )H ) 

+ C 2 n f ( 1 6( 1 + x) ( -3H i C 2 - 6H_ ! ! , + 3H_ ! A o + 6H ,o^2 - 4H ,o,o,o - 2H 2 ^ 2 
+ 2H 2 , ,o - 2H 3 , - 6H 4 ) + 16(1 - jc) (-2H_ 2 C 2 - 4H_ 2 _ 1)0 + 2H_ 2 ,o, + 2H C 3 

- 3H1C2 + 3Hi A o) - 16/3 (3 - Ax- 1 - 8x 2 )H 2i o + 16/3 ^ 3 (3 + Ax- 1 + 36x) 

+ 16/3 (3 - lx + 24jc 2 )H ,o,o - 16/3 (9 - 4jT 1 + 12jc - 16jc 2 )H 3 + 16/3 ^ 2 (9 - Ax~ l 
+ 16x- 16x 2 )H + 32/3 (9 + 2x)H_ 2i0 + 16/5 Cf(l 1 + 6x) - 32/3 (12 - 5x~ l - 12x 
+ 5x 2 )Hi, - 8/9 (234 - 80x~ l - Tlx + 32x 2 )H 2 + 16/45 (585 + Ax~ z + 625x 

- 36x 3 )H_ 1;0 - 4/9 (657 - 149jT 1 - 450x - 58x 2 )Hi + 8/45 ^ 2 (1 170 - AOOx' 1 
+ 1 1 15jc + 160jc 2 - 72x 3 ) - 2/45 ( 1484 - 628jT 1 - 791 lx - 1588jc 2 )H 

- 4/45 (3015 + 1825x + 920x 2 - 144jc 3 )H 0j0 - 1/45 (6931 - 5184jc 1 - 7391* 
+ 5644x 2 )+64H_ 3i o) 

+ C F nj ( 1 6/3 ( 1 + x) ( -2^3 + 3H ,0,0 - H 2 ,o) - 64/9 (3 + x~ 1 + 3x + x 2 )H 1 ,0 

- 8/9 (3 + 4jT 1 - 3x - 4x 2 )H lj0 + 8/9 (23 - 12jc 1 + 95x - 20jc 2 )H , 

+ 4/27 (35 1 - 2 1 2x~ 1 + 1 1 Ax - 2A0x 2 ) H + 2/8 1 (5 1 39 - 2 1 55x _ 1 - 3420x + 436x 2 ) 
-64/9^ 2 (x- 1 +3x)) . (B.2) 

(2) (2) 

The off-diagonal contributions K qg and K gq can be written as 

4 2) w = 

C%n f (s Pqg (x) [-HC3H1 -8^2^,0+4^,0,0,0 -8C2Hi ) i+12H M) o+4H 1)1)0) o 
+ 14/3Hi, 2 +4Hi, 3 +61/3 H 2 ,i +4H 2 ,i ;0 +4H 2;2 +8H 3 ,i] + A Pqg (-x) [-16H_ 2 , 2 
-6^3 H_i - 16H_i -2,0 +8C2H_i _i +32H_i _i _i )0 -411 ,. ,,,,, + 12£ 2 H_i, 
-8H_i, ,o,o -20H_ 1]3 ] + 8/3^ 2 (-16jT 1 + 1 +llx+18x 2 )H_i 
+ 32(1 -4x-6x 2 )H+ 1;2 +8(l -2x) [-7^3 Hi -4C 2 H 1; o + 2H 1A0 , -4 C 2 H M 
+ 2Hi,i, 0) o+2Hi, 3 ] +8(1 +2x) [4H_ 2 , 2 -7£ 3 H_i -4H_i,_ 2 , + 8C 2 H_i,_i 
+ 8H_i_i _i,o -6H_i_i, 0) o -4C 2 H_i, +2H_i, OA o +2H_i, 3 ] -8/3(1 + 345* 
- 36x 2 ) H ,o,o + 64^ (2 + 9x) H - 8 (3 + 2x) H_ 3j0 - 12 (3 + 14* - Sx 2 ) H 4 
+ 16 (4 - 15 x) Ho,o )0) o - 8 (5 + 6x + 16x 2 ) H_ 2 , , + 8 ^ 2 (7 + 3Sx) H 2 + 4^ 2 (9 
+ 38jt-24x 2 )H ,o + 16(11 -30x + 8jc 2 )H_ 2 ,_i, + 16/3 (-16.XT 1 + 13 -37* 
-54x 2 )H_ li _i i0 -8/3(-8x _1 + 15 + 420x-436jc 2 )H 1 + 8^ 2 (15 -22* 
+ 24x 2 )H_ 2 +4/3 (32x^ + 17 -691x + 308x 2 )H 3 -4/3 £2 (17 -441x + 308x 2 )H 
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+ 8/3(17 - 184jc+116jc 2 )H 2 , -8/3(17 + 127 jc+150jc 2 )H_i A0 +8/3 ^(^jT 1 
+ 23 + 65jc-82jc 2 )Hi -4(25 + 90jc + 8jc 2 )H 2A0 +8/3 (37 -222jc + 210jc 2 )H m 

- 8/3 (-16x- 1 +46 + 125 x- 10jc 2 )H_ 2 ,o -4/3 (32x^+99 + 159 jc- 280jc 2 )Hi, , 
+ 8/9 (244JT 1 + 99 - 264* - 179 jc 2 ) H_ 1;0 + 4/3 ^ 3 (80 jc" 1 + 133 + 134jc - 286jc 2 ) 
+ 2/5 ^ 2 2 (149 + 174jc + 352jc 2 ) - 2/9 (- 192;T 1 + 163 + 2758jc - 2824 jc 2 ) Hi )0 

- 4/9 (-200 JT 1 + 375 + 2514jc - 1852jc 2 ) H 2 +4/9 ^ 2 (288JT 1 + 375 + 1986* 

- 1852jc 2 ) - 2/9 (372JT 1 + 1589 + 6362jc - 2308jc 2 ) H + 2/9 (56JT 1 + 2039 

- 4942 jc + 6340 x 2 ) H ,o + 1/27 (- 19297 jc" 1 + 25545 - 58161 jc + 46270 jc 2 ) 
-64((3x-x 2 )H 3 ,o -j^H-i-i^ +Hi,i |0 +H lj2 + 2H 2 ,i)) 

+ C}n f (l6 Pqe (x) [-4^ 2 H ,o -9C 3 Hi -8H1 _ 2 ,o - 10^Hi j0 + 4Hi, ,o,o -4C 2 H M 
-3Hi,i, +2Hi,i,o,o -3Hi )2 + 4Hi, 2 ,o +6Hi, 3 -4^ 2 H 2 + 4H 3 , +4H 4 ] 
+ I6p qg (-x) [-8C2H_ 2 -8H_ 2j _ lj0 +4H_ 2A0 -7C3H_i + 8£ 2 H_i_i 
+ 8H_i _i _i |0 +2H_i_i A0 -4H_i -1,2 -6^ 2 H_i, + 2H_i, 3 +4H ,o A o] 
+ 16(1 -2x) [-7C3H1 -4^2^,0+2^,0,0,0 -4^ 2 H U +2H 1)1A o + 2H 1;3 -2H 3 , 
-2H4] +16(1 +2x) [-7£ 3 H_i -4H_i,_2,o+8C2H-i,_i + 8H_i,_i,_i, 

- 6H_i -i A o — 4H_i,_i,2 — 4^2H_i,o +2H_i A o,o + 2H_i, 3 — 3Ho,o,o,o] 

- 16^3 (1 — 34jc + 4jc 2 ) H + 128 (1 +x 2 ) H 22 + 8 (3 - 2x - 12x 2 ) H 2 ,o - 32 (3 

- 2x + 8 jc 2 ) H_ 3 , - 32/5 ^ 2 (5 - 80jc - 15jc 2 - 24 jc 3 ) H - 8/5 ^ 2 (5 + 90* - 24 jc 2 ) 
-16(7 +2jc-12jc 2 )Hi A0 +8(13 -42jc + 34jc 2 )H m 

+ 16/15 ^ 2 (~^ 2 +15 + 70jc-180jc 2 + 36jc 3 )Hi - 16/3 J; 3 (15 - 175jc-36jc 2 

-36jc 3 ) + 32/15 (15 -55jc-45jc 2 -36jc 3 )H 3 + 4(19 -68jc + 52jc 2 )Hi, 

+ 32/15 (-x- 2 +30- 160jc-270jc 2 - 36jc 3 )H_! _ 1;0 - 32/15 (-jc~ 2 + 30+50jc 

-45jc 2 -36jc 3 )H_i A0 -16/5£ 2 (-*+ 2 +50 + 80jc + 90jc 2 + 36jc 3 )H_i 

+ 32/15 (jc~ 2 +90 +40jc + 36jc 3 )H_i, 2 -32/15 (150 - 185jc-270jc 2 -36jc 3 )H_ 2 ,o 

-4/15 (8x -1 + 161 + 1986jc-732jc 2 )H 2 -1/45 (64jc -1 + 217 +3528* 

- 13320 x 2 ) H +2/45 (-48JT 1 +249 + 2494jc- 1728jc 2 

- 2880 jc 3 ) H ,o +4/45 ^2 (48jT 1 +483 + 1426jc - 2196jc 2 + 1440 jc 3 ) 

-2/15 (-16JT 1 + 1849 + 756jc-2604jc 2 )Hi - 1/45 (-64 jT 1 +2572 + 13953 jc 

- 16596jc 2 ) -8/45 (-20x 2 - 12x 1 +3651 + 2266 jc- 1602 jc 2 -720jc 3 )H_i j0 
_ 64 (4 jc — jc 2 ) H 2)0 ,o - 16/15 (85jc- 90jc 2 +72jc 3 ) H oao - 192jc 2 H_ 2 ,o,o 
-32(2C2H_2-C2Ho,o-2Hi,i,o-2Hi, 2 ) + 128x[-2H_2,-i,o+C2H2]) 

+ C A C F n f (s Pqg (x) [23^ 3 Hi + 12Hi _ 2 ,o + 22£ 2 Hi, -4Hi A o,o + 16£ 2 Hi,i 
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+ 7/3Hi,i,o +29/3Hi, 2 -8Hi, 3 - 12H 2 ,i +4H 2 ,i, +4H 2 , 2 ] 
+ I6p qg (-x) [13 £ 2 H_ 2 + 8H_i,_ 2 , + 9H_i _i, 0|0 + 10H_i _i, 2 + 12C 2 H_i >0 
-6H_i, 0) o ; o -2H_i, 2) o -5H_i,3 -4^ 3 H ] +24(1 -2x) [7£ 3 Hi +4C 2 Hi, 
-2Hi, ,o,o+4C2Hi,i -2H M)0 ,o-2Hi )3 ] +16(1 + 2x) [3£ 2 H_ 2 +6H_i,_ 2 ,o 
+ 9H_i -1,0,0 +6H_i_i, 2 +6£ 2 H_i j o — 3H_i,o,o,o — 3H_i ;3 — 2H 3j o] 

- 32 (1 + 3x) Ho.0,0,0 + 64 ( 1 + 14x + 5x 2 ) H_ 2 ,-i,o + 8 C2 (3 - 2x - Ax 2 ) H ,o 
-8(3 +10x-4x 2 )H 4 -96^ 2 (5 +10x + 6x 2 )H_ 1 ,_ 1 +4/15 (5 + 1660x + 840x 2 
+ 288x 3 )H 3 -4/15 £ 2 (5 + 2960x + 840x 2 + 576x 3 )H -32(7 + 2x + 8x 2 )H_ 2;2 
-64(7 +14x + 8jc 2 )H_i,_i,_i, - 16(7 + 26x + 10jc 2 )H_ 2A o + 16(13 -6x 

+ 8jc 2 )H_ 3; o +8(17 +38jc+16jc 2 )H 2A o + 16^3 (25 + 50x + 29x 2 )H_! 

- 4/3 (-16X" 1 + 33 - 42x - 60x 2 ) H 2 , - 8/15 ^ 2 (-2x~ 2 + 60x 1 + 35 + 370x 
-465x 2 +72x 3 )H! +4/5^ 2 2 (59 + 314x- 120x 2 ) + 16/15 (2^ 2 + 60x- ! + 110 
+ 660x + 755x 2 + 72x 3 )H_i,_i, -4/3(131 -654x + 582x 2 )H u -4/3 £ 3 (175 
+ 934x+114x 2 + 144x 3 ) +4/15 (195 +550x+1800x 2 + 288x 3 )Ho,o, 

- 16/15 (2x _2 + 40x _1 +205 + 355x + 320x 2 + 72x 3 )H_i A0 

+ 4/3 (16x^+235 - 116x-72x 2 )H li0 , +4/15 (128X" 1 +261 + 3881x 

- 2142x 2 ) H 2 + 16/15 (265 - 325* - 755 x 2 - 72x 3 ) H_ 2 , + 2/15 (784x~ ! + 399 
+ 7156* - 8064 x 2 ) Hi - 16/15 (2x~ 2 + 20 x" 1 + 420 + 350* + 65x 2 + 72x 3 ) H_i, 2 
_ 1/45 (-848;T 1 + 695 - 25786x + 30216x 2 ) H - 4/45 & (408X" 1 + 783 

+ 253 1 x - 6426x 2 + 2016x 3 ) - 2/9 (-256x~ 1 + 83 1 - 3264* + 2596x 2 ) H 1;0 

- 4/45 (-24 jT 1 +927 + 1897jc + 4726jc 2 - 2016jc 3 ) H ,o + 8/15 ^ 2 (6x~ 2 + lOOx" 1 

+ 950 + 1360x + 885x 2 + 216jc 3 ) H_i - 1/45 (-4732x _1 + 3042 - 41977 x 

+ 46292 x 2 ) + 8/45 (-28x~ 2 - 12x _1 +4281 + 4556x-522x 2 - 1008x 3 )H_i j0 

-96^ 2 (4x-x 2 )H 2 -800^ 3 xH +128H 2;1 ) 

+ C A n 2 (s/3p qg (x) [-3C2H1 +6Hi,i +4Hi, 2 +2H 2 ,i +2H 3 ] + 8/3p qg (-x) [C 2 H_i 
+2H_i,_i,o+3H_i, 0)0 ] +8/3(1 +2x) [-2H_ 2)0 +H_i C2 + 2H_i,_i, -H_i )0)0 ] 
+ 8/3 ^3 (1 -46x+10x 2 ) +8/3^(1 -2x)Hi -16/3(1 + 4x)H 2)0 

- 16/3 ^2(1 +2x 2 )H +8(2 -3x+x 2 )H! +8/9 (-4x~ ! +5 +38x + 39x 2 )H , 
+ 16/9 (-2x~ ! +6 -15x+llx 2 )H 2 - 16/9 £ 2 (2x _1 +6 +3x+llx 2 ) 

+ 8/9(-4x" 1 +8 -46x + 53x 2 )H li0 -16/9(4x^ + 15 +18x + 4x 2 )H_ li0 
-4/27(114x _1 + 69 +90x-278x 2 )H +2/81 (-1328X" 1 +3105 -3114x 
+ 3524x 2 ) + 16/3 (x 2 H 1 ,o, - H ,o,o - 2H M ) ) 
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+ C F n 2 (s/3p qg (x) [-2H 1)0)0 -3H U +2Hi,i )0 -2H 1;2 -3H 2 + 3H 2;0 -2H 3 ] 
+ 32/3p qg (-x) [-C2H-1 -2H_i,-i,o+H_ 1AO ] +(1 -2x) [-32H ,o A o 
-16/3H^ 2 ] +16/3^2(1 -14x + 2x 2 )H -4/9 (8jc _1 +3 + 48x-44x 2 )Hi, 

- 16 (3 - 1 1 x - 2jc 2 ) H ,o,o + 32/3 (5 - 6x + Ax 2 ) H_ 2;0 + 16/3 ^3 (7 - 6x) 
-8/45^(40x^ + 15 +10x-90x 2 -72x 3 ) + 16/45 (jc" 2 -20x" 1 +225 +40jc 

- 155x 2 +36x 3 )H_i, -8/45 (40x^ + 885 + 470x + 460x 2 + 72x 3 )H ,o 

- 4/135 (712x^+5304 - 6771 x- 9078x 2 ) H - 1 /405 (9796x~ ! +7731 

- 15 1506* + 105494 x 2 ) - 8/3 (2x _1 + 6x - 1 x 2 ) Hi + 32/3 (Hi,i +H 2 ) 
+ 16/3x 2 [2^ 2 H_j +4H_i_i, -2H_i A0 -2H 1A0 -H 2 , ]) 

+ nj (%/9p qg (x) [-H 0j0 -Hi )0 ] +8/27 (2x ^ 3 +6jc)H +4/81 (13x~ ! -45 

+ 81x-49x 2 )) (B.3) 

nd 

K$(x) = 

C 2 A C F (8pgq(*) [-8H1 _2,o + 2Hi A o,o +8£ 2 Hi,i +13/3H M ,o -4H UA0 + 13/3 Hi, 2 
+ 8H 1A0 +4H 1;3 ] + 16/? gq (-x) [-4H_ 2 . 2 +8H_ l! _ 2 , - 14£ 2 H_i _! - 8H_i _i _i, 
+ 13H-i,-i A o + lOH-1,-1,2 +7 C 2 H_ 1)0 -5H_i,o A o -2H_i, 2 ,o -5H_i >3 ] 
+ 24(2 -x)(-Hi,i, -Hi, 2 )-104C3(-2^ 1 +2 -x)Hi - 168^3 (2x _1 +2 + x)H_i 
+ 8 (32x _1 + 2 + 25x) H_ 2A0 + -4 & (2 + 1 1 x) H ,o + 4 (2 + 13 x) H 4 - 16 ^3 (-8 x~ ! 
+ 4 - 9x) H + 32 (2x -1 + 4 + 3x) H 3 , - 8 ^ 2 (-8X" 1 + 6 - 5x) H 2 - 16 (8x _1 + 6 
+ 5x) - 8 (16 - 15 x) Ho.0,0,0 + 4 (18 +x) H 2 , , - 8 C 2 (24x~ ! + 30 + 13*) H 2 

-8(34-x)H_ 3 , -4/3 ^(-6x^+56 -7x + 16x 2 )Hi - 16/3 (-20X" 1 +65 -Ix 
+ 12jc 2 )H 2 , -4(-92x~ ! +75 -32x)Hi A o + 16/3 (75jT ! + 106 + 50x + 8x 2 )H_i ;2 
+ 8/3 (22JT 1 + 110 + 19*- 24x 2 ) H 0A0 - 8/3 (58X" 1 + 144 + 51x+ 16x 2 ) H_i _i )0 
+ 4/3 ^(-44x^ + 187 +145x + 32x 2 )H -4/3 (44x^ + 187 +63x + 32x 2 )H 3 

- 4/9 (- 17X" 1 + 242 - 121 x) H M - 2/5 £ 2 2 (384x" J + 242 + 247 x) +8/3 (-44 x" 1 
+ 353 +41x + 48x 2 ) H_ 2 , + 8/3 £3 (-347 x" 1 +444 - 24x+ 16x 2 ) - 4/3 £ 2 (358X" 1 
+ 568 + 25 1 x + 48x 2 ) H 1 + 8/9 £ 2 (-269 jT 1 + 704 + 24x + 96x 2 ) + 4/3 (632x _1 

+ 724+215x + 64x 2 )H_i A0 - 4/9 (671 x" 1 + 1408 +205x+ 192x 2 )H 2 

- 4/9 (1209 x- 1 + 1484 + 157x + 472x 2 ) H_i ;0 +4/27 (-5419X" 1 + 1687 +709x 
+ 1446x 2 ) Hi - 2/9 (-1946X" 1 +2190 - 1269x - 400x 2 ) Hi )0 - 2/9 (-1892X" 1 

+ 10782 - 1667x+ 1632x 2 ) H , + 1/27 (28529X" 1 + 27762 + 10188x+ 15400 x 2 ) H 

- 1/54 (-223601 x _1 + 172270 -72332x-7424x 2 ) + 64H_ 2)2 ) 
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+ C A C 2 (Sp gq (x) [4Hi _2, - 14^2 Hi, + 10H 1;0) o,o - 16£ 2 H U + ll/3Hi,i )0 
+ 16H U! o,o + H/3 Hi, 2 +8H lj3 + 10H 2A0 +5H 2 j +4H 2i1j0 + 4H 2;2 ] 
+ I6p gq (-x) [-8H_i _ 2) o + 14^ 2 H_i,_i + I6H-1 -1-1,0 -7H_i -i A o -6H_i _i i2 
-2C2H_i ) o-H_i ) o,o ) o-H_i )3 ] +16(2 -x)(-3H 2 ,i -H 3)0 ) -72^(2 +x)H , 
+ 264^ 3 (-2x- 1 + 2 -x)Hi +184£ 3 (2x~ 1 + 2+x)H_i +16(4+3x)H OA o,o 

- 32 £ 3 (3x~ 1 + 5 + 9x) H + 24 ^ 2 (4x~ 1 + 6 + 5 x) H 2 - 8 ( 1 2x~ 1 + 6 + 1 1 x) H_ 2 A0 
+ 16 (12x _1 + 10 + 15*) H-2,-1,0 + 8 1,2 (-20x~ ! + 18 - 19 x) H 2 - 32/3 (1 Ijc 1 + 18 
+ 9x + 2x 2 )H_i ;2 +8(18 + 7x)H 4 -8(6x^+44 + 17x + 8x 2 ) H_ 2;0 

+ 16/3 ^(-18x^+56 -49x + 6x 2 ) + 8(30x~ 1 +60 + 13x + 8x 2 )H_i,_i,o 

- 4/5 J; 2 2 (-44x^+72 + 27x)-8/3(-45x~ 1 +94 -17x + 8x 2 )Hi A0 

-2/3 (-210x 1 + 116 -23jc)Hi j0 -4(134 -67x + 16x 2 )H OA o -4/3(-96* _1 + 142 

- 95 x) H 2;0 + 4/ 1 5 ( 1 6 2 - 1 35 x~ 1 + 1 80 + 205 x + 36x 3 ) H i >0 - 4/3 ( 1 34x~ 1 + 204 
+ 63 x + 32x 2 ) H-1,0,0 + 4/3 £ 2 (- 134x _ 1 + 224 - 61 x + 24 x 2 ) Hi + 4/3 £ 2 ( 178x~ 1 

+ 324 + 1 1 1 x + 40 jc 2 ) H i - 8/9 (- 1043x~ 1 + 420 - 201 x + 65 x 2 ) Hi 

+ 4/3 & (-256X- 1 + 556 - 299x + 48x 2 ) H - 4/3 (-220JT 1 + 556 - 197x + 48x 2 ) H 3 

-4/9 (-1497x^+607 -683x-208x 2 )H 2 +4/45 ^ 2 (-7890X" 1 +3035 -2800x 

- 1040jc 2 + 108 jc 3 ) + 2/45 (4450 + 2785* + 2080x 2 -216x 3 )H ,o - 2/45 (-7269 x" 1 
+ 10022 -6078x + 1756x 2 ) H - 1/90 (-52479X" 1 + 133506 - 13606x-3656x 2 ) 

+ 68^^1,1 + 16x(-H_ 3;0 + 6H 2> o,o) - 64H_ 2 , 2 ) 

+ Cj (Sp m (x) [14^ 3 Hi +8C 2 Hi j0 -4Hi A0j0 +8C 2 Hi ; i -4H UA0 -4Hi >3 +3H 2 , 
+ 4H 2) i i0 +4H 2 , 2 +4H 3 , +8H 3 ,i] + I6p m (-x) [l^H-i +4H_ 1; _ 2! o -8£ 2 H_i _i 
-8H_i_i _i,o +6H_i _i,o,o +4H_i _i ;2 + 4C 2 H_i j0 -2H_i, ,o,o -2H_i, 3 ] 
+ 64(1 +x)H_i 5 2 + 16(2-x)[+2H 2)0 +H 3 -H 3 , -3H 4 ] - 16(2 +x)H OA o,o 

- 16& (-6*" 1 + 2 - 1 1 x) H + 16 (2 - 9x) H OA o + 16 (2 - 5x) H 2A0 + 32^ 3 Ox' 1 
+ 3 +8x)+32C 3 (-2x" 1 +4+5x)H + 24(2x~ ! +4 + 3x)H_i A0 + 32(3x^+4 

+ 5x)H_ 2 ,o +8(4 -x)Hi A0 +32(4 + x)H_ 3 ,o - 16^ 2 (6 -5x)H 2 + 16^ 2 (6 -x)H ,o 

- 16(6x _1 + 8 +5x)H_i_i j0 -8^ 2 (8 -5x)Hi -32(4x~ 1 + 10 +7x)H_ 2 -i,o 
+ 16 (4X- 1 + 14 +7x) H_ 2A o - 8/5 & (-16X" 1 + 16 +x) - 8 & (6*" 1 + 16 

+ 13x)H_i -4(-14x _1 + 18 +x)Hi, -4(-10x" 1 + 18 +25x)H 2 - 16^ 2 (4x" 1 + 18 
+ 7x)H_ 2 -2(+27x _1 +42 - 16x)Hi + 8/15 (-16x~ 2 + 285x _1 + 180 - 85x 

- 36x 3 ) H i ;0 + 4/ 1 5 £ 2 (420 x~ 1 + 270 + 205 x - 72 x 3 ) - 4/ 1 5 (360 + 85 x - 72 x 3 ) H 0;0 

- 2/15 (716* -1 + 2132 +597x- 144x 2 )H +2/15 (-2734x^ + 2716 +369x+ 144 x 2 ) 
+ 16x- 1 (2Hij +3Hi,i )0 + 3H 1;2 + 6H 2 ,i ) + 128H_ 2 , 2 ) 
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+ C 2 n / (l6/3p gq (x) [7^ 2 H -2^ 2 H 1 + 2Hi, ,o -H M ,o -H 1)2 -3H 2 ,o -3H 2 ,i -7H 3 ] 
+ 1 6 p gq ( -x) [- 1,2 H i - 2 H i _ i , + H i ,o,o ] - 32 ( 1 + x) H i , + 8/3 (2 - x) [-2 ^ H 
+ 12H OAO ,o +H 2 , +2H 3 ] +4/3(-18x- ! +4 - x)H 1;0 +8(8 + 3x)H OA o 
+ 4/3 (-91jc 1 +30 -7x)Hi - 16/3 £ 3 (-26*" 1 +38 - 13 x) +8/9 (-161-xT 1 + 110 
-67jc)H 2 -8/9^ 2 (-161jc" 1 + 110-31x)+4/9(104^ 1 +506+47x)H ,o 
+ 4/27 (55;T 1 + 1401 - 198x- 112x 2 )H + 1/27 (-5779*" 1 + 13788 - 8223* 
+ 304x 2 ) -24jc _1 Hi,i -64H_ 2 ,o) 

+ C A C F n f (4/9p gq (x) [-12^ 2 H +42^ 2 Hi -90Hi |0 ,o - 17Hi,i -24Hi,i |0 -24Hi, 2 ] 
+ 8/3p gq (-x)K 2 H_! +6H_i_i,o+7H_i i o,o+4H_i ) 2+2H 3 ] + 32/3 (-3*- 1 + 1 

- 2x) H 2 , + 32/3 (-2x~ l + 1 - 2x) H_ 2 ,o + 12 (2 -x) Hi,i + 8/9 (-45.XT 1 + 4 + 17x 
+ 8x 2 )H_ 1; o + 16C 3 (-4x" 1 +4 -3jc) -32/3 (jc -1 +5 +4x)H ,o,o +8/9(-17x~ 1 

+ 160 - 26x+ 16x 2 ) H 2 - 8/9 £ 2 (28x~ ! + 160 - 43x+ 16x 2 ) + 8/9 (-212*" 1 + 232 

- 101x + 4x 2 )H lj0 +8/9 (-229x^+368 -201x + 44x 2 )H , +4/27 (26x^+418 
-407x- 192x 2 )H! +4/27 (-2579x _1 +4572 - 1050x+ 112x 2 )H 

+ 2/81 (-21551X" 1 + 11598 - 12585x-6364x 2 ) -32/3xH 3 ) 

+ C F nj (8/27 p gq (x) [39 H ,o + 29 H i + 1 5 H i , + 6 Hi , i + 1 2 H 2 ] + 32/9 ^ 2 ( -2 x~ 1 + 2 
-x)-4/9(-119x" 1 + 106 -52x)H -2/27(-1003x" 1 + 814 -422x)) . (B.4) 
Finally the second diagonal NNLO entry is given by 

K$(x) = 

Cj (p gg (x) [18974/27 -440^ 3 - 3008/9 £ 2 -24/5£f -64H_ 3 , +96^ 2 H_ 2 
+ 64H_ 2; _ l!0 -352/3 H_ 2 ,o -64H_ 2 , 0) o -64H_ 2;2 + 11185/27 H - 112^ 3 H 

- 880/3^2 Ho +3008/9 H ,o - 128^ 2 H ,o +440/3 H OA o +64H 0) o A o +9782/27Hi 
-96H^ 3 - 176/3£ 2 Hi -64Hi_ 2 ,o +1360/3 Hi j0 -96£ 2 Hi |0 +616/3H 1A0 

+ 128H 1AO ,o +968/9Hi,i + 176Hi,i |0 + 128H MA0 + 176H lj2 + 128H 1A0 

+ 128Hi, 3 + 1360/3 H 2 - 64^ 2 H 2 + 176H 2 , + 160H 2A0 + 176H 2il + 128H 2 , li0 

+ 128H 2;2 +968/3 H 3 +160H 3 , + 128H 3 ,i +128H 4 ] + Sp gg (-x) [-77/3£ 3 

- 134/9 £ 2 + 11 Cl - 8H_ 3 ,o + 32^ 2 H_ 2 + 16H_ 2 ,_ l!0 - 22H_ 2 ,o - 36H_ 2 , 0| o 
-24H_ 2j2 +24C 3 H_i +88/3C2H_i+16H_i_2,o-32C2H_i_i + 88/3H_i,-i i0 
+ 48H_i _i,o,o +32H_i _i,2 - 268/9 H_i, +36^2 H_i, - 110/3H_i A o 

- 32H_i, OA o - 44/3 H u - 8H_i A o - 32H_i, 3 - 10£ 3 H - 1 1 £ 2 H + 134/9 H ,o 

- 16^ 2 H ,o + HHo.0,0 +8H ,o A o +4^H 2 +22/3H 3 +4H 3 , + 16H 4 ] 

+ 32(1 +x) [-2C2H2 +5H 2A0 +8H 3 , + 13H 4 ] +64(1 -x) K 2 H_ 2 -2H_ 2 ,-i,o 
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-3H_ 2) o,o] — 16 (— IIjc 1 +3 +3jc-11jc 2 )H_i_i ;0 -8^ 2 (11x" 1 + 3 -3jc 

- 1 1 x 2 ) Hj - 128 ^ 3 (3 + x) H + 8 (12JT 1 + 5 - 5x - \2x 2 ) H 1>0 - 64 (5 - x) H 3j0 

- 32^2 (13 + 1 1 x) H ,o + 16/3(19 +117x- 132x 2 ) H ,o,o - 16/3 (-22JT 1 +21 
-51x + 66x 2 )H 2 , -8(-33x _1 +25 -25x + 33x 2 )Hi A o + 16 (IIjc 1 +25 +25* 
+ 1 lx 2 ) H_!, 2 - 16/5 C| (33 + 43jc) - 8 & (1 1 *~ l + 53 + 53x+ 1 lx 2 ) Hi 

+ 32/3 (44jT 1 + 57 + 51 x + 44 x 2 ) H_i A o + 8 ^ 3 (-55 x -1 + 59 + 22* + 66x 2 ) 
+ 16/3 (-22x^+62 + 19x + 77x 2 )H_ 2 ,o + 8/3 £ 2 (-33JT 1 +75 -7x + 231x 2 )H 
-8/3(llx _1 + 75 -45jc + 231jc 2 )H 3 + 8/9 (-134jT 1 + 137 -316x-229jc 2 )H 2 
-8/9^ 2 (631^ 1 + 137 +461jc-229jc 2 ) -8/3 (255.XT 1 +259 + 259jc + 255x 2 )H_ 1 , 

- 8/9 (-242JT 1 + 1593 - 642x + 2217x 2 ) H ,o +2/27 (-3774 jtT 1 + 8194 - 6797* 
+ 5105x 2 )H! + 1/27 (23364JT 1 +65830 +21391 x + 53565x 2 ) H 

- 2/8 1 (- 1295 17 jT 1 + 83236 - 96986* + 107847 x 2 ) + 448xH ,o A o - 1 /486 (592399 
+ 38880^ 5 )8(l-x) +(1988/3^ + 11185/27^2 -16^2^3 -88£|)8(1-jc)) 

+ C\ n f (l p gg (x) [- 10429/8 1 + 24 ^3 + 48 & + 32/3 H_ 2 ,o - 228 1 /27 H + 80/3 & H 
-48H ,o -40/3 Ho,o,o- 1882/27Hi +16/3 ^ 2 Hi - 512/9 Hi, -56/3 Hi A o 

- 176/9 Hi,i -16H Ui0 -16H 1)2 -512/9 H 2 - 16H 2 , - 16H 2 j -88/3H 3 ] 
+ 16/9p gg (-x) [21^3 +10^ 2 + 18H_ 2 ,o -24^9 H_i -24H_i _i |0 +20H_i |0 

+ 30H_i )0 ,o + 12H_i, 2 +9 ^2 H - 101% -9H OA o -6H3] +8 (1 +x) [-6^ 2 H 2 
-4H 2 ,o +3H 2 ,o,o +3H 4 ] +48(1 -x) K2H-2 -2H_ 2 ,-i,o +H_ 2)0)0 ] 
+ 32/3 (jT 1 + 3 + 3 x + x 2 ) H_i )2 + 2/27 (374JT 1 + 3 - 807* - 1 100x 2 ) Hj 
-8^2(3 +5x)H ,o + 16^3 (-jT 1 +5 +5x) + 16(5 -jc)H_ 3 ,o - 16(-2jc _1 +7 +7 x 
-2x 2 )H_i_i )0 -8^ 2 (2x" 1 +7 -7x-2x 2 )Hi +40/3 (-2x~ l + 17 + 11 x 
-2jc 2 )H_i ;0 -8/3 (20+51x)H ,o,o +4/3 (4x _1 +21 -21jc-4x 2 )Hi A0 
+ 12/5 C 2 (23 + 13*) + 16/9 (-35jT 1 +27 -27jc + 35jc 2 )H!, + 8/3 {-Ax' 1 +27 
+ 27* - 4x 2 ) H_i A o - 4/9 (88;T 1 + 32 + 761 x - 396x 2 ) H ,o - 8/3 ^ 2 (-2JT 1 + 33 
+ 33x-2x 2 )H_ 1 + 4/3^ 2 (-4jT 1 +39 +71x+ 12x 2 )H -4/3 {-Ax' 1 +39 +81x 
+ 12x 2 )H 3 +8/3 (47 -5x-12x 2 )H_ 2i0 +4/9(-52^ 1 + 145 +193x+140jc 2 )H 2 
-4/9£2(8jc -1 + 145 -317x+ 140 x 2 ) -2/27 (2008-xT 1 +586 + 3196*- 1757x 2 )H 
-2/81 (7616X" 1 + 13155 -5298x + 2473x 2 ) + 16x[-6^ 3 H +H ,o A o] 
+ (89027/162 - 80/3 ^3 - 4562/27 £ 2 + 16^1 ) 8 (1 -*)) 

+ Cpn f (32(1 +x) £ 2 H , -Ho.0,0,0 " C2H2 +H 2A0 -H 3 , -H4] 
+ 32 (1 -x) K2 H_2 - 2 H 2 , 1 , + H_ 2 ,o,o ] - 1 6 £ 3 (4x~ 1 + 1 + Ax) 
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- 32^ (1 + 3 x) H + (-4jT 1 + 3 + 3x - Ax 2 ) [-16/3 ^ 2 H_i - 32/3 H_i _i )0 
+ 16/3 H_i, 0|0 ] + (4X- 1 + 3 - 3x - Ax 2 ) [-16/3 £ 2 Hi + 16/3 H 1A o] 

+ 16/5 C| (3 - 2x) + 32/3 (3 -2x- Ax 2 ) H_ 2 ,o - 16/3 (3 + 2x) H ,o,o - 16/3 (Ax~ l 

+ 6 + 3x) H 2 , + 16/3 ^ 2 (Ax' 1 +9- 10*) H - 16/3 (4jT 1 +9 - 6x) H 3 

+ 8/9(-44x" 1 +27 -27x + 44x 2 )Hi, + 64/45 (-jT 2 + 105 + 95x + 9jc 3 )H_i, 

-4/45 (135 + 1835jc-440jc 2 + 144jc 3 )H , +4/27 (-691x _1 +675 - 5Ax 

+ 70x 2 )Hi - 2/135 (1884x _1 + 1368 +6933*+ 164 x 2 )H + -1/135 (7632JT 1 

+ 4807 - 3427 x - 901 2x 2 ) - 8/9 (80x _1 + 27 x) H 2 + 8/45 & (400 jT 1 + 895x + 72x 3 ) 

+ 64H_3,o+5(l-*)) 

+ C A C F n f (&p gg (x) [-22/3 + 4£ 3 -H -Hj ] + 16 (1 +x) [6^ 2 H ,o +7^ 2 H 2 -7H 2 , , 
-2H 3j o-6H 4 ] + 16(1-*) [7H_ 2 C2+14H_2-i,o-3H_2,o l0 ] + 128^(1 +2x)H 
+ 64(2 -3jc)Ho,o,o,o -64/3 (x^ 1 +3 +3x+x 2 )H_ li2 -4/3 (32JT 1 + 11 -59x 
+ 16x 2 ) Hi, - 16/3 (4jT 1 + 15-18*- 4jc 2 ) H 3 + 16/3 £ 2 (15 - 3x- Ax 2 ) H 

- 8/3 (8jT 1 + 23 - 7jc) H 2;0 - 8/3 ^ 3 (-20jT 1 + 25 +43 x + Sx 2 ) - 16/3 (Ax~ l + 30 

- 15*- 12x 2 )H_ 2 , -8/3 (20x _1 +33 -33x-20x 2 )Hi, ,o -8/3 (43 + 31x)H ,o,o 
+ 16/3(-16j»T 1 +45 +45x-16x 2 )H_ 1 ,_ 1 ,o -8/3 (Ax- l +A5 + 45x + 4x 2 ) H_i )0)0 
+ 8/3£ 2 (16jT 1 +45 -45jc-16jc 2 )Hi + 8/3 ^ 2 (-8*" 1 +69 + 69x-8jc 2 )H_i 

+ 16/9 ^ 2 (87 -216jc+14jc 2 ) - 16/9 (87 - 3x+ Ux 2 ) H 2 - 16/3 (95 + 71* 

- 2Ax 2 ) H_i,o + 4/9 (- 1 16X" 1 + 335 + 1223* + 360 x 2 ) H ,o - 4/27 (-241 * _1 
+ 2196 -2223x + 268x 2 )H 1 - 2/27 (1284JT 1 + 2463 - 2286* + 8464x 2 )H 

- 1/81 (24445 jiT 1 +5079 +31074jc- 52282x 2 ) — 8 (7 C2 + 12H -3,o) + (465/2 
-176^3 -8^2)5(1-*)) 

+ C A nj (4/81 p gg (x) [571 - 144^ 2 +417H + 144H , +318H! + 144H 1;0 + 72H M 
+ 144 H 2 ] +32/9(1 + x)H 0;0 +32/9(2 -x + x 2 ) [-£3 +H 2 ] +2/27 (52x^ + 225 
-81jc + 80x 2 )Hi + 2/27 (192X" 1 +404 +5x+146x 2 )H +2/81 (1540x~ 1 +891 
+ 303x-202x 2 ) + (-1139/18 - 16^ 3 +556/27^)8(1 -x)\ 

+ C F n 2 (Ap gg (x) +16/3 (1+*) [2^3 +5H ,o,o +H 2 ,o] + 8/9 (Ax~ l +3 - 3x 

- 4x 2 )H 1:0 + 8/9(20x _1 + 47 + 5x-4x 2 )H , +4/9 (IOOjT 1 + 155 + 2x) H 
+ 2/27 (1079x 1 +309 - 1314x+ 178x 2 ) + (-43 + 32 ^ 3 ) 8( 1 - x) \ 

+ nj ^8/81 p gg (x) [-10 - 9Hq - 6H1 ] + 8/81 (2 -x + x 2 ) [-10 - 9H - 6H1 ] 

+ 8(25/243 -1/9 ^ 2 ) 8(1 -x)) . (B.5) 
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Appendix C: The gluon coefficient function for F 2 at large x 



(n) 

In this final appendix we provide the large-x coefficients of the gluon coefficient function c\ l(x) 
for the photon-exchange structure function F 2 , which were not written down in Ref . [9] for brevity. 
As their counterparts for c^ 1 ' Jx) in Eqs. (13.221) - (13.291) . these coefficients contribute to the large-* 

(n) in) 

behaviour of the off-diagonal physical kernels and K± 2 presented above. 



The second-order coefficients are (recall L\ = In (1 — x)) 

c (2) 

C 2,g 
.(2) 



C A n f + — C F n f , (C.l) 



C 2,g 



c (2) 

C 2,g 



The corresponding three-loop results read 



3 - 2 A f T F f 

2 = -4C A n f -9C F n f , (C.2) 
1 

(14 -S^ 2 )C A n f - (2 + 16&)C F n f . (C.3) 



c 2,g 



c (3) 

C 2,g 

.(3) 
•2,g 

c (3) 

C 2,g 

c (3) 

C 2,g 



? 10 9 

■Cln f + — C}n f , (C.4) 



293 2 206 83 2 17 2 7 2 

L 4 = -54 C A n f-^j- C FC A n f - — C F n f + — C F n f + — C A n f , (C.5) 

/3056 136 N 2 /3442 136 \ 

127 376 x 2 496 2 152 

1,2 UHf C F n f 

9 9 S / F f 81 7 81 



+ (— ^^)C}n f -—C F nj -—C A nj , (C.6) 



12043 232 52 \ 

13789 652 . 176 . \ « / 1096 40 . \ „ . 

C2 C3 Ctrif + t 2 )C A n; 

81 9 b 3 b V A f \ 81 9 b / y 

205 .„ \ / 1565 16 
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142 C 2 - 72 C 2 » r + ( - ^ C2) C F , (C.7) 



c (3) 

C 2,g 



/ 32968 2380,, 1250 . 64 „ 2 \ „ „ 

= 1 £2 C3 U)C F C A n f 

L\ V 243 9 S 9 s 5 s2 / 7 

/ 69526 6052 . 952 . cn . 2 \ - 

+ C2 C3 +50Co )Cin f 

\ 243 27 S 9 S S2 / A r 

/254 46 . 632 . 2 \ 2 

+ (— -y^2 -456^3+ — Ci • J Q «/ 

/ 10798 496 . 16 . \ 2 /3076 328 . 64 . \ 2 ^ on 

Note the close similarity between Eqs. (13.221) , (13.251 ), (IC.ll ) and (IC.4I ) for the leading logarithms at 
NNLO and N 3 LO which may point to yet another general structure. 
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